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

Molecular mechanisms for stabilizing biologics in the solid state

Protein drugs exhibit challenges of biophysical and biochemical instability due to their structural complexity and rich dynamics. Solid-state biologics aim to enhance stability by increasing molecular rigidity within the formulation matrix, representing a primary category of drug products alongside sterile liquid formulations. Understanding the molecular mechanisms behind the stabilization and destabilization of protein drugs, influenced by formulation composition and drying processes, provides scientific rationale for drug product design. This review aims to elaborate on the two primary models of water-to-sugar substitution and matrix vitrification, respectively, via thermodynamic and kinetic stabilization. It offers an up-to-date review of experimental investigations into these hypotheses, specifically elucidating protein structure and protein-excipient interactions at the molecular level, molecular dynamics across a broad range of motion regimes, and microscopic attributes such as protein-sugar and protein-salt miscibility and microenvironmental acidity, in relevant liquid, frozen, and solid states, using advanced biophysical techniques for solid-state analysis. Moreover, we discuss how these mechanistic understandings facilitate the investigation and prediction of critical stability behaviors and enables the design of solid biological drug products.

Structural features of the glassy state and their impact on the solid-state properties of organic molecules in pharmaceutical systems

This paper reviews the structure and properties of amorphous active pharmaceutical ingredients (APIs), including small molecules and proteins, in the glassy state (below the glass transition temperature, Tg). Amorphous materials in the neat state and formulated with excipients as miscible amorphous mixtures are included, and the role of absorbed water in affecting glass structure and stability has also been considered. We defined the term "structure" to indicate the way the various molecules in a glass interact with each other and form distinctive molecular arrangements as regions or domains of varying number of molecules, molecular packing, and density. Evidence is presented to suggest that such systems generally exist as heterogeneous structures made up of high-density domains surrounded by a lower density arrangement of molecules, termed the microstructure. It has been shown that the method of preparation and the time frame for handling and storage can give rise to variable glass structures and varying physical properties. Throughout this paper, examples are given of theoretical, computer simulation, and experimental studies which focus on the nature of intermolecular interactions, the size of heterogeneous higher density domains, and the impact of such systems on the relative physical and chemical stability of pharmaceutical systems.

Dynamics of Hydrogen Bonds between Water and Intrinsically Disordered and Structured Regions of Proteins

Recent studies indicate more restricted dynamics of water around intrinsically disordered proteins (IDPs) than structured proteins. We examine here the dynamics of hydrogen bonds between water molecules and two proteins, small ubiquitin-related modifier-1 (SUMO-1) and ubiquitin-conjugating enzyme E2I (UBC9), which we compare around intrinsically disordered regions (IDRs) and structured regions of these proteins. It has been recognized since some time that excluded volume effects, which influence access of water molecules to hydrogen-bonding sites, and the strength of hydrogen bonds between water and protein affect hydrogen bond lifetimes. While we find those two properties to mediate lifetimes of hydrogen bonds between water and protein residues in this study, we also find that the lifetimes are affected by the concentration of charged groups on other nearby residues. These factors are more important in determining the hydrogen bond lifetimes than whether a residue hydrogen bonding with water belongs to an IDR or to a structured region.

Quantifying the Effect of Guest Binding on Host Environment

The environment around a host-guest complex is defined by intermolecular interactions between the complex, solvent molecules, and counterions. These interactions govern both the solubility of these complexes and the rates of reactions occurring within the host molecules and can be critical to catalytic and separation applications of host-guest systems. However, these interactions are challenging to detect using standard analytical chemistry techniques. Here, we quantify the hydration and ion pairing of a FeII4L4 coordination cage with a set of guest molecules having widely varying physicochemical properties. The impact of guest properties on host ion pairing and hydration was determined through microwave microfluidic measurements paired with principal component analysis (PCA). This analysis showed that introducing guest molecules into solution displaced counterions that were bound to the cage, and that the solvent solubility of the guest has the greatest impact on the solvent and ion-pairing dynamics surrounding the host. Specifically, we found that when we performed PCA of the measured equivalent circuit parameters and the solubility and dipole moment, we observed a high (>90%) explained variance for the first two principal components for each circuit parameter. We also observed that cage-counterion pairing is well-described by a single ion-pairing type, with a one-step reaction model independent of the type of cargo, and that the ion-pairing association constant is reduced for cargo with higher water solubility. Quantifying hydration and cage-counterion interactions is a critical step to building the next generation of design criteria for host-guest chemistries.

