Terahertz spectroscopy of bacteriorhodopsin and rhodopsin: similarities and differences

Selectivity of Protein Interactions Stimulated by Terahertz Signals

It has been established that Terahertz (THz) band signals can interact with biomolecules through resonant modes. Specifically, of interest here, protein activation. Our research goal is to show how directing the mechanical signaling inside protein molecules using THz signals can control changes in their structure and activate associated biochemical and biomechanical events. To establish that, we formulate a selectivity metric that quantifies the system performance and captures the capability of the nanoantenna to induce a conformational change in the desired protein molecule/population. The metric provides a score between -1 and 1 that indicates the degree of control we have over the system to achieve targeted protein interactions. To develop the selectivity measure, we first use the Langevin stochastic equation driven by an external force to model the protein behavior. We then determine the probability of protein folding by computing the steady-state energy of the driven protein and then generalize our model to account for protein populations. Our numerical analysis results indicate that a maximum selectivity score is attained when only the targeted population experiences a folding behavior due to the impinging THz signal. From the achieved selectivity values, we conclude that the system response not only depends on the resonant frequency but also on the system controlling parameters namely, the nanoantenna force, the damping constant, and the abundance of each protein population. Based on the selectivity metric, the nanoantenna must be tuned to a frequency that is not necessarily the resonant frequency of the protein. The presented work sheds light on the potential associated with the electromagnetic-based control of protein networks, which could lead to a plethora of applications in the medical field ranging from bio-sensing to targeted therapy.

Epithelial Sodium Channel Inhibition by Amiloride Addressed with THz Spectroscopy and Molecular Modeling

THz spectroscopy is important for the study of ion channels because it directly addresses the low frequency collective motions relevant for their function. Here we used THz spectroscopy to investigate the inhibition of the epithelial sodium channel (ENaC) by its specific blocker, amiloride. Experiments were performed on A6 cells' suspensions, which are cells overexpressing ENaC derived from Xenopus laevis kidney. THz spectra were investigated with or without amiloride. When ENaC was inhibited by amiloride, a substantial increase in THz absorption was noticed. Molecular modeling methods were used to explain the observed spectroscopic differences. THz spectra were simulated using the structural models of ENaC and ENaC-amiloride complexes built here. The agreement between the experiment and the simulations allowed us to validate the structural models and to describe the amiloride dynamics inside the channel pore. The amiloride binding site validated using THz spectroscopy agrees with previous mutagenesis studies. Altogether, our results show that THz spectroscopy can be successfully used to discriminate between native and inhibited ENaC channels and to characterize the dynamics of channels in the presence of their specific antagonist.

Hydration Dynamics in Biological Membranes: Emerging Applications of Terahertz Spectroscopy

Water drives the spontaneous self-assembly of lipids and proteins into quasi two-dimensional biological membranes that act as catalytic scaffolds for numerous processes central to life. However, the functional relevance of hydration in membrane biology is only beginning to be addressed, predominantly because of challenges associated with direct measurements of hydration microstructure and dynamics in a biological milieu. Our recent work on the novel interplay of membrane electrostatics and crowding in shaping membrane hydration dynamics utilizing terahertz (THz) spectroscopy represents an important step in this context. In this Perspective, we provide a glimpse into the ever-broadening functional landscape of hydration dynamics in biological membranes in the backdrop of the unique physical chemistry of water molecules. We further highlight the immense (and largely untapped) potential of the THz toolbox in addressing contemporary problems in membrane biology, while emphasizing the adaptability of the analytical framework reported recently by us to such studies.

Broadband terahertz signatures and vibrations of dopamine

Dopamine (DA) is an essential neurotransmitter and hormone of the nervous system, its structural and conformational properties play critical roles in biological functions and signal transmission processes. Although this neuroactive molecule has been studied extensively, the low-frequency vibration features that are closely related to the conformation and molecular interactions in the terahertz (THz) band still remain unclear. In this study, a broadband THz time-domain spectroscopy (THz-TDS) system in the frequency band of 0.5-18 THz was used to characterize the unique THz fingerprint of DA. In addition, density functional theory (DFT) calculations were performed to analyze the vibrational properties of DA. The results suggest that each THz resonant absorption peak of DA corresponds to specific vibrational modes, and the collective vibration also exists in the broadband THz range. Moreover, the interactions between the DA ligand and the D2 and D3 receptors were investigated by docking, and the simulated THz spectra were obtained. The results indicate the dominant role of hydrogen bonding interactions and the specificity of molecular conformation. This work may help to understand the resonance coupling between THz electromagnetic waves and neurotransmitters.

