Influence of water and trehalose on α- and β-relaxation of freeze-dried lysozyme formulations

Molecular mobility in form of alpha and beta relaxations is considered crucial for characterization of amorphous lyophilizates and reflected in the transition temperatures Tgα and Tgβ. Based on an overview of applied methods to study beta relaxations, Dynamic Mechanical analysis was used to measure Tgα and Tgβ in amorphous freeze-dried samples. Lysozyme and trehalose as well as their mixtures in varying ratios were investigated. Three different residual moisture levels, ranging from roughly 0.5-7 % (w/w), were prepared via equilibration of the freeze-dried samples. Known plasticising effects of water on Tgα were confirmed, also via differential scanning calorimetry. In addition and contrary to expectations, an influence of water on the Tgβ also was observed. On the other hand, an increasing amount of trehalose lowered Tgα but increased Tgβ showing that Tgα and Tgβ are not paired. The findings were interpreted with regard to their underlying molecular mechanisms and a correlation with the known influences of water and trehalose on stability. The results provide encouraging hints for future stability studies of freeze-dried protein formulations, which are urgently needed, not least for reasons of sustainability.

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.

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.

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.

Similarities and Differences among Protein Dynamics Studied by Variable Temperature Nuclear Magnetic Resonance Relaxation

Understanding and describing the dynamics of proteins is one of the major challenges in biology. Here, we use multifield variable-temperature NMR longitudinal relaxation (R₁) measurements to determine the hierarchical activation energies of motions of four different proteins: two small globular proteins (GB1 and the SH3 domain of α-spectrin), an intrinsically disordered protein (the C-terminus of the nucleoprotein of the Sendai virus, Sendai Ntail), and an outer membrane protein (OmpG). The activation energies map the motions occurring in the side chains, in the backbone, and in the hydration shells of the proteins. We were able to identify similarities and differences in the average motions of the proteins. We find that the NMR relaxation properties of the four proteins do share similar features. The data characterizing average backbone motions are found to be very similar, the same for methyl group rotations, and similar activation energies are measured. The main observed difference occurs for the intrinsically disordered Sendai Ntail, where we observe much lower energy of activation for motions of protons associated with the protein-solvent interface as compared to the others. We also observe variability between the proteins regarding side chain ¹⁵N relaxation of lysine residues, with a higher activation energy observed in OmpG. This hints at strong interactions with negatively charged lipids in the bilayer and provides a possible mechanistic clue for the "positive-inside" rule for helical membrane proteins. Overall, these observations refine the understanding of the similarities and differences between hierarchical dynamics in proteins.

Hydrophobic Collapse of Ubiquitin Generates Rapid Protein-Water Motions

We report time-resolved measurements of the coupled protein-water modes of solvated ubiquitin during protein folding. Kinetic terahertz absorption (KITA) spectroscopy serves as a label-free technique for monitoring large scale conformational changes and folding of proteins subsequent to a sudden T-jump. We report here KITA measurements at an unprecedented time resolution of 500 ns, a resolution 2 orders of magnitude better than those of any previous KITA measurements, which reveal the coupled ubiquitin-solvent dynamics even in the initial phase of hydrophobic collapse. Complementary equilibrium experiments and molecular simulations of ubiquitin solutions are performed to clarify non-equilibrium contributions and reveal the molecular picture upon a change in structure, respectively. On the basis of our results, we propose that in the case of ubiquitin a rapid (<500 ns) initial phase of the hydrophobic collapse from the elongated protein to a molten globule structure precedes secondary structure formation. We find that these very first steps, including large-amplitude changes within the unfolded manifold, are accompanied by a rapid (<500 ns) pronounced change of the coupled protein-solvent response. The KITA response upon secondary structure formation exhibits an opposite sign, which indicates a distinct effect on the solvent-exposed surface.

Aqueous TMAO solutions as seen by theoretical THz spectroscopy: hydrophilic versus hydrophobic water

All THz resonances of aqueous TMAO solutions are computed and assigned based on ab initio molecular dynamics simulations.

Thermal Gradient Mid- and Far-Infrared Spectroscopy as Tools for Characterization of Protein Carbohydrate Lyophilizates