Contrasting Changes in Strongly and Weakly Bound Hydration Water of a Protein upon Denaturation

Water is considered integral for the stabilization and function of proteins, which has recently attracted significant attention. However, the microscopic aspects of water ranging up to the second hydration shell, including strongly and weakly bound water at the sub-nanometer scale, are not yet well understood. Here, we combined terahertz spectroscopy, thermal measurements, and infrared spectroscopy to clarify how the strongly and weakly bound hydration water changes upon protein denaturation. With denaturation, that is, the exposure of hydrophobic groups in water and entanglement of hydrophilic groups, the number of strongly bound hydration water decreased, while the number of weakly bound hydration water increased. Even though the constraint of water due to hydrophobic hydration is weak, it extends to the second hydration shell as it is caused by the strengthening of hydrogen bonds between water molecules, which is likely the key microscopic mechanism for the destabilization of the native state due to hydration.

Terahertz spectroscopic monitoring and analysis of citrus leaf water status under low temperature stress

Low temperature stress, in the form of chilling and freezing, is one of the major environmental factors impacting on citrus yield, which changes plant's water state and results in the crops' sub-health or injury. The innovative terahertz (THz) spectroscopy and imaging based sensing technology has been shown to be a suitable tool for plant leaf water status determination, due to THz radiation's innate sensitivity to hydrogen bond vibration in aqueous solutions, which is usually related to plant phenotype change. We demonstrate experimentally that the THz absorption coefficient of leaf could be used for distinguishing plant's physiological stress status, exhibiting clear decreasing or increasing trend under chilling or freezing stress respectively. The underlying rationale might be that membrane damage shows a diverse pattern, changing the intra- or extra-cellular liquid environments, likely being linked to the various THz spectral characteristics. There were different adaptations in leaf morphology, leading to different leaf density, which in turn affects the water volume fraction. Moreover, different patterns of the dynamic equilibrium state of free water and bound water under chilling and freezing treatment were revealed by THz spectroscopy. Here, THz spectroscopic monitoring has shown unique potential for judging citrus's low temperature stress state through bio-water detection and discrimination.

Label-free laser spectroscopy for respiratory virus detection: A review

Infectious diseases are among the most severe threats to modern society. Current methods of virus infection detection based on genome tests need reagents and specialized laboratories. The desired characteristics of new virus detection methods are noninvasiveness, simplicity of implementation, real-time, low cost and label-free detection. There are two groups of methods for molecular biomarkers' detection and analysis: (i) a sample physical separation into individual molecular components and their identification, and (ii) sample content analysis by laser spectroscopy. Variations in the spectral data are typically minor. It requires the use of sophisticated analytical methods like machine learning. This review examines the current technological level of laser spectroscopy and machine learning methods in applications for virus infection detection.

Towards super-resolved terahertz microscopy for cellular imaging

Biomedical imaging includes the use of a variety of techniques to study organs and tissues. Some of the possible imaging modalities are more spread at clinical level (CT, MRI, PET), while others, such as light and electron microscopy are preferred in life sciences research. The choice of the imaging modalities can be based on the capability to study functional aspects of an organism, the delivered radiation dose to the patient, and the achievable resolution. In the last few decades, spectroscopists and imaging scientists have been interested in the use of terahertz (THz) frequencies (30 μm to 3 mm wavelength) due to the low photon energy associated (E∼1 meV, not causing breaking of the molecular bonds but still interacting with some vibrational modes) and the high penetration depth that is achievable. THz has been already adopted in security, quality control and material sciences. However, the adoption of THz frequencies for biological and clinical imaging means to face, as a major limitation, the very scarce resolution associated with the use of such long wavelengths. To address this aspect and reconcile the benefit of minimal harmfulness for bioimaging with the achievable resolving power, many attempts have been made. This review summarises the state-of-the-art of THz imaging applications aimed at achieving super-resolution, describing how practical aspects of optics and quasi-optics may be treated to efficaciously implement the use of THz as a new low-dose and versatile modality in biomedical imaging and clinical research.