Quantitative analysis of bisphenol analogue mixtures by terahertz spectroscopy using machine learning method

Quantitative analysis of complex mixtures is a great challenge for spectral analysis. Bisphenol A (BpA) is a chemical predominantly used in manufacturing and is being replaced by other analogs due to its potential toxicity. Reliability methods is hence crucial for identification and quantification of bisphenol mixtures. In this study we present an attractive strategy for composition determination of BpA incorporated in its analogue mixtures. Terahertz spectra of four bisphenol components are analyzed using machine learning method (SVR) to learn the underlying model of the frequency against the target concentration of BpA in mixtures. The learned mode predicts the concentrations of the unknown samples with decision coefficient R² = 0.98. Absorption spectra for bisphenols mixtures were successfully reconstructed by a hold-out validation scheme. The results indicate the terahertz spectroscopy in combination with SVR is robust and accurate in mixture quantitative analysis and should play a significant role for industrial applications in the future.

Information Rates of Controlled Protein Interactions Using Terahertz Communication

In this work, we present a paradigm bridging electromagnetic (EM) and molecular communication through a stimuli-responsive intra-body model. It has been established that protein molecules, which play a key role in governing cell behavior, can be selectively stimulated using Terahertz (THz) band frequencies. By triggering protein vibrational modes using THz waves, we induce changes in protein conformation, resulting in the activation of a controlled cascade of biochemical and biomechanical events. To analyze such an interaction, we formulate a communication system composed of a nanoantenna transmitter and a protein receiver. We adopt a Markov chain model to account for protein stochasticity with transition rates governed by the nanoantenna force. Both two-state and multi-state protein models are presented to depict different biological configurations. Closed form expressions for the mutual information of each scenario is derived and maximized to find the capacity between the input nanoantenna force and the protein state. The results we obtain indicate that controlled protein signaling provides a communication platform for information transmission between the nanoantenna and the protein with a clear physical significance. The analysis reported in this work should further research into the EM-based control of protein networks.

Collective excitations in α-helical protein structures interacting with the water environment

Low-frequency vibrational excitations of protein macromolecules in the terahertz frequency region are suggested to contribute to many biological processes such as enzymatic catalysis, intra-protein energy/charge transport, recognition, and allostery. To explain high effectiveness of these processes, two possible mechanisms of the long-lived excitation were proposed by H. Fröhlich and A.S. Davydov, which relate to either vibrational modes or solitary waves, respectively. In this paper, we developed a quantum dynamic model of vibrational excitation in α-helical proteins interacting with the aqueous environment. In the model, we distinguished three coupled subsystems, i.e., (i) a chain of hydrogen-bonded peptide groups (PGs), interacting with (ii) the subsystem of the side-chain residuals which in turn interact with (iii) the environment, surrounding water responsible for dissipation and fluctuation in the system. It was shown that the equation of motion for phonon variables of the PG chain can be transformed to nonlinear Schrodinger equation which admits bifurcation into the solution corresponding to the weak-damped vibrational modes (Fröhlich-type regime) and Davydov solitons. A bifurcation parameter is derived through the strength of phonon-phonon interaction between the side-chains and hydration-shell water molecules. As shown, the energy of these excited states is pumped through the interaction of the side-chains with fluctuating water environment of the proteins. The suggested mechanism of the collective vibrational mode excitation is discussed in connection with the recent experiments on the long-lived collective protein excitations in the terahertz frequency region and vibrational energy transport pathways in proteins.

Terahertz Spectroscopy and Imaging: A Cutting-Edge Method for Diagnosing Digestive Cancers