Protein drugs play an important role in modern day medicine. Typically, these proteins are formulated as liquids requiring cold chain processing. To circumvent the cold chain and achieve better storage stability, these proteins can be dried in the presence of carbohydrates. We demonstrate that thermal gradient mid- and far-infrared spectroscopy (FTIR and THz-TDS, respectively) can provide useful information about solid-state protein carbohydrate formulations regarding mobility and intermolecular interactions. A model protein (BSA) was lyophilized in the presence of three carbohydrates with different size and protein stabilizing capacity. A gradual increase in mobility was observed with increasing temperature in formulations containing protein and/or larger carbohydrates (oligo- or polysaccharides), lacking a clear onset of fast mobility as was observed for smaller molecules. Furthermore, both techniques are able to identify the glass transition temperatures (T_g) of the samples. FTIR provides additional information as it can independently monitor changes in protein and carbohydrate bands at the T_g. Lastly, THz-TDS confirms previous findings that protein–carbohydrate interactions decrease with increasing molecular weight of the carbohydrate, which results in decreased protein stabilization.

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

We have performed dielectric spectral measurements of lysozyme in a solid state to understand the effects of hydration and thermal excitation on the low-frequency dynamics of protein. Dielectric measurements were performed under changing hydration conditions at room temperature in the frequency region of 0.5 GHz to 1.8 THz. We also studied the temperature dependence (83 to 293 K) of the complex dielectric spectra in the THz frequency region (0.3 THz to 1.8 THz). Spectral analyses were performed using model functions for the complex dielectric constant. To reproduce the spectra, we found that two relaxational modes and two underdamped modes are necessary together with an ionic conductivity term in the model function. At room temperature, the two relaxational modes have relaxation times of ∼20 ps and ∼100 ps. The faster component has a major spectral intensity and is suggested to be due to coupled water-protein motion. The two underdamped modes are necessary to reproduce the temperature dependence of the spectra in the THz region satisfactorily. The protein dynamical transition is a well-known behavior in the neutron-scattering experiment for proteins, where the atomic mean-square displacement shows a sudden change in the temperature dependence at approximately 200 K, when the samples are hydrated. A similar behavior has also been observed in the temperature dependence of the absorption spectra of protein in the THz frequency region. From our broadband dielectric spectroscopic measurements, we conclude that the increase in the spectral intensities in the THz region at approximately 200 K is due to a spectral blue-shift of the fast relaxational mode.

Direct measurement of molecular mobility and crystallisation of amorphous pharmaceuticals using terahertz spectroscopy

Despite much effort in the area, no comprehensive understanding of the formation and behaviour of amorphous solids has yet been achieved. This severely limits the industrial application of such materials, including drug delivery where, in principle, amorphous solids have demonstrated their great usefulness in increasing the bioavailability of poorly aqueous soluble active pharmaceutical ingredients. Terahertz time-domain spectroscopy is a relatively novel analytical technique that can be used to measure the fast molecular dynamics of molecules with high accuracy in a non-contact and non-destructive fashion. Over the past decade a number of applications for the characterisation of amorphous drug molecules and formulations have been developed and it has been demonstrated how this technique can be used to determine the onset and strength in molecular mobility that underpins the crystallisation of amorphous drugs. In this review we provide an overview of the history, fundamentals and future perspective of pharmaceutical applications related to the terahertz dynamics of amorphous systems.

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.

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.

Protein surface and core dynamics show concerted hydration-dependent activation

By specifically labeling leucine/valine methyl groups and lysine side chains "inside" and "outside" dynamics of proteins on the nanosecond timescale are compared using neutron scattering. Surprisingly, both groups display similar dynamics as a function of temperature, and the buried hydrophobic core is sensitive to hydration and undergoes a dynamical transition.

Computer Simulation of Water Sorption on Flexible Protein Crystals

The first simulation study of water sorption on a flexible protein crystal is presented, along with a new computational approach for calculating sorption isotherms on compliant materials. The flexible ubiquitin crystal examined in the study exhibits appreciable sorption-induced swelling during fluid uptake, similar to that reported in experiments on protein powders. A completely rigid ubiquitin crystal is also examined to investigate the impact that this swelling behavior has on water sorption. The water isotherms for the flexible crystal exhibit Type II-like behavior with sorption hysteresis, which is consistent with experimental measurements on protein powders. Both of these behaviors, however, are absent in the rigid crystal, indicating that modeling flexibility is crucial for predicting water sorption behavior in protein systems. Changes in the enthalpy of adsorption, specific volume, and internal protein fluctuations that occur during sorption in the flexible crystal are also shown to compare favorably with experiment.