Broadband Dielectric Spectroscopy on Lysozyme in the Sub-Gigahertz to Terahertz Frequency Regions: Effects of Hydration and Thermal Excitation

Nonthermal acceleration of protein hydration by sub-terahertz irradiation

The collective intermolecular dynamics of protein and water molecules, which overlap in the sub-terahertz (THz) frequency region, are relevant for expressing protein functions but remain largely unknown. This study used dielectric relaxation (DR) measurements to investigate how externally applied sub-THz electromagnetic fields perturb the rapid collective dynamics and influence the considerably slower chemical processes in protein-water systems. We analyzed an aqueous lysozyme solution, whose hydration is not thermally equilibrated. By detecting time-lapse differences in microwave DR, we demonstrated that sub-THz irradiation gradually decreases the dielectric permittivity of the lysozyme solution by reducing the orientational polarization of water molecules. Comprehensive analysis combining THz and nuclear magnetic resonance spectroscopies suggested that the gradual decrease in the dielectric permittivity is not induced by heating but is due to a slow shift toward the hydrophobic hydration structure in lysozyme. Our findings can be used to investigate hydration-mediated protein functions based on sub-THz irradiation.

Incoherent Neutron Scattering and Terahertz Time-Domain Spectroscopy on Protein and Hydration Water

Incoherent inelastic and quasi-elastic neutron scattering (INS) and terahertz time-domain spectroscopy (THz-TDS) are spectroscopy methods that directly detect molecular dynamics, with an overlap in the measured energy regions of each method. Due to the different characteristics of their probes (i.e., neutron and light), the information obtained and the sample conditions suitable for each method differ. In this review, we introduce the differences in the quantum beam properties of the two methods and their associated advantages and disadvantages in molecular spectroscopy. Neutrons are scattered via interaction with nuclei; one characteristic of neutron scattering is a large incoherent scattering cross-section of a hydrogen atom. INS records the auto-correlation functions of atomic positions. By using the difference in neutron scattering cross-sections of isotopes in multi-component systems, some molecules can be selectively observed. In contrast, THz-TDS observes the cross-correlation function of dipole moments. In water-containing biomolecular samples, the absorption of water molecules is particularly large. While INS requires large-scale experimental facilities, such as accelerators and nuclear reactors, THz-TDS can be performed at the laboratory level. In the analysis of water molecule dynamics, INS is primarily sensitive to translational diffusion motion, while THz-TDS observes rotational motion in the spectrum. The two techniques are complementary in many respects, and a combination of the two is very useful in analyzing the dynamics of biomolecules and hydration water.

Water dynamics in the hydration shell of hyper-branched poly-ethylenimine

We performed THz and GHz dielectric relaxation spectroscopy to investigate the reorientational dynamics of water molecules in the hydration shell of amphiphilic hyper-branched poly-ethylenimine (HPEI). Four Debye equations were employed to describe four types of water in the hydration shell, including bulk-like water, under-coordinated water, slow water (water molecules hydrating the hydrophobic groups and water molecules accepting hydrogen bonds from the NH₂ groups) and super slow water (water molecules donating hydrogen bonds to and accepting hydrogen bonds from NH groups). The time scales of undercoordinated and bulk-like water show a slight decline from 0.4 to 0.1 ps and from 8 to 2 ps, respectively. Because of hydrophilic amino groups, HPEI molecules exhibit a strong retardation effect, where the time scales of slow and super slow water increase with concentration from 17 to 39.9 ps and from 88 to 225 ps, respectively.

Stochastic Analysis of Molecular Dynamics Reveals the Rotation Dynamics Distribution of Water around Lysozyme