The Terahertz's wavelength is located between the microwave and the infrared region of the electromagnetic spectrum. Because it is non-ionizing and non-invasive, Terahertz (THz)-based detection represents a very attractive tool for repeated assessments, patient monitoring, and follow-up. Cancer acts as the second leading cause of death in many regions, and current predictions estimate a continuous increasing trend. Of all types of tumors, digestive cancers represent an important percentage and their incidence is expected to increase more rapidly than other tumor types due to unhealthy lifestyle habits. Because it can precisely differentiate between different types of molecules, depending on water content, the information obtained through THz-based scanning could have several uses in the management of cancer patients and, more importantly, in the early detection of different solid tumors. The purpose of this manuscript is to offer a comprehensive overview of current data available on THz-based detection for digestive cancers. It summarizes the characteristics of THz waves and their interaction with tissues and subsequently presents available THz-based technologies (THz spectroscopy, THz-tomography, and THZ-endoscope) and their potential for future clinical use. The third part of the review is focused on highlighting current in vitro and in vivo research progress in the field, for identifying specific digestive cancers known as oral, esophageal, gastric, colonic, hepatic, and pancreatic tumors.

Broadband terahertz recognizing conformational characteristics of a significant neurotransmitter γ-aminobutyric acid

γ-Aminobutyric acid (GABA) is the chief inhibitory neurotransmitter in the central nervous system, its conformational behavior is critical for selective biological functions and the process of signal transmission. Although this neuroactive molecule has been extensively studied, its vibrational properties related to the conformation and intermolecular interactions in the terahertz (THz) band have not been identified experimentally yet. In this study, we applied a broadband THz time-domain spectroscopy (THz-TDS) system from 0.5 to 18 THz to characterize a unique THz fingerprint of GABA. The density functional theory calculation results agree well with the THz experimental spectrum. The study shows that the vibrational modes of GABA at 1.15 and 1.39 THz originate from distinct collective vibrations. The absorptions at the higher THz frequencies also carry part of collective vibrations, but more reflect the specific and local vibrational information, including the skeleton deformation and the rocking of the functional groups, which are closely associated with the conformation and flexibility of GABA. This study may help to understand the conformational transitions of neurotransmitter molecules and the resonant response to THz waves.

γ-Aminobutyric acid (GABA) is the chief inhibitory neurotransmitter in the central nervous system, its conformational behavior is critical for selective biological functions and the process of signal transmission.

Terahertz-infrared spectroscopy of Shewanella oneidensis MR-1 extracellular matrix

Employing optical spectroscopy we have performed a comparative study of the dielectric response of extracellular matrix and filaments of electrogenic bacteria Shewanella oneidensis MR-1, cytochrome c, and bovine serum albumin. Combining infrared transmission measurements on thin layers with data of the terahertz spectra, we obtain the dielectric permittivity and AC conductivity spectra of the materials in a broad frequency band from a few cm-1 up to 7000 cm-1 in the temperature range from 5 to 300 K. Strong absorption bands are observed in the three materials that cover the range from 10 to 300 cm-1 and mainly determine the terahertz absorption. When cooled down to liquid helium temperatures, the bands in Shewanella oneidensis MR-1 and cytochrome c reveal a distinct fine structure. In all three materials, we identify the presence of liquid bound water in the form of librational and translational absorption bands at ≈ 200 and ≈ 600 cm-1, respectively. The sharp excitations seen above 1000 cm-1 are assigned to intramolecular vibrations.

Effect of Temperature and Hydration Level on Purple Membrane Dynamics Studied Using Broadband Dielectric Spectroscopy from Sub-GHz to THz Regions

To investigate the effects of temperature and hydration on the dynamics of purple membrane (PM), we measured the broadband complex dielectric spectra from 0.5 GHz to 2.3 THz using a vector network analyzer and terahertz time-domain spectroscopy from 233 to 293 K. In the lower temperature region down to 83 K, the complex dielectric spectra in the THz region were also obtained. The complex dielectric spectra were analyzed through curve fitting using several model functions. We found that the hydrated states of one relaxational mode, which was assigned as the coupled motion of water molecules with the PM surface, began to overlap with the THz region at approximately 230 K. On the other hand, the relaxational mode was not observed for the dehydrated state. On the basis of this result, we conclude that the protein-dynamical-transition-like behavior in the THz region is due to the onset of the overlap of the relaxational mode with the THz region. Temperature hysteresis was observed in the dielectric spectrum at 263 K when the hydration level was high. It is suggested that the hydration water behaves similarly to supercooled liquid at that temperature. The third hydration layer may be partly formed to observe such a phenomenon. We also found that the relaxation time is slower than that of a globular protein, lysozyme, and the microscopic environment in the vicinity of the PM surface is suggested to be more heterogeneous than lysozyme. It is proposed that the spectral overlap of the relaxational mode and the low-frequency vibrational mode is necessary for the large conformational change of protein.

