Collective Dynamics of Lysozyme in Water: Terahertz Absorption Spectroscopy and Comparison with Theory

Quantitative Measurements of DP in Cellulose Paper Based on Terahertz Spectroscopy

The power transformer is vital to the reliability of the power grid which is most commonly insulated with Kraft paper and immersed in mineral oil, among which the aged state of the paper is mainly correlated to the operating life of the transformer. Degree of polymerization (DP) is a direct parameter to assess the aged condition of insulating paper, but existing DP measurement by viscosity methods are destructive and complicated. In this paper, terahertz time-domain spectroscopy (THz-TDS) was introduced to reach rapid, non-destructive detection of the DP of insulating paper. The absorption spectra of insulating paper show that characteristic peak regions at 1.8 and 2.23 THz both exhibit a log-linear quantitative relationship with DP, and their universalities are confirmed by conducting the above relationship on different types of insulating paper. Fourier transform infrared spectroscopy (FTIR) analysis and molecular dynamics modeling further revealed that 1.8 and 2.23 THz were favorably associated with the growth of water-cellulose hydrogen bond strength and amorphous cellulose, respectively. This paper demonstrates the viability of applying THz-TDS to the non-destructive detection of DP in insulating paper and assigned the vibration modes of the characteristic absorption peaks.

Toward in-process technology-aided automation for enhanced microbial food safety and quality assurance in milk and beverages processing

Ensuring the safety of food products is critical to food production and processing. In food processing and production, several standard guidelines are implemented to achieve acceptable food quality and safety. This notwithstanding, due to human limitations, processed foods are often contaminated either with microorganisms, microbial byproducts, or chemical agents, resulting in the compromise of product quality with far-reaching consequences including foodborne diseases, food intoxication, and food recall. Transitioning from manual food processing to automation-aided food processing (smart food processing) which is guided by artificial intelligence will guarantee the safety and quality of food. However, this will require huge investments in terms of resources, technologies, and expertise. This study reviews the potential of artificial intelligence in food processing. In addition, it presents the technologies and methods with potential applications in implementing automated technology-aided processing. A conceptual design for an automated food processing line comprised of various operational layers and processes targeted at enhancing the microbial safety and quality assurance of liquid foods such as milk and beverages is elaborated.

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.

Kinetic coherence underlies the dynamics of disordered proteins

The dynamics of two proteins of similar size, the globular lysozyme and the intrinsically disordered Huntingtin interacting protein, has been simulated in three states resembling a globule, a pre-molten globule, and a molten globule. A coherence time τ has been defined, measuring the delay in the display of a stochastic behaviour after a perturbation of the system. This time has been computed for two sets of collective variables: the projection of the phase point onto the positions and momenta subspaces (τ_r and τ_p), and the principal components (PCs) of positions q and momenta π produced by a covariance analysis in these subspaces (τ_q and τ_π). In all states τ_p ≈ 3.5τ_r, and τ_π ≈ 3.5τ_q. The coherence times of individual PCs, τ^(l)_q and τ^(l)_π, have also been computed, and τ^(l)_π > τ^(l)_q in all states. The prevalence of τ_p over τ_r, or of τ_π over τ_q, drives the dynamics of the protein over a time range of ≈1–2 ps; moreover, a hidden synchronism appears to raise the momenta subspace's coherence above that of its individual PCs. In the transition of lysozyme to the molten globule the τ^(l)_q decrease but, unexpectedly, the τ^(l)_π increase; after this transition τ_p ≈ 5τ_r and τ_π ≈ 5τ_q. A gain of kinetic coherence accompanies thus the loss of structural coherence caused by the denaturation of the protein in the transition from globule to molten globule. The increase of the τ^(l)_π does not take place in the analogous transition of the Huntingtin protein. These results are compared with those of a similar analysis performed on three pseudo-proteins designed by scrambling the primary sequence of the Huntingtin interacting protein, and on two oligopeptides. The hidden synchronism appears to be a generic property of these polypeptides. The τ^(l)_π spectrum is similar in denaturated and in intrinsically disordered biomolecules; but the gain of kinetic coherence as a result of denaturation seems to be a specific property of the biologically functional lysozyme.

In the phase space of a globular or intrinsically disordered protein, the momenta's dynamics is less chaotic than the coordinates' dynamics. When a protein is denaturated, a gain in kinetic coherence accompanies the loss of structural coherence.

Solvent dynamics play a decisive role in the complex formation of biologically relevant redox proteins