Probing terahertz dynamics of multidomain protein in cell-like confinement

The development of meaningful descriptions of multidomain proteins exhibiting complex inter-domain dynamics modes is a key challenge for understanding their roles in molecular recognition and signalling processes. Here we developed a generally applicable approach for probing the low frequency collective hydration dynamics of multidomain proteins that uses terahertz spectroscopy of a protein molecule confined in a phospholipid reverse micelles environment (named Droplet THz). With the combination of normal mode analysis, we demonstrated the binding of calcium ions modulates the local inter-domain motion of the human coagulant factor VIII protein in a concentration-dependent manner. These findings highlight the Droplet THz as a valuable tool for dissecting the ultrafast dynamics of domain motion in the multidomain proteins and suggest a modulating mechanism of calcium ions on the structural flexibility and function of human coagulant proteins.

Roles of interfacial water states on advanced biomedical material design

When biomedical materials come into contact with body fluids, the first reaction that occurs on the material surface is hydration; proteins are then adsorbed and denatured on the hydrated material surface. The amount and degree of denaturation of adsorbed proteins affect subsequent cell behavior, including cell adhesion, migration, proliferation, and differentiation. Biomolecules are important for understanding the interactions and biological reactions of biomedical materials to elucidate the role of hydration in biomedical materials and their interaction partners. Analysis of the water states of hydrated materials is complicated and remains controversial; however, knowledge about interfacial water is useful for the design and development of advanced biomaterials. Herein, we summarize recent findings on the hydration of synthetic polymers, supramolecular materials, inorganic materials, proteins, and lipid membranes. Furthermore, we present recent advances in our understanding of the classification of interfacial water and advanced polymer biomaterials, based on the intermediate water concept.

Probing Adaptation of Hydration and Protein Dynamics to Temperature

Protein dynamics is strongly influenced by the surrounding environment and physiological conditions. Here we employ broadband megahertz-to-terahertz spectroscopy to explore the dynamics of water and myoglobin protein on an extended time scale from femto- to nanosecond. The dielectric spectra reveal several relaxations corresponding to the orientational polarization mechanism, including the dynamics of loosely bound, tightly bound, and bulk water, as well as collective vibrational modes of protein in an aqueous environment. The dynamics of loosely bound and bulk water follow non-Arrhenius behavior; however, the dynamics of water molecules in the tightly bound layer obeys the Arrhenius-type relation. Combining molecular simulations and effective-medium approximation, we have determined the number of water molecules in the tightly bound hydration layer and studied the dynamics of protein as a function of temperature. The results provide the important impact of water on the biochemical functions of proteins.

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.

A rationale for non-linear responses to strong electric fields in molecular dynamics simulations

Many approaches for calculation of the field-dependent electric properties of water solutions rely on the Onsager and Kirkwood theories of polar dielectrics. Such basic theories implicitly consider the electric field intensity to fulfill the so-called 'weak field conditions', i.e. to produce a linear response in the system. In this work we made use of molecular dynamics simulations to investigate possible non-linear effects induced by high intensity electric fields, specifically continuous wave bursts with nanosecond duration, comparing them with the ones predicted by the theory. We found that field intensities above 0.15 V nm^-1 produce remarkable nonlinear responses in the whole 100 MHz-100 GHz frequency window considered, with the onset of higher order polarization signals, which are the clear fingerprint of harmonic distorsions. That non-linear response turned out to depend on the considered frequency. We finally show that MD outcomes are consistent with a modelization based on an extended formulation of the Langevin function including a frequency-dependent parameter.

Enhanced Mobility during Diels-Alder Reaction: Results of Molecular Simulations

Recent measurements indicate enhanced mobility of solvent molecules during Diels-Alder (DA) and other common chemical reactions. We present results of molecular dynamics simulations of the last stages of the DA cycloaddition reaction, from the transition state configuration to product, of furfurylamine and maleimide in acetonitrile at reactant concentrations studied experimentally. We find enhanced mobility of solvent and reactant molecules up to at least a nanometer from the DA product over hundreds of picoseconds. Local heating is ruled out as a factor in the enhanced mobility observed in the simulations, which is instead found to be due to solvent relaxation following the formation of the DA product.