Water dynamics is essential to biochemical processes by mediating all such reactions, including biomolecular degeneration in solutions. To disentangle the molecular-scale distribution of water dynamics around a solute biomolecule, we investigated here the rotational dynamics of water around lysozyme by combining molecular dynamics (MD) simulations and broadband dielectric spectroscopy (BDS). A statistical analysis using the relaxation times and trajectories of every single water molecule was proposed, and the two-dimensional probability distribution of water at a distance from the lysozyme surface with a rotational relaxation time was given. For the observed lysozyme solutions of 34-284 mg/mL, we discovered that the dielectric relaxation time obtained from this distribution agrees well with the measured γ relaxation time, which suggests that rotational self-correlation of water molecules underlies the gigahertz domain of the dielectric spectra. Regardless of protein concentration, water rotational relaxation time versus the distance from the lysozyme surface revealed that the water rotation is severely retarded within 3 Å from the lysozyme surface and is nearly comparable to pure water when farther than 10 Å. The dimension of the first hydration layer was subsequently identified in terms of the relationship between the acceleration of water rotation and the distance from the protein surface.

Functional-State Dependence of Picosecond Protein Dynamics

We examine temperature-dependent picosecond dynamics of two benchmarking proteins lysozyme and cytochrome c using temperature-dependent terahertz permittivity measurements. We find that a double Arrhenius temperature dependence with activation energies E₁ ∼ 0.1 kJ/mol and E₂ ∼ 10 kJ/mol fits the folded and ligand-free state response. The higher activation energy is consistent with the so-called protein dynamical transition associated with beta relaxations at the solvent-protein interface. The lower activation energy is consistent with correlated structural motions. When the structure is removed by denaturing, the lower-activation-energy process is no longer present. Additionally, the lower-activation-energy process is diminished with ligand binding but not for changes in the internal oxidation state. We suggest that the lower-energy activation process is associated with collective structural motions that are no longer accessible with denaturing or binding.

Correlating solvation with conformational pathways of proteins in alcohol-water mixtures: a THz spectroscopic insight

Water, being an active participant in most of the biophysical processes, is important to trace how protein solvation changes as its conformation evolves in the presence of solutes or co-solvents. In this study, we investigate how the secondary structures of two diverse proteins - lysozyme and β-lactoglobulin - change in the aqueous mixtures of two alcohols - ethanol and 2,2,2-trifluoroethanol (TFE) using circular dichroism measurements. We observe that these alcohols change the secondary structures of these proteins and the changes are protein-specific. Subsequently, we measure the collective solvation dynamics of these two proteins both in the absence and in the presence of alcohols by measuring the frequency-dependent absorption coefficient (α(ν)) in the THz (0.1-1.2 THz) frequency domain. The alcohol-water mixtures exhibit a non-ideal behaviour with the highest absorption difference (Δα) obtained at X_alcohol = 0.2. The protein solvation in the presence of the alcohols shows an oscillating behaviour in which Δα_protein changes with X_alcohol. Such an oscillatory behaviour of protein solvation results from a delicate interplay between the protein-water, protein-alcohol and water-alcohol associations. We attempt to correlate the various structural conformations of the proteins with the associated solvation.

THz irradiation inhibits cell division by affecting actin dynamics

Biological phenomena induced by terahertz (THz) irradiation are described in recent reports, but underlying mechanisms, structural and dynamical change of specific molecules are still unclear. In this paper, we performed time-lapse morphological analysis of human cells and found that THz irradiation halts cell division at cytokinesis. At the end of cytokinesis, the contractile ring, which consists of filamentous actin (F-actin), needs to disappear; however, it remained for 1 hour under THz irradiation. Induction of the functional structures of F-actin was also observed in interphase cells. Similar phenomena were also observed under chemical treatment (jasplakinolide), indicating that THz irradiation assists actin polymerization. We previously reported that THz irradiation enhances the polymerization of purified actin in vitro; our current work shows that it increases cytoplasmic F-actin in vivo. Thus, we identified one of the key biomechanisms affected by THz waves.

Photo-Switching of Protein Dynamical Collectivity

We examine changes in the picosecond structural dynamics with irreversible photobleaching of red fluorescent proteins (RFP) mCherry, mOrange2 and TagRFP-T. Measurements of the protein dynamical transition using terahertz time-domain spectroscopy show in all cases an increase in the turn-on temperature in the bleached state. The result is surprising given that there is little change in the protein surface, and thus, the solvent dynamics held responsible for the transition should not change. A spectral analysis of the measurements guided by quasiharmonic calculations of the protein absorbance reveals that indeed the solvent dynamical turn-on temperature is independent of the thermal stability/photostate however the protein dynamical turn-on temperature shifts to higher temperatures. This is the first demonstration of switching the protein dynamical turn-on temperature with protein functional state. The observed shift in protein dynamical turn-on temperature relative to the solvent indicates an increase in the required mobile waters necessary for the protein picosecond motions, that is, these motions are more collective. Melting-point measurements reveal that the photobleached state is more thermally stable, and structural analysis of related RFP’s shows that there is an increase in internal water channels as well as a more uniform atomic root mean squared displacement. These observations are consistent with previous suggestions that water channels form with extended light excitation providing O2 access to the chromophore and subsequent fluorescence loss. We report that these same channels increase internal coupling enhancing thermal stability and collectivity of the picosecond protein motions. The terahertz spectroscopic characterization of the protein and solvent dynamical onsets can be applied generally to measure changes in collectivity of protein motions.

Nonthermal excitation effects mediated by sub-terahertz radiation on hydrogen exchange in ubiquitin

Water dynamics in the hydration layers of biomolecules play crucial roles in a wide range of biological functions. A hydrated protein contains multiple components of diffusional and vibrational dynamics of water and protein, which may be coupled at ∼0.1-THz frequency (10-ps timescale) at room temperature. However, the microscopic description of biomolecular functions based on various modes of protein-water-coupled motions remains elusive. A novel approach for perturbing the hydration dynamics in the subterahertz frequency range and probing them at the atomic level is therefore warranted. In this study, we investigated the effect of klystron-based, intense 0.1-THz excitation on the slow dynamics of ubiquitin using NMR-based measurements of hydrogen-deuterium exchange. We demonstrated that the subterahertz irradiation accelerated the hydrogen-deuterium exchange of the amides located in the interior of the protein and hydrophobic surfaces while decelerating this exchange in the amides located in the surface loop and short 3₁₀ helix regions. This subterahertz-radiation-induced effect was qualitatively contradictory to the increased-temperature-induced effect. Our results suggest that the heterogeneous water dynamics occurring at the protein-water interface include components that are nonthermally excited by the subterahertz radiation. Such subterahertz-excited components may be linked to the slow function-related dynamics of the protein.

Freezable and Unfreezable Hydration Water: Distinct Contributions to Protein Dynamics Revealed by Neutron Scattering

Hydration water plays a crucial role for activating the protein dynamics required for functional expression. Yet, the details are not understood about how hydration water couples with protein dynamics. A temperature hysteresis of the ice formation of hydration water is a key phenomenon to understand which type of hydration water, unfreezable or freezable hydration water, is crucial for the activation of protein dynamics. Using neutron scattering, we observed a temperature-hysteresis phenomenon in the diffraction peaks of the ice of freezable hydration water, whereas protein dynamics did not show any temperature hysteresis. These results show that the protein dynamics is not coupled with freezable hydration water dynamics, and unfreezable hydration water is essential for the activation of protein dynamics. Decoupling of the dynamics between unfreezable and freezable hydration water could be the cause of the distinct contributions of hydration water to protein dynamics.

Hemoglobin Dynamics in Solution vis-à-vis Under Confinement: An Electrochemical Perspective

Confining heme protein in silico often leads to beneficial functionalities such as an enhanced electrochemical response from the heme center. This can be harnessed to design effective biosensors for medical diagnostics. Proteins under confinement, surface confinement on the electrode to be precise, have more ordered and monodisperse structure compared to the protein in bulk solution. As the electrochemical response of a protein comes from those protein molecules that are confined within the electrical double layer across the electrode-electrolyte interface, it is expected that restriction of conformational fluctuations of the polymeric protein will help in enhancement of the electrochemical response. This is probably the prima facie reason for electrochemical response enhancement under confinement. We examine the dynamic features of hemoglobin under confinement vis-à-vis that in bulk solution. We use a variety of spectroscopic techniques across a wide time-space window to establish the following facts: (a) hardening of the protein polypeptide backbone, (b) slowing down of protein diffusion, (c) increase in relaxation times in NMR, and (d) slowing down of dielectric relaxation times under confinement. This indicates an overall quenching of protein dynamics when the protein is confined inside silica matrix. Thus, we hypothesize that along with retention of secondary structure, this quenching of dynamics contributes to the enhancement of electrochemical response observed.

Length-scale dependence of protein hydration-shell density

Here we present a computational approach based on molecular dynamics (MD) simulation to study the dependence of the protein hydration-shell density on the size of the protein molecule. The hydration-shell density of eighteen different proteins, differing in size, shape and function (eight of them are antifreeze proteins), is calculated. The results obtained show that an increase in the hydration-shell density, relative to that of the bulk, is observed (in the range of 4-14%) for all studied proteins and that this increment strongly correlates with the protein size. In particular, a decrease in the density increment is observed for decreasing protein size. A simple model is proposed in which the basic idea is to approximate the protein molecule as an effective ellipsoid and to partition the relevant parameters, i.e. the solvent-accessible volume and the corresponding solvent density, into two regions: inside and outside the effective protein ellipsoid. It is found that, within the model developed here, almost all of the hydration-density increase is located inside the protein ellipsoid, basically corresponding to pockets within, or at the surface of the protein molecule. The observed decrease in the density increment is caused by the protein size only and no difference is found between antifreeze and non-antifreeze proteins.

Broadband dielectric spectroscopy from sub GHz to THz of hydrated lipid bilayer of DMPC

In order to study the dynamics of a phospholipid and its hydration water, we measured complex dielectric spectra of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) from the sub-GHz to the THz frequency region with varying the temperature and hydration level of the sample. Spectra obtained from a vector network analyzer and two terahertz time-domain spectrometers are adjusted, which enables us to analyze the dielectric spectra from the sub-GHz region to the THz region by a model function. We confirmed a fast relaxational mode in the sub-THz region, which was suggested by the previous work which only used a THz spectrometer.

Surface Topography Effects of Globular Biomolecules on Hydration Water

Hydration water serves as a microscopic manifestation of structural stability and functions of biomolecules. To develop bio-nanomaterials in applications, it is important to study how the surface topography and heterogeneity of biomolecules result in their diversity of the hydration dynamics and energetics. We here performed molecular dynamics simulations combined with the steered molecular dynamics and umbrella sampling to investigate the dynamics and escape process associated with the free energy change of water molecules close to a globular biomolecule, i.e., hemoglobin (Hb) and G-quadruplex DNA (GDNA). The residence time, power of long-time tail, and dipole relaxation time were found to display drastic changes within the averaged hydration shell of 3.0-5.0 Å. Compared with bulk water, in the inner hydration shell, the water dipole moment displays a slower relaxation process and is more oriented toward GDNA than toward Hb, forming a hedgehog-like structure when it surrounds GDNA. In particular, a spine water structure is observed in the GDNA narrow groove. The water isotope effect not only prolongs the dynamic time scales of libration motion in the inner hydration shell and the dipole relaxation processes in the bulk but also strengthens the DNA spine water structure. The potential of the mean force profile reflects the integrity of the hydration shell structure and enables us to obtain detailed insights into the structures formed by water, such as the caged H-bond network and the edge bridge structures; it also reveals that local hydration shell free energy (LHSFE) depends on H-bond rupture processes and ranges from 0.2 to 4.2 kcal/mol. Our results demonstrate that the surface topography of a biomolecule influences the integrity of the hydration shell structure and LHSFE. Our studies are able to identify various further applications in the areas of microfluid devices and nano-dewetting on bioinspired surfaces.

Role of Disulfide Bonds and Topological Frustration in the Kinetic Partitioning of Lysozyme Folding Pathways

Disulfide bonds in proteins can strongly influence the folding pathways by constraining the conformational space. Lysozyme has four disulfide bonds and is widely studied for its antibacterial properties. Experiments on lysozyme infer that the protein folds through a fast and a slow pathway. However, the reasons for the kinetic partitioning in the folding pathways are not completely clear. Using a coarse-grained protein model and simulations, we show that two out of the four disulfide bonds, which are present in the α-domain of lysozyme, are responsible for the slow folding pathway. In this pathway, a kinetically trapped intermediate state, which is close to the native state, is populated. In this state, the orientations of α-helices present in the α-domain are misaligned relative to each other. The protein in this state has to partially unfold by breaking down the interhelical contacts between the misaligned helices to fold to the native state. However, the topological constraints due to the two disulfide bonds present in the α-domain make the protein less flexible, and it is trapped in this conformation for hundreds of milliseconds. On disabling these disulfide bonds, we find that the kinetically trapped intermediate state and the slow folding pathway disappear. Simulations mimicking the folding of protein without disulfide bonds under oxidative conditions show that the native disulfide bonds are formed as the protein folds, indicating that folding guides the formation of disulfide bonds. The sequence of formation of the disulfide bonds is Cys64-Cys80 → Cys76-Cys94 → Cys30-Cys115 → Cys6-Cys127. Any disulfide bond that forms before its precursor in the sequence has to break and follow the sequence for the protein to fold. These results show that lysozyme also serves as a very good model system to probe the role of disulfide bonds and topological frustration in protein folding. The predictions from the simulations can be verified by single-molecule fluorescence resonance energy transfer or single-molecule pulling experiments, which can probe heterogeneity in the folding pathways.

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.

High-Precision Megahertz-to-Terahertz Dielectric Spectroscopy of Protein Collective Motions and Hydration Dynamics

The low-frequency collective vibrational modes in proteins as well as the protein-water interface have been suggested as dominant factors controlling the efficiency of biochemical reactions and biological energy transport. It is thus crucial to uncover the mystery of the hydration structure and dynamics as well as their coupling to collective motions of proteins in aqueous solutions. Here, we report dielectric properties of aqueous bovine serum albumin protein solutions as a model system using an extremely sensitive dielectric spectrometer with frequencies spanning from megahertz to terahertz. The dielectric relaxation spectra reveal several polarization mechanisms at the molecular level with different time constants and dielectric strengths, reflecting the complexity of protein-water interactions. Combining the effective-medium approximation and molecular dynamics simulations, we have determined collective vibrational modes at terahertz frequencies and the number of water molecules in the tightly bound and loosely bound hydration layers. High-precision measurements of the number of hydration water molecules indicate that the dynamical influence of proteins extends beyond the first solvation layer, to around 7 Å distance from the protein surface, with the largest slowdown arising from water molecules directly hydrogen-bonded to the protein. Our results reveal critical information of protein dynamics and protein-water interfaces, which determine biochemical functions and reactivity of proteins.

Boson peak dynamics of natural polymer starch investigated by terahertz time-domain spectroscopy and low-frequency Raman scattering

Terahertz time-domain spectroscopy and low-frequency Raman scattering were performed on the natural polymer starch to investigate the boson peak (BP) dynamics. In the infrared spectrum, the BP was observed at 0.99THz at the lowest temperature. Compared to the result from a previous study for vitreous glucose, both the frequency of the BP and absorption coefficient show lower values than those of the vitreous glucose. These behaviors originate from the longer correlation length of the medium-range order and lower concentration of hydroxyl groups in the starch. In the Raman spectrum, the BP was observed at 1.1THz at room temperature, although the BP was not observed around room temperature due to the excess wing of the fast relaxation modes in the infrared spectrum. The temperature dependence of ε″(ν) during the heating process and cooling process shows a hysteresis below 230K. During the heating process, kinks were observed at 140K and 230K. These kinks are attributed to the β-relaxation and the β_wet-relaxation, respectively.

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.

Time-Domain THz Spectroscopy Reveals Coupled Protein-Hydration Dielectric Response in Solutions of Native and Fibrils of Human Lysozyme

Here we reveal details of the interaction between human lysozyme proteins, both native and fibrils, and their water environment by intense terahertz time domain spectroscopy. With the aid of a rigorous dielectric model, we determine the amplitude and phase of the oscillating dipole induced by the THz field in the volume containing the protein and its hydration water. At low concentrations, the amplitude of this induced dipolar response decreases with increasing concentration. Beyond a certain threshold, marking the onset of the interactions between the extended hydration shells, the amplitude remains fixed but the phase of the induced dipolar response, which is initially in phase with the applied THz field, begins to change. The changes observed in the THz response reveal protein-protein interactions mediated by extended hydration layers, which may control fibril formation and may have an important role in chemical recognition phenomena.