Electron transfer processes between proteins are vital in many biological systems. Yet, the role of the solvent in influencing these redox reactions remains largely unknown. In this study, terahertz-time domain spectroscopy (THz-TDS) is used to probe the collective hydration dynamics of flavoenzyme ferredoxin-NADP+-reductase (FNR), electron transfer protein ferredoxin-1 (PetF), and the transient complex that results from their interaction. Results reveal changes in the sub-picosecond hydration dynamics that are dependent upon the surface electrostatic properties of the individual proteins and the transient complex. Retarded solvent dynamics of 8-9 ps are observed for FNR, PetF, and the FNR:PetF transient complex. Binding of the FNR:PetF complex to the substrate NADP+ results in bulk-like solvent dynamics of 7 ps, showing that formation of the ternary complex is entropically favored. Our THz measurements reveal that the electrostatic interaction of the protein surface with water results in charge sensitive changes in the solvent dynamics. Complex formation between the positively charged FNR:NADP+ pre-complex and the negatively charged PetF is not only entropically favored, but in addition the solvent reorganization into more bulk-like water assists the molecular recognition process. The change in hydration dynamics observed here suggests that the interaction with the solvent plays a significant role in mediating electron transfer processes between proteins.

A Nondestructive Eggshell Thickness Measurement Technique Using Terahertz Waves

Eggshells play a number of important roles in the avian and reptile kingdom: protection of internal contents and as a major source of minerals for developing embryos. However, when researching these respective roles, eggshell thickness measurement remains a bottleneck due to the lack of a non-destructive measurement techniques. As a result, many avian and reptile research protocols omit consideration of eggshell thickness bias on egg or embryo growth and development. Here, we validate a non-destructive method to estimate eggshell thickness based on terahertz (THz) reflectance spectroscopy using chicken white coloured eggs. Since terahertz waves are reflected from outer air-eggshell interface, as well as the inner eggshell-membrane boundary, the resulting interference signals depend on eggshell thickness. Thus, it is possible to estimate shell thickness from the oscillation distance in frequency-domain. A linear regression-based prediction model for non-destructive eggshell thickness measurement was developed, which had a coefficient of determination (R²) of 0.93, RMSEP of 0.009, RPD of 3.45 and RER 13.67. This model can estimate eggshell thickness to a resolution of less than 10 μm. This method has the potential to expand the protocols in the field of avian and reptile research, as well as be applied to industrial grading of eggs.

A comparative study of the low-frequency vibrations of l-histidine molecule in different solid states

Low-frequency vibrations of l-histidine in a neat solid state and its monohydrochloride monohydrate are both investigated using THz spectroscopy and DFT calculations. The molecular motions in those modes are further quantitatively decomposed into a number of submotions and discussed in terms of their contributions to a mode. The results show significant differences in the averaged contribution percentage of intermolecular motions and the dihedral angle distortions of the imidazole ring between these two crystals. Those phenomena are interpreted from the viewpoint of their hydrogen-bond configurations.

A Comprehensive Review on Food Applications of Terahertz Spectroscopy and Imaging

Food product safety is a public health concern. Most of the food safety analytical and detection methods are expensive, labor intensive, and time consuming. A safe, rapid, reliable, and nondestructive detection method is needed to assure consumers that food products are safe to consume. Terahertz (THz) radiation, which has properties of both microwave and infrared, can penetrate and interact with many commonly used materials. Owing to the technological developments in sources and detectors, THz spectroscopic imaging has transitioned from a laboratory-scale technique into a versatile imaging tool with many practical applications. In recent years, THz imaging has been shown to have great potential as an emerging nondestructive tool for food inspection. THz spectroscopy provides qualitative and quantitative information about food samples. The main applications of THz in food industries include detection of moisture, foreign bodies, inspection, and quality control. Other applications of THz technology in the food industry include detection of harmful compounds, antibiotics, and microorganisms. THz spectroscopy is a great tool for characterization of carbohydrates, amino acids, fatty acids, and vitamins. Despite its potential applications, THz technology has some limitations, such as limited penetration, scattering effect, limited sensitivity, and low limit of detection. THz technology is still expensive, and there is no available THz database library for food compounds. The scanning speed needs to be improved in the future generations of THz systems. Although many technological aspects need to be improved, THz technology has already been established in the food industry as a powerful tool with great detection and quantification ability. This paper reviews various applications of THz spectroscopy and imaging in the food industry.

Spectroscopic methods for assessing the molecular origins of macroscopic solution properties of highly concentrated liquid protein solutions

In cases of subcutaneous injection of therapeutic monoclonal antibodies, high protein concentrations (>50 mg/ml) are often required. During the development of these high concentration liquid formulations (HCLF), challenges such as aggregation, gelation, opalescence, phase separation, and high solution viscosities are more prone compared to low concentrated protein formulations. These properties can impair manufacturing processes, as well as protein stability and shelf life. To avoid such unfavourable solution properties, a detailed understanding about the nature of these properties and their driving forces are required. However, the fundamental mechanisms that lead to macroscopic solution properties, as above mentioned, are complex and not fully understood, yet. Established analytical methods for assessing the colloidal stability, i.e. the ability of a native protein to remain dispersed in solution, are restricted to dilute conditions and provide parameters such as the second osmotic virial coefficient, B₂₂, and the diffusion interaction coefficient, k_D. These parameters are routinely applied for qualitative estimations and identifications of proteins with challenging solution behaviours, such as high viscosities and aggregation, although the assays are prepared for low protein concentration conditions, typically between 0.1 and 20 mg/ml ("ideal" solution conditions). Quantitative analysis of samples of high protein concentration is difficult and it is hard to obtain information about the driving forces of such solution properties and corresponding protein-protein self-interactions. An advantage of using specific spectroscopic methods is the potential of directly analysing highly concentrated protein solutions at different solution conditions. This allows for collecting/gaining valuable information about the fundamental mechanisms of solution properties of the high protein concentration regime. In addition, the derived parameters might be more predictive as compared to the parameters originating from assays which are optimized for the low protein concentration range. The provided information includes structural data, molecular dynamics at various timescales and protein-solvent interactions, which can be obtained at molecular resolution. Herein, we provide an overview about spectroscopic techniques for analysing the origins of macroscopic solution behaviours in general, with a specific focus on pharmaceutically relevant high protein concentration and formulation conditions.

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.

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.

Using THz Spectroscopy, Evolutionary Network Analysis Methods, and MD Simulation to Map the Evolution of Allosteric Communication Pathways in c-Type Lysozymes

It is now widely accepted that protein function is intimately tied with the navigation of energy landscapes. In this framework, a protein sequence is not described by a distinct structure but rather by an ensemble of conformations. And it is through this ensemble that evolution is able to modify a protein's function by altering its landscape. Hence, the evolution of protein functions involves selective pressures that adjust the sampling of the conformational states. In this work, we focus on elucidating the evolutionary pathway that shaped the function of individual proteins that make-up the mammalian c-type lysozyme subfamily. Using both experimental and computational methods, we map out specific intermolecular interactions that direct the sampling of conformational states and accordingly, also underlie shifts in the landscape that are directly connected with the formation of novel protein functions. By contrasting three representative proteins in the family we identify molecular mechanisms that are associated with the selectivity of enhanced antimicrobial properties and consequently, divergent protein function. Namely, we link the extent of localized fluctuations involving the loop separating helices A and B with shifts in the equilibrium of the ensemble of conformational states that mediate interdomain coupling and concurrently moderate substrate binding affinity. This work reveals unique insights into the molecular level mechanisms that promote the progression of interactions that connect the immune response to infection with the nutritional properties of lactation, while also providing a deeper understanding about how evolving energy landscapes may define present-day protein function.

Salt effects on the picosecond dynamics of lysozyme hydration water investigated by terahertz time-domain spectroscopy and an insight into the Hofmeister series for protein stability and solubility

The addition of salts into protein aqueous solutions causes changes in protein solubility and stability, the ability of which is known to be ordered in the Hofmeister series.

Optical measurements of long-range protein vibrations

Protein biological function depends on structural flexibility and change. From cellular communication through membrane ion channels to oxygen uptake and delivery by haemoglobin, structural changes are critical. It has been suggested that vibrations that extend through the protein play a crucial role in controlling these structural changes. While nature may utilize such long-range vibrations for optimization of biological processes, bench-top characterization of these extended structural motions for engineered biochemistry has been elusive. Here we show the first optical observation of long-range protein vibrational modes. This is achieved by orientation-sensitive terahertz near-field microscopy measurements of chicken egg white lysozyme single crystals. Underdamped modes are found to exist for frequencies >10 cm(-1). The existence of these persisting motions indicates that damping and intermode coupling are weaker than previously assumed. The methodology developed permits protein engineering based on dynamical network optimization.

Using THz time-scale infrared spectroscopy to examine the role of collective, thermal fluctuations in the formation of myoglobin allosteric communication pathways and ligand specificity

In this investigation we use THz time-scale spectroscopy to conduct an initial set of studies on myoglobin with the aim of providing further insight into the global, collective thermal fluctuations in the protein that have been hypothesized to play a prominent role in the dynamic formation of transient ligand channels as well as in shaping the molecular level basis for ligand discrimination. Using the two ligands O2 and CO, we have determined that the perturbation from the heme-ligand complex has a strong influence on the characteristics of the myoglobin collective dynamics that are excited upon binding. Further, the differences detected in the collective protein motions in Mb-O2 compared with those in Mb-CO appear to be intimately tied with the pathways of long-range allosteric communication in the protein, which ultimately determine the trajectories selected by the respective ligands on the path to and from the heme-binding cavity.

Aligning experimental and theoretical anisotropic B-factors: water models, normal-mode analysis methods, and metrics

The strength of X-ray crystallography in providing the information for protein dynamics has been under appreciated. The anisotropic B-factors (ADPs) from high-resolution structures are invaluable in studying the relationship among structure, dynamics, and function. Here, starting from an in-depth evaluation of the metrics used for comparing the overlap between two ellipsoids, we applied normal-mode analysis (NMA) to predict the theoretical ADPs and then align them with experimental results. Adding an extra layer of explicitly treated water on protein surface significantly improved the energy minimization results and better reproduced the anisotropy of experimental ADPs. In comparing experimental and theoretical ADPs, we focused on the overlap in shape, the alignment of dominant directions, and the similarity in magnitude. The choices of water molecules, NMA methods, and the metrics for evaluating the overlap of ADPs determined final results. This study provides useful information for exploring the physical basis and the application potential of experimental ADPs.

THz time scale structural rearrangements and binding modes in lysozyme-ligand interactions

Predicting the conformational changes in proteins that are relevant for substrate binding is an ongoing challenge in the aim of elucidating the functional states of proteins. The motions that are induced by protein-ligand interactions are governed by the protein global modes. Our measurements indicate that the detected changes in the global backbone motion of the enzyme upon binding reflect a shift from the large-scale collective dominant mode in the unbound state towards a functional twisting deformation that assists in closing the binding cleft. Correlated motion in lysozyme has been implicated in enzyme function in previous studies, but detailed characterization of the internal fluctuations that enable the protein to explore the ensemble of conformations that ultimately foster large-scale conformational change is yet unknown. For this reason, we use THz spectroscopy to investigate the picosecond time scale binding modes and collective structural rearrangements that take place in hen egg white lysozyme (HEWL) when bound by the inhibitor (NAG)3. These protein thermal motions correspond to fluctuations that have a role in both selecting and sampling from the available protein intrinsic conformations that communicate function. Hence, investigation of these fast, collective modes may provide knowledge about the mechanism leading to the preferred binding process in HEWL-(NAG)3. Specifically, in this work we find that the picosecond time scale hydrogen-bonding rearrangements taking place in the protein hydration shell with binding modify the packing density within the hydrophobic core on a local level. These localized, intramolecular contact variations within the protein core appear to facilitate the large cooperative movements within the interfacial region separating the α- and β- domain that mediate binding. The THz time-scale fluctuations identified in the protein-ligand system may also reveal a molecular mechanism for substrate recognition.

Protein and water confined in nanometer-scale reverse micelles studied by near infrared, terahertz, and ultrafast visible spectroscopies

Protein-containing reverse (PCR) micelles are suitable systems to study the properties of proteins and waters in a cell-like environment. A model for determining the structural parameters of PCR micelles, such as the aqueous cavity size and molecule number of water within the reverse micelle, is presented. The model is based on an important hypothesis that the structural parameters of the protein-unfilled reverse micelle do not change after solubilization of protein. I describe a procedure using near infrared spectroscopy of OH stretching vibration band of water to verify the hypothesis. Further, the terahertz (THz) absorption spectrum of myoglobin is derived from THz time-domain spectroscopy of the PCR micellar solution, and the states of waters in reverse micelles with and without protein are discussed on the basis of the structural parameters. The last topic is on internal dynamics of PCR micelles on timescales from femtoseconds to nanoseconds studied by femtosecond time-resolved fluorescence spectroscopy.

Importance of protein conformational motions and electrostatic anchoring sites on the dynamics and hydrogen bond properties of hydration water

The microscopic dynamic properties of water molecules present in the vicinity of a protein are expected to be sensitive to its local conformational motions and the presence of polar and charged groups at the surface capable of anchoring water molecules through hydrogen bonds. In this work, we attempt to understand such sensitivity by performing detailed molecular dynamics simulations of the globular protein barstar solvated in aqueous medium. Our calculations demonstrate that enhanced confinement at the protein surface on freezing its local motions leads to increasingly restricted water mobility with long residence times around the secondary structures. It is found that the inability of the surface water molecules to bind with the protein residues by hydrogen bonds in the absence of protein-water (PW) electrostatic interactions is compensated by enhanced water-water hydrogen bonds around the protein with uniform bulklike behaviors. Importantly, it is further noticed that in contrast to the PW hydrogen bond relaxation time scale, the kinetics of the breaking and formation of such bonds are not affected on freezing the protein's conformational motions.

Thickness of the hydration layer of a protein from molecular dynamics simulation

Water molecules around a protein exhibit slow dynamics with respect to that of pure bulk water. One important issue in protein hydration is the thickness of the hydration layer (i.e., the distance from the protein surface up to which the water dynamics is influenced by the protein). Estimation of thickness is crucial to understand better the properties of "biological water" and the role that it plays in guiding the protein's function. We have performed an atomistic molecular dynamics simulation of an aqueous solution of the protein villin headpiece subdomain or HP-36 to estimate the thickness of its hydration water. In particular, several dynamical properties of water around different segments (three alpha-helices) of the protein have been calculated by varying the thickness of the hydration layers. It is found that in general the influence of the helices on water properties extends beyond the first hydration layer. However, the heterogeneous nature of water among the first hydration layers of the three helices diminishes as the thickness is increased. It indicates that, for a small protein such as HP-36, the thickness of "biological water" is uniform for different segments of the protein.

Effects of surface water on protein dynamics studied by a novel coarse-grained normal mode approach

Normal mode analysis (NMA) has received much attention as a direct approach to extract the collective motions of macromolecules. However, the stringent requirement of computational resources by classical all-atom NMA limits the size of the macromolecules to which the method is normally applied. We implemented a novel coarse-grained normal mode approach based on partitioning the all-atom Hessian matrix into relevant and nonrelevant parts. It is interesting to note that, using classical all-atom NMA results as a reference, we found that this method generates more accurate results than do other coarse-grained approaches, including elastic network model and block normal mode approaches. Moreover, this new method is effective in incorporating the energetic contributions from the nonrelevant atoms, including surface water molecules, into the coarse-grained protein motions. The importance of such improvements is demonstrated by the effect of surface water to shift vibrational modes to higher frequencies and by an increase in overlap of the coarse-grained eigenvector space (the motion directions) with that obtained from molecular dynamics simulations of solvated protein in a water box. These results not only confirm the quality of our method but also point out the importance of incorporating surface structural water in studying protein dynamics.

Applications of terahertz spectroscopy in biosystems

Terahertz (THz) spectroscopic investigations of condensed-phase biological samples are reviewed ranging from the simple crystalline forms of amino acids, carbohydrates and polypeptides to the more complex aqueous forms of small proteins, DNA and RNA. Vibrationally resolved studies of crystalline samples have revealed the exquisite sensitivity of THz modes to crystalline order, temperature, conformational form, peptide sequence and local solvate environment and have given unprecedented measures of the binding force constants and anharmonic character of the force fields, properties necessary to improve predictability but not readily obtainable using any other method. These studies have provided benchmark vibrational data on extended periodic structures for direct comparisons with classical (CHARMm) and quantum chemical (density functional theory) theories. For the larger amorphous and/or aqueous phase samples, the THz modes form a continuum-like absorption that arises because of the full accessibility to conformational space and/or the rapid time scale for inter-conversion in these environments. Despite severe absorption by liquid water, detailed investigations have uncovered the photo- and hydration-induced conformational flexibility of proteins, the solvent shell depth of the water/biomolecule boundary layers and the solvent reorientation dynamics occurring in these interfacial layers that occur on sub-picosecond time scales. As such, THz spectroscopy has enhanced and extended the accessibility to intermolecular forces, length- and timescales important in biological structure and activity.

Terahertz emission from collapsing field domains during switching of a gallium arsenide bipolar transistor

Broadband pulsed THz emission with peak power in the sub-mW range has been observed experimentally during avalanche switching in a gallium arsenide bipolar junction transistor at room temperature, while significantly higher total generated power is predicted in simulations. The emission is attributed to very fast oscillations in the conductivity current across the switching channels, which appear as a result of temporal evolution of the field domains generated in highly dense electron-hole plasma. This plasma is formed in turn by powerful impact ionization in multiple field domains of ultrahigh amplitude.