Structural changes of Bombyx mori fibroin from silk gland to fiber as evidenced by Terahertz spectroscopy and other methods

Here we investigated the structural changes of silk fibroin during Bombyx mori silkworm spinning and reconstitution process. X-ray diffraction, Fourier transform infrared spectroscopy, polarized optical microscopy, and terahertz (THz) spectroscopy were applied to monitor the structural features of silk fibroin from posterior, middle silk glands, to cocoons, and then to reconstituted silk. Results show that from silk gland to cocoon, fibroin experiences a significant transformation in crystal structure from a typical silk I, to a silk I-rich mixed structure, and finally to a typical silk II state, accompanied with a change in secondary structure from α-helix and random coil structures to preferential orientation β-sheets. Compared with natural silk fibroins, the reconstituted silk fibroin lacks β-sheet conformation and orientation crystallization. Terahertz spectroscopy readily follows these silk fibroin structural changes. Two characteristic peaks for silk fibroin is observed in 2-10 THz. Their strength ratio is strongly correlated with the β-sheet conformation. The absorbance properties in 0.2-2.0THz also significantly change as a function of changing their crystal structures caused by diverse sources. All of these observations will help in the study of overall structure in silk fibroin to understand more completely the fibroin assembly process in natural spinning and reconstitution process.

Terahertz Frequency Spectroscopy to Determine Cold Shock Protein Stability upon Solvation and Evaporation - A Molecular Dynamics Study

Infrared (IR) and Terahertz (THz) spectroscopy simulations were carried out using CHARMM35b2 to determine protein stability. The stabilities of three bacterial cold shock proteins (Csps) originating from mesophiles, thermophiles and hyper- thermophiles respectively were investigated in this study. The three different Csps were investigated by Normal-Mode analysis and Molecular Dynamics simulation of THz spectra using the Hessian matrix for solvated systems, interpreted in the harmonic approximation at optimum near-melting temperatures of each homologue, by incorporating differences in the hydrous and anhydrous states of the Csps. The results show slight variations in the large scale protein motion. However, the IR spectra of Csps observed at the low frequency saddle surface region, clearly distinguishes the thermophilic and mesophilic proteins based on their stability. Further studies on protein stability employing low-frequency collective modes have the potential to reveal functionally important conformational changes that are biologically significant.

Vibrational resonance, allostery, and activation in rhodopsin-like G protein-coupled receptors

G protein-coupled receptors are a large family of membrane proteins activated by a variety of structurally diverse ligands making them highly adaptable signaling molecules. Despite recent advances in the structural biology of this protein family, the mechanism by which ligands induce allosteric changes in protein structure and dynamics for its signaling function remains a mystery. Here, we propose the use of terahertz spectroscopy combined with molecular dynamics simulation and protein evolutionary network modeling to address the mechanism of activation by directly probing the concerted fluctuations of retinal ligand and transmembrane helices in rhodopsin. This approach allows us to examine the role of conformational heterogeneity in the selection and stabilization of specific signaling pathways in the photo-activation of the receptor. We demonstrate that ligand-induced shifts in the conformational equilibrium prompt vibrational resonances in the protein structure that link the dynamics of conserved interactions with fluctuations of the active-state ligand. The connection of vibrational modes creates an allosteric association of coupled fluctuations that forms a coherent signaling pathway from the receptor ligand-binding pocket to the G-protein activation region. Our evolutionary analysis of rhodopsin-like GPCRs suggest that specific allosteric sites play a pivotal role in activating structural fluctuations that allosterically modulate functional signals.

Surface-enhanced terahertz spectroscopy using gold rod structures resonant with terahertz waves

Terahertz (THz) spectroscopy is a promising method to measure the spectrum of low-frequency modes of molecules or ensembles, such as crystals and polymers, including proteins. However, the main drawback of THz spectroscopy is its extremely low sensitivity. In the present study, we report on signal enhancement in THz spectroscopy achieved by depositing amino acid molecules or their derivatives on a gold rod structured silicon substrate whose localized surface plasmon resonance is exhibited in the THz frequency region. The distinct peaks derived from the enhancement of the inherent spectrum based on a molecular crystal were clearly observed when a longitudinal plasmon resonance mode of the gold rod structure was excited and the plasmon resonance band overlapped the molecular/intermolecular vibrational mode. We discuss the mechanism by which surface-enhanced THz spectroscopy was induced from the viewpoint of the enhancement of light-matter coupling due to plasmon excitation and the modulation of the plasmon band by dipole coupling between the plasmon dipole and molecular/intermolecular vibrational modes.

Terahertz spectroscopy: a powerful new tool for the chemical sciences?

Terahertz spectroscopy is only now beginning to make its transition from initial development by physicists and engineers to broader use by chemists, materials scientists and biologists, thanks to the increasing availability of commercial terahertz spectrometers. With the unique insights that terahertz spectroscopy can provide into intermolecular bonding and crystalline matter, it could prove to be an invaluable addition to the chemist's analytical toolset. This tutorial review aims to give an introduction to terahertz spectroscopy, its techniques, equipment, current applications and potential for the chemical sciences to a broad readership.

Dielectric relaxation dynamics of water in model membranes probed by terahertz spectroscopy

We study hydrated model membranes, consisting of stacked bilayers of 1,2-dioleoyl-sn-glycero-3-phosphocholine lipids, using terahertz time-domain spectroscopy and infrared spectroscopy. Terahertz spectroscopy enables the investigation of water dynamics, owing to its sensitivity to dielectric relaxation processes associated with water reorientation. By controlling the number of water molecules per lipid molecule in the system, we elucidate how the interplay between the model membrane and water molecules results in different water dynamics. For decreasing hydration levels, we observe the appearance of new types of water dynamics: the collective bulklike dynamics become less pronounced, whereas an increased amount of both very slowly reorienting (i.e., irrotational) and very rapidly reorienting (i.e., fast) water molecules appear. Temperature-dependent measurements reveal the interconversion between the three distinct types of water present in the system.

Rhodopsin: the functional significance of asn-linked glycosylation and other post-translational modifications

Rhodopsin, the G-protein coupled receptor in retinal rod photoreceptors, is a highly conserved protein that undergoes several types of post-translational modifications. These modifications are essential to maintain the protein's structure as well as its proper function in the visual transduction cycle. Rhodopsin is N-glycosylated at Asn-2 and Asn-15 in its extracellular N-terminal domain. Mutations within the glycosylation consensus sequences of rhodopsin cause autosomal dominant retinitis pigmentosa, a disease that leads to blindness. Several groups have studied the role of rhodopsin's N-linked glycan chains in protein structure and function using a variety of approaches. These include the generation of a transgenic mouse model, study of a naturally occurring mutant animal model, in vivo pharmacological inhibition of glycosylation, and in vitro analyses using transfected COS-1 cells. These studies have provided insights into the possible role of rhodopsin glycosylation, but have yielded conflicting results.

A dynamics study of the A-chain of ricin by terahertz vibrational calculation and normal modes analysis

We studied the terahertz (THz) spectroscopy and low frequency normal modes of both apo- and holo- (adenosine monophosphate (AMP)-bound) ricin-A-chain (RTA) as a means to understand the dynamical changes that RTA undergoes upon substrate binding. The calculated THz spectra of apo- and holo-RTAs demonstrated a general intensity suppression upon substrate binding, which is attributed to the reduced number of collective motion in THz region. In normal mode analysis of RTA we find a shearing motion that is shared by both the apo- and holo-RTAs, whereas a breathing motion, and an upward hinge rising and an alpha-G bending characteristic motion are dampened significantly upon AMP binding, suggesting these motions are involved in the necessary flexibility of the active site. In contrast, we find a normal mode motion that separates domains I and II of RTA at the interface that is more common in the holo-protein. We hypothesized that the flexibility of the entrance of RTA can facilitate the entry of rRNA and allow the substrate to adjust its conformation and orientation prior to depurination. This process suggests an rRNA binding pathway which is supplemental the current RTA depurination mechanism.