Quantification of solid-phase chemical reactions using the temperature-dependent terahertz pulsed spectroscopy, sum rule, and Arrhenius theory: thermal decomposition of α-lactose monohydrate

Transformations of the low-energy vibrational spectra are associated with structural changes in an analyte and closely related to the instability of weak chemical bounds. Terahertz (THz)/far-infrared optical spectroscopy is commonly used to probe such transformation, aimed at characterization of the underlying solid-phase chemical reactions in organic compounds. However, such studies usually provide quite qualitative information about the temperature- and time-dependent parameters of absorption peaks in dielectric spectra of an analyte. In this paper, an approach for quantitative analyses of the solid-phased chemical reactions based on the THz pulsed spectroscopy was developed. It involves studying an evolution of the sample optical properties, as a function of the analyte temperature and reaction time, and relies on the classical oscillator model, the sum rule, and the Arrhenius theory. The method allows one to determine the temperature-dependent reaction rate V₁(T) and activation energy E_a. To demonstrate the practical utility of this method, it was applied to study α-lactose monohydrate during its temperature-induced molecular decomposition. Analysis of the measured THz spectra revealed the increase of the reaction rate in the range of V₁ ≃ ~9 × 10^-4-10^-2 min^-1, when the analyte temperature rises from 313 to 393 K, while the Arrhenius activation energy is E_a ≃ ~45.4 kJ/mol. Thanks to a large number of obtained physical and chemical parameters, the developed approach expands capabilities of THz spectroscopy in chemical physics, analytical chemistry, and pharmaceutical industry.

The origin and impact of bound water around intrinsically disordered proteins

Proteins and water couple dynamically over a wide range of time scales. Motivated by their central role in protein function, protein-water dynamics and thermodynamics have been extensively studied for structured proteins, where correspondence to structural features has been made. However, properties controlling intrinsically disordered protein (IDP)-water dynamics are not yet known. We report results of megahertz-to-terahertz dielectric spectroscopy and molecular dynamics simulations of a group of IDPs with varying charge content along with structured proteins of similar size. Hydration water around IDPs is found to exhibit more heterogeneous rotational and translational dynamics compared with water around structured proteins of similar size, yielding on average more restricted dynamics around individual residues of IDPs, charged or neutral, compared with structured proteins. The on-average slower water dynamics is found to arise from excess tightly bound water in the first hydration layer, which is related to greater exposure to charged groups. The more tightly bound water to IDPs correlates with the smaller hydration shell found experimentally, and affects entropy associated with protein-water interactions, the contribution of which we estimate based on the dielectric measurements and simulations. Water-IDP dynamic coupling at terahertz frequencies is characterized by the dielectric measurements and simulations.

Newly identified C–H⋯O hydrogen bond in histidine

Histidine C–H bonds observed in protein structures include (clockwise from top left): myoglobin, β-lactamase, and photoactive yellow protein; calculations indicate that tautomeric/protonation state influences H-bonding ability (bottom left).

Long-range DNA-water interactions

DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ∼8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ∼18 Å from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.

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.

The Structural and Dynamical Properties of the Hydration of SNase Based on a Molecular Dynamics Simulation

The dynamics of protein–water fluctuations are of biological significance. Molecular dynamics simulations were performed in order to explore the hydration dynamics of staphylococcal nuclease (SNase) at different temperatures and mutation levels. A dynamical transition in hydration water (at ~210 K) can trigger larger-amplitude fluctuations of protein. The protein–water hydrogen bonds lost about 40% in the total change from 150 K to 210 K, while the Mean Square Displacement increased by little. The protein was activated when the hydration water in local had a comparable trend in making hydrogen bonds with protein– and other waters. The mutations changed the local chemical properties and the hydration exhibited a biphasic distribution, with two time scales. Hydrogen bonding relaxation governed the local protein fluctuations on the picosecond time scale, with the fastest time (24.9 ps) at the hydrophobic site and slowest time (40.4 ps) in the charged environment. The protein dynamic was related to the water’s translational diffusion via the relaxation of the protein–water’s H-bonding. The structural and dynamical properties of protein–water at the molecular level are fundamental to the physiological and functional mechanisms of SNase.

High precision dual-modulation differential terahertz ATR sensor for liquid measurements

We describe a highly sensitive and stable quantum-cascade-laser-based attenuated total reflection (ATR) terahertz sensor for the detection of very low concentration solutions, using a dual-modulation differential approach and ATR geometry. This sensor offers a very high dynamic range and a long-term stability of 40 dB, which extends the potential of terahertz radiation for the analysis of liquid and biological samples. The performance is illustrated by measurements on standard solutions of ions, sugars, and proteins, for concentrations down to 1 µM.

Fluorescent bacterial biosensor E. coli/pTdcR-TurboYFP sensitive to terahertz radiation

A fluorescent biosensor E. coli/pTdcR-TurboYFP sensitive to terahertz (THz) radiation was developed via transformation of Escherichia coli (E. coli) cells with plasmid, in which the promotor of the tdcR gene controls the expression of yellow fluorescent protein TurboYFP. The biosensor was exposed to THz radiation in various vessels and nutrient media. The threshold and dynamics of fluorescence were found to depend on irradiation conditions. Heat shock or chemical stress yielded the absence of fluorescence induction. The biosensor is applicable to studying influence of THz radiation on the activity of tdcR promotor that is involved in the transport and metabolism of threonine and serine in E. coli.

Hydration and dynamics of L-glutamate ion in aqueous solution

Aqueous solutions of sodium l-glutamate (NaGlu) in the concentration range 0 < c/M ≤ 1.90 at 25 °C were investigated by dielectric relaxation spectroscopy (DRS) and statistical mechanics (1D-RISM and 3D-RISM calculations) to study the hydration and dynamics of the l-glutamate (Glu-) anion. Although at c → 0 water molecules beyond the first hydration shell are dynamically affected, Glu- hydration is rather fragile and for c ⪆ 0.3 M apparently restricted to H2O molecules hydrogen bonding to the carboxylate groups. These hydrating dipoles are roughly parallel to the anion moment, leading to a significantly enhanced effective dipole moment of Glu-. However, l-glutamate dynamics is determined by the rotational diffusion of individual anions under hydrodynamic slip boundary conditions. Thus, the lifetime of the hydrate complexes, as well as of possibly formed [Na+Glu-]0 ionpairs and l-glutamate aggregates, cannot exceed the characteristic timescale for Glu- rotation.

Low-Frequency Protein Motions Coupled to Catalytic Sites

This review examines low-frequency vibrational modes of proteins and their coupling to enzyme catalytic sites. That protein motions are critical to enzyme function is clear, but the kinds of motions present in proteins and how they are involved in function remain unclear. Several models of enzyme-catalyzed reaction suggest that protein dynamics may be involved in the chemical step of the catalyzed reaction, but the evidence in support of such models is indirect. Spectroscopic studies of low-frequency protein vibrations consistently show that there are underdamped modes of the protein with frequencies in the tens of wavenumbers where overdamped behavior would be expected. Recent studies even show that such underdamped vibrations modulate enzyme active sites. These observations suggest that increasingly sophisticated spectroscopic methods will be able to unravel the link between low-frequency protein vibrations and enzyme function.

From adaptive resolution to molecular dynamics of open systems

ABSTRACT

We provide an overview of the Adaptive Resolution Simulation method (AdResS) based on discussing its basic principles and presenting its current numerical and theoretical developments. Examples of applications to systems of interest to soft matter, chemical physics, and condensed matter illustrate the method's advantages and limitations in its practical use and thus settle the challenge for further future numerical and theoretical developments.

Terahertz Spectroscopy Tracks Proteolysis by a Joint Analysis of Absorptance and Debye Model

Terahertz waves have attracted great attention in biomolecule research because of the fact that they cover the range of energy levels of weak interactions, skeleton vibrations, and dipole rotations during inter- and intramolecular interactions in biomacromolecules. In this study, we validated the feasibility of employing terahertz time-domain spectroscopy (THz-TDS) for the nondestructive and label-free monitoring of protein digestion. The acid protease, pepsin, was used at its optimal pH to hydrolyze bovine serum albumin. Correspondingly, the control group experiment was also conducted by adjusting the pH value to inactivate pepsin. The progress of these two experiments was tracked by a compact commercial THz-TDS for 1 h. On one hand, the reaction-time-dependent absorption coefficient was calculated, and a direct absorption coefficient analysis was completed. The results indicate that protein hydrolysis can be easily monitored over time by focusing on the variation tendency of the absorption coefficient from a macroscopic perspective. On the other hand, we explored the use of the Debye model to analyze the dielectric properties of the solution during protein hydrolysis. The results of the Debye analysis prove that it is possible to investigate in detail the microscopic dynamics of biomacromolecule solutions at the molecular level by THz-TDS. Our research examined the process of protein hydrolysis by a combination of absorption spectra and Debye analysis and demonstrated that terahertz spectroscopy is a powerful technology for the investigation of biomolecular reactions, with potential applications in variety of fields.

Terahertz pulse-altered gene networks in human induced pluripotent stem cells

Terahertz (THz) irradiation has been exploited in biomedical applications involving non-invasive manipulation of living cells. We developed an apparatus for studying the effects of THz pulse irradiation on living human induced pluripotent stem cells. The THz pulse of the maximum electric field reached 0.5 MV/cm and was applied for one hour with 1 kHz repetition to the entire cell-culture area, a diameter of 1 mm. RNA sequencing of global gene-expression revealed that many THz-regulated genes were driven by zinc-finger transcription factors. Combined with a consideration of the interactions of metal ions and a THz electric field, these results imply that the local intracellular concentration of metal ions, such as Zn²⁺, was changed by the effective electrical force of our THz pulse.

The effect of the degree of substitution on the solubility of cellulose acetoacetates in water: A molecular dynamics simulation and density functional theory study

The effect of the degree of substitution (DS) on the aqueous solubility of cellulose acetoacetates (CAA) was investigated by molecular dynamics simulations and density functional theory calculations. Using average non-covalent interaction and the electrostatic potential analyses done on cellobiose as the model, it was showed both polar and non-polar areas of the system increased as the more hydroxyls were replaced by acetoacetate groups. Analyses of the solvation free energies of a celludecose (glucan containing 10 monosaccharide sugar units) at constant pressure and temperature showed the polar solvation free energies and the number of decose-water hydrogen bonds increased as DS was varied from 0.3 to 0.8, which contributes to higher solubility in water. When the DS of CAA increased from 0.8 to 1.5, it became insoluble again because the plateaued increase in solvation free energy could no longer compensate for the decreasing CAA-water hydrogen bonding interactions. The growing van der Waals interactions among CAA molecules as the molecule grows bigger with each attached AA group also contributes to the decreasing water solubility.

Entropic Nature of the Debye Relaxation in Glass-Forming Monoalcohols

The dynamics and thermodynamics of the Debye and structural (α) relaxations in isomeric monoalcohols near the glass transition temperature T_g are explored using dielectric and calorimetric techniques. The α relaxation strength at T_g is found to correlate with the heat capacity increment, but no thermal signals can be detected to link to the Debye relaxation. We also observed that the activation energy of the Debye relaxation in monoalcohols is quantitatively correlated with that of the α relaxation at the kinetic T_g, sharing the dynamic behavior of the Rouse modes found in polymers. The experimental results together with the analogy to the Rouse modes in polymers suggest that the Debye process in monoalcohols is an entropic process manifested by the total dipole fluctuation of the supramolecular structures, which is triggered and driven by the α relaxation.

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.

Water Structure and Properties at Hydrophilic and Hydrophobic Surfaces

The properties of water on both molecular and macroscopic surfaces critically influence a wide range of physical behaviors, with applications spanning from membrane science to catalysis to protein engineering. Yet, our current understanding of water interfacing molecular and material surfaces is incomplete, in part because measurement of water structure and molecular-scale properties challenges even the most advanced experimental characterization techniques and computational approaches. This review highlights progress in the ongoing development of tools working to answer fundamental questions on the principles that govern the interactions between water and surfaces. One outstanding and critical question is what universal molecular signatures capture the hydrophobicity of different surfaces in an operationally meaningful way, since traditional macroscopic hydrophobicity measures like contact angles fail to capture even basic properties of molecular or extended surfaces with any heterogeneity at the nanometer length scale. Resolving this grand challenge will require close interactions between state-of-the-art experiments, simulations, and theory, spanning research groups and using agreed-upon model systems, to synthesize an integrated knowledge of solvation water structure, dynamics, and thermodynamics.

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.

Protein Solvent Shell Structure Provides Rapid Analysis of Hydration Dynamics

The solvation layer surrounding a protein is clearly an intrinsic part of protein structure-dynamics-function, and our understanding of how the hydration dynamics influences protein function is emerging. We have recently reported simulations indicating a correlation between regional hydration dynamics and the structure of the solvation layer around different regions of the enzyme Candida antarctica lipase B, wherein the radial distribution function (RDF) was used to calculate the pairwise entropy, providing a link between dynamics (diffusion) and thermodynamics (excess entropy) known as Rosenfeld scaling. Regions with higher RDF values/peaks in the hydration layer (the first peak, within 6 Å of the protein surface) have faster diffusion in the hydration layer. The finding thus hinted at a handle for rapid evaluation of hydration dynamics at different regions on the protein surface in molecular dynamics simulations. Such an approach may move the analysis of hydration dynamics from a specialized venture to routine analysis, enabling an informatics approach to evaluate the role of hydration dynamics in biomolecular function. This paper first confirms that the correlation between regional diffusive dynamics and hydration layer structure (via water center of mass around protein side-chain atom RDF) is observed as a general relationship across a set of proteins. Second, it seeks to devise an approach for rapid analysis of hydration dynamics, determining the minimum amount of information and computational effort required to get a reliable value of hydration dynamics from structural data in MD simulations based on the protein-water RDF. A linear regression model using the integral of the hydration layer in the water-protein RDF was found to provide statistically equivalent apparent diffusion coefficients at the 95% confidence level for a set of 92 regions within five different proteins. In summary, RDF analysis of 10 ns of data after simulation convergence is sufficient to accurately map regions of fast and slow hydration dynamics around a protein surface. Additionally, it is anticipated that a quick look at protein-water RDFs, comparing peak heights, will be useful to provide a qualitative ranking of regions of faster and slower hydration dynamics at the protein surface for rapid analysis when investigating the role of solvent dynamics in protein function.

How Does Solvation Layer Mobility Affect Protein Structural Dynamics?

Solvation is critical for protein structural dynamics. Spectroscopic studies have indicated relationships between protein and solvent dynamics, and rates of gas binding to heme proteins in aqueous solution were previously observed to depend inversely on solution viscosity. In this work, the solvent-compatible enzyme Candida antarctica lipase B, which functions in aqueous and organic solvents, was modeled using molecular dynamics simulations. Data was obtained for the enzyme in acetonitrile, cyclohexane, n-butanol, and tert-butanol, in addition to water. Protein dynamics and solvation shell dynamics are characterized regionally: for each α-helix, β-sheet, and loop or connector region. Correlations are seen between solvent mobility and protein flexibility. So, does local viscosity explain the relationship between protein structural dynamics and solvation layer dynamics? Halle and Davidovic presented a cogent analysis of data describing the global hydrodynamics of a protein (tumbling in solution) that fits a model in which the protein's interfacial viscosity is higher than that of bulk water's, due to retarded water dynamics in the hydration layer (measured in NMR τ2 reorientation times). Numerous experiments have shown coupling between protein and solvation layer dynamics in site-specific measurements. Our data provides spatially-resolved characterization of solvent shell dynamics, showing correlations between regional solvation layer dynamics and protein dynamics in both aqueous and organic solvents. Correlations between protein flexibility and inverse solvent viscosity (1/η) are considered across several protein regions and for a rather disparate collection of solvents. It is seen that the correlation is consistently higher when local solvent shell dynamics are considered, rather than bulk viscosity. Protein flexibility is seen to correlate best with either the local interfacial viscosity or the ratio of the mobility of an organic solvent in a regional solvation layer relative to hydration dynamics around the same region. Results provide insight into the function of aqueous proteins, while also suggesting a framework for interpreting and predicting enzyme structural dynamics in non-aqueous solvents, based on the mobility of solvents within the solvation layer. We suggest that Kramers' theory may be used in future work to model protein conformational transitions in different solvents by incorporating local viscosity effects.