Description of Hydration Water in Protein (Green Fluorescent Protein) Solution

Spontaneous Curvature Induction in an Artificial Bilayer Membrane

Maintaining lipid asymmetry across membrane leaflets is critical for functions like vesicular traffic and organelle homeostasis. However, a lack of molecular-level understanding of the mechanisms underlying membrane fission and fusion processes in synthetic systems precludes their development as artificial analogs. Here, we report asymmetry induction of a bilayer membrane formed by an extended π-conjugated molecule with oxyalkylene side chains bearing terminal tertiary amine moieties (BA1) in water. Autogenous protonation of the tertiary amines in the periphery of the bilayer by water induces anisotropic curvature, resulting in membrane fission to form vesicles and can be monitored using time-dependent spectroscopy and microscopy. Interestingly, upon loss of the induced asymmetry by extensive protonation using an organic acid restored bilayer membrane. The mechanism leading to the compositional asymmetry in the leaflet and curvature induction in the membrane is validated by density functional theory (DFT) calculations. Studies extended to control molecules having changes in hydrophilic (BA2) and hydrophobic (BA3) segments provide insight into the delicate nature of molecular scale interactions in the dynamic transformation of supramolecular structures. The synergic effect of hydrophobic interaction and the hydrated state of BA1 aggregates provide dynamicity and unusual stability. Our study unveils mechanistic insight into the dynamic transformation of bilayer membranes into vesicles.

An experimentally representative in-silico protocol for dynamical studies of lyophilised and weakly hydrated amorphous proteins

Characterization of biopolymers in both dry and weakly hydrated amorphous states has implications for the pharmaceutical industry since it provides understanding of the effect of lyophilisation on stability and biological activity. Atomistic Molecular Dynamics (MD) simulations probe structural and dynamical features related to system functionality. However, while simulations in homogenous aqueous environments are routine, dehydrated model assemblies are a challenge with systems investigated in-silico needing careful consideration; simulated systems potentially differing markedly despite seemingly negligible changes in procedure. Here we propose an in-silico protocol to model proteins in lyophilised and weakly hydrated amorphous states that is both more experimentally representative and routinely applicable. Since the outputs from MD align directly with those accessed by neutron scattering, the efficacy of the simulation protocol proposed is shown by validating against experimental neutron data for apoferritin and insulin. This work also highlights that without cooperative experimental and simulative data, development of simulative procedures using MD alone would prove most challenging.

Effects of Crowding and Cosolutes on Biomolecular Function at Extreme Environmental Conditions

The extent of the effect of cellular crowding and cosolutes on the functioning of proteins and cells is manifold and includes the stabilization of the biomolecular systems, the excluded volume effect, and the modulation of molecular dynamics. Simultaneously, it is becoming increasingly clear how important it is to take the environment into account if we are to shed light on biological function under various external conditions. Many biosystems thrive under extreme conditions, including the deep sea and subseafloor crust, and can take advantage of some of the effects of crowding. These relationships have been studied in recent years using various biophysical techniques, including neutron and X-ray scattering, calorimetry, FTIR, UV-vis and fluorescence spectroscopies. Combining knowledge of the structure and conformational dynamics of biomolecules under extreme conditions, such as temperature, high hydrostatic pressure, and high salinity, we highlight the importance of considering all results in the context of the environment. Here we discuss crowding and cosolute effects on proteins, nucleic acids, membranes, and live cells and explain how it is possible to experimentally separate crowding-induced effects from other influences. Such findings will contribute to a better understanding of the homeoviscous adaptation of organisms and the limits of life in general.

Variations in the Protein Hydration and Hydrogen-Bond Network of Water Molecules Induced by the Changes in the Secondary Structures of Proteins Studied through Near-Infrared Spectroscopy

This study investigated how the secondary structural changes of proteins in aqueous solutions affect their hydration and the hydrogen-bond network of water molecules using near-infrared (NIR) spectroscopy. The aqueous solutions of three types of proteins, i.e., ovalbumin, β-lactoglobulin, and bovine serum albumin, were denatured by heating, and changes in the NIR bands of water reflecting the states of hydrogen bonds induced via protein secondary structural changes were investigated. On heating, the intermolecular hydrogen bonds between water molecules as well as between water and protein molecules were broken, and protein molecules were no longer strongly bound by the surrounding water molecules. Consequently, the denaturation was observed to proceed depending on the thermodynamic properties of the proteins. When the aqueous solutions of proteins were cooled after denaturation, the hydrogen-bond network was reformed. However, the state of protein hydration was changed owing to the secondary structural changes of proteins, and the variation patterns were different depending on the protein species. These changes in protein hydration may be derived from the differences in the surface charges of proteins. The elucidation of the mechanism of protein hydration and the formation of the hydrogen-bond network of water molecules will afford a comprehensive understanding of the protein functioning and dysfunctioning derived from the structural changes in proteins.

New insights into the protein stabilizing effects of trehalose by comparing with sucrose

Disaccharides are well known to be efficient stabilizers of proteins, for example in the case of lyophilization or cryopreservation. However, although all disaccharides seem to exhibit bioprotective and stabilizing properties, it is clear that trehalose is generally superior compared to other disaccharides. The aim of this study was to understand this by comparing how the structural and dynamical properties of aqueous trehalose and sucrose solutions influence the protein myoglobin (Mb). The structural studies were based on neutron and X-ray diffraction in combination with empirical potential structure refinement (EPSR) modeling, whereas the dynamical studies were based on quasielastic neutron scattering (QENS) and molecular dynamics (MD) simulations. The results show that the overall differences in the structure and dynamics of the two systems are small, but nevertheless there are some important differences which may explain the superior stabilizing effects of trehalose. It was found that in both systems the protein is preferentially hydrated by water, but that this effect is more pronounced for trehalose, i.e. trehalose forms less hydrogen bonds to the protein surface than sucrose. Furthermore, the rotational motion around dihedrals between the two glucose rings of trehalose is slower than in the case of the dihedrals between the glucose and fructose rings of sucrose. This leads to a less perturbed protein structure in the case of trehalose. The observations indicate that an aqueous environment closest to the protein molecules is beneficial for an efficient bioprotective solution.

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.

Spatially confined protein assembly in hierarchical mesoporous metal-organic framework

Immobilization of biomolecules into porous materials could lead to significantly enhanced performance in terms of stability towards harsh reaction conditions and easier separation for their reuse. Metal-Organic Frameworks (MOFs), offering unique structural features, have emerged as a promising platform for immobilizing large biomolecules. Although many indirect methods have been used to investigate the immobilized biomolecules for diverse applications, understanding their spatial arrangement in the pores of MOFs is still preliminary due to the difficulties in directly monitoring their conformations. To gain insights into the spatial arrangement of biomolecules within the nanopores. We used in situ small-angle neutron scattering (SANS) to probe deuterated green fluorescent protein (d-GFP) entrapped in a mesoporous MOF. Our work revealed that GFP molecules are spatially arranged in adjacent nanosized cavities of MOF-919 to form "assembly" through adsorbate-adsorbate interactions across pore apertures. Our findings, therefore, lay a crucial foundation for the identification of proteins structural basics under confinement environment of MOFs.

Dynamic Electric Field Alignment Determines the Water Rotational Motion around Protein

Water rotational dynamics in biomolecular solution is crucial to evaluating and controlling biomolecule stability. In this molecular dynamics simulation (MD) study on lysozyme solutions, we present how the exerted internal electric field determines water rotational dynamics. We find that the relaxation time of water rotation is equivalent to that of the reorientation of the exerted overall electric field for every single water molecule, regardless of its translation mode. Namely, water molecular rotation synchronizes with the exerted field reorientation. We also map the reorientation process of the electric field at fixed points relative to protein in the solution, which displays the local hydration dynamics commensurate with the reported time-dependent fluorescence Stokes shift (TDFSS) measurements. Comparing the spatial distribution of local field reorientation relaxation time with that of rotational relaxation time, we further suggest that water rotation dynamics are subject to the reorientation of the local overall field within the hydration layer. While outside the hydration layer, the relaxation time of the local electric field reorientation is short enough (subpicosecond) to assume the δ function, showing the electric force with randomly changing orientation is applied to each water molecule.

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.

CHESS: The future direct geometry spectrometer at the second target station

CHESS, chopper spectrometer examining small samples, is a planned direct geometry neutron chopper spectrometer designed to detect and analyze weak signals intrinsic to small cross sections (e.g., small mass, small magnetic moments, or neutron absorbing materials) in powders, liquids, and crystals. CHESS is optimized to enable transformative investigations of quantum materials, spin liquids, thermoelectrics, battery materials, and liquids. The broad dynamic range of the instrument is also well suited to study relaxation processes and excitations in soft and biological matter. The 15 Hz repetition rate of the Second Target Station at the Spallation Neutron Source enables the use of multiple incident energies within a single source pulse, greatly expanding the information gained in a single measurement. Furthermore, the high flux grants an enhanced capability for polarization analysis. This enables the separation of nuclear from magnetic scattering or coherent from incoherent scattering in hydrogenous materials over a large range of energy and momentum transfer. This paper presents optimizations and technical solutions to address the key requirements envisioned in the science case and the anticipated uses of this instrument.

Diffusions in Aqueous Solutions with Multivalent Cations and Especially in Cationic First Hydration Shell

Compared with univalent cationic diffusion, little is known about the diffusion behavior of multivalent cations let alone the diffusion of water in their first hydration shell. Here, we show that all published translational diffusion coefficients of multivalent cations and water measured at room temperature exhibit the same concentration dependence when plotted as a function of the mass fraction of free water or of hydrated solute. This behavior is held until only hydration and confined water remain in solutions, wherein their concentration dependences become cationic charge number-dependent. We also found that the iceberg model can well describe the diffusions of multivalent cations with decreasing water content until only hydration water is present. However, ¹H-pulsed-field-gradient nuclear magnetic resonance measurements confirmed that ¹H in the first hydration shell diffuses faster than Al³⁺ at room temperature and they have the same diffusion coefficient only with decreasing the temperature down to about 243 K. Therefore, iceberg model only equivalently describes the effect of strong ion-water interaction on multivalent cationic diffusion. These results will also help us reconceive our current understanding of the pathway for hydration water affecting the diffusion behavior of solutes with relatively weak solute-water interactions.

Hydration Dynamics of Model Peptides with Different Hydrophobic Character

The multi-scale dynamics of aqueous solutions of the hydrophilic peptide N-acetyl-glycine-methylamide (NAGMA) have been investigated through extended frequency-range depolarized light scattering (EDLS), which enables the broad-band detection of collective polarizability anisotropy fluctuations. The results have been compared to those obtained for N-acetyl-leucinemethylamide (NALMA), an amphiphilic peptide which shares with NAGMA the same polar backbone, but also contains an apolar group. Our study indicates that the two model peptides induce similar effects on the fast translational dynamics of surrounding water. Both systems slow down the mobility of solvating water molecules by a factor 6–8, with respect to the bulk. Moreover, the two peptides cause a comparable far-reaching spatial perturbation extending to more than two hydration layers in diluted conditions. The observed concentration dependence of the hydration number is explained considering the random superposition of different hydration shells, while no indication of solute aggregation phenomena has been found. The results indicate that the effect on the dynamics of water solvating the amphiphilic peptide is dominated by the hydrophilic backbone. The minor impact of the hydrophobic moiety on hydration features is consistent with structural findings derived by Fourier transform infrared (FTIR) measurements, performed in attenuated total reflectance (ATR) configuration. Additionally, we give evidence that, for both systems, the relaxation mode in the GHz frequency range probed by EDLS is related to solute rotational dynamics. The rotation of NALMA occurs at higher timescales, with respect to the rotation of NAGMA; both processes are significantly slower than the structural dynamics of hydration water, suggesting that solute and solvent motions are uncoupled. Finally, our results do not indicate the presence of super-slow water (relaxation times in the order of tens of picoseconds) around the peptides investigated.

Influence of hydration on segmental chain librations and dynamical transition in lipid bilayers

Continuous wave electron paramagnetic resonance spectroscopy of chain-labeled phospholipids is used to investigate the effects of hydration on the librational oscillations and the dynamical transition of phospholipid membranes in the low-temperature range 120-270 K. Bilayers of dipalmitoylphostatidiycholine (DPPC) spin-labeled at the first acyl chain segments and at the methyl ends and prepared at full, low, and very low hydration are considered. The segmental mean-square angular amplitudes of librations, 〈α²〉, are larger in the bilayer interior than at the polar/apolar interface and larger in the fully and low hydrated than in the very low hydrated membranes. For chain segments at the beginning of the hydrocarbon region, 〈α²〉-values are markedly restricted and temperature independent in DPPC with the lowest water content, whereas they increase with temperature in the low and fully hydrated bilayers, particularly at the highest temperatures. For chain segments at the chain termini, the librational amplitudes increase progressively, first slowly and then more rapidly with temperature in bilayers at any level of hydration. From the temperature dependence of the mean-square librational amplitude, the dynamical transition is detected around 240 K at the polar/apolar interface in fully and low hydrated DPPC and at around 225 K at the inner hydrocarbon region for bilayers at any hydration condition. At the dynamical transition the bilayers cross low energy barriers of activation energy in the range 10-20 kJ/mol. The results highlight biophysical properties of DPPC bilayers at low-temperature and provide evidence of the effects of the hydration on the dynamical transition in bilayers.

Self-Diffusive Properties of the Intrinsically Disordered Protein Histatin 5 and the Impact of Crowding Thereon: A Combined Neutron Spectroscopy and Molecular Dynamics Simulation Study

Intrinsically disordered proteins (IDPs) are proteins that, in comparison with globular/structured proteins, lack a distinct tertiary structure. Here, we use the model IDP, Histatin 5, for studying its dynamical properties under self-crowding conditions with quasi-elastic neutron scattering in combination with full atomistic molecular dynamics (MD) simulations. The aim is to determine the effects of crowding on the center-of-mass diffusion as well as the internal diffusive behavior. The diffusion was found to decrease significantly, which we hypothesize can be attributed to some degree of aggregation at higher protein concentrations, (≥100 mg/mL), as indicated by recent small-angle X-ray scattering studies. Temperature effects are also considered and found to, largely, follow Stokes–Einstein behavior. Simple geometric considerations fail to accurately predict the rates of diffusion, while simulations show semiquantitative agreement with experiments, dependent on assumptions of the ratio between translational and rotational diffusion. A scaling law that previously was found to successfully describe the behavior of globular proteins was found to be inadequate for the IDP, Histatin 5. Analysis of the MD simulations show that the width of the distribution with respect to diffusion is not a simplistic mirroring of the distribution of radius of gyration, hence, displaying the particular features of IDPs that need to be accounted for.

Site-Specific Water Dynamics in the First Hydration Layer of an Anti-Freeze Glyco-Protein: A Simulation Study

Antifreeze glycoproteins (AFGPs) inhibit ice recrystallization by a mechanism remaining largely elusive. Dynamics of AFGPs’ hydration water and its involvement in the antifreeze activity, for instance, have not been identified...

Universal dynamical onset in water at distinct material interfaces

Interfacial water remains liquid and mobile much below 0 °C, imparting flexibility to the encapsulated materials to ensure their diverse functions at subzero temperatures. However, a united picture that can describe the dynamical differences of interfacial water on different materials and its role in imparting system-specific flexibility to distinct materials is lacking. By combining neutron spectroscopy and isotope labeling, we explored the dynamics of water and the underlying substrates independently below 0 °C across a broad range of materials. Surprisingly, while the function-related anharmonic dynamical onset in the materials exhibits diverse activation temperatures, the surface water presents a universal onset at a common temperature. Further analysis of the neutron experiment and simulation results revealed that the universal onset of water results from an intrinsic surface-independent relaxation: switching of hydrogen bonds between neighboring water molecules with a common energy barrier of ∼35 kJ mol⁻¹.

We demonstrated that the dynamical onset of interfacial water is an intrinsic property of water itself, resulting from a surface independent relaxation process in water with an approximately universal energy barrier of ∼35 kJ mol⁻¹.

Broadband Wide-Angle VElocity Selector (BWAVES) neutron spectrometer designed for the SNS Second Target Station

A recently proposed wide-angle velocity selector (WAVES) device for choosing the velocity of detected neutrons after they have been scattered by the sample paves the way for inverted geometry neutron spectrometers with continuously adjustable final neutron wavelength. BWAVES broadband inverted geometry spectrometer proposed for the Second Target Station at the Spallation Neutron Source at Oak Ridge National Laboratory is designed using WAVES to simultaneously probe dynamic processes spanning 4.5 decades in time (energy transfer). This makes BWAVES a uniquely flexible instrument which can be viewed as either a quasielasitc neutron scattering (QENS) spectrometer with a practically unlimited (overlapping with the vibrational excitations) range of energy transfers, or a broadband inelastic vibrational neutron spectrometer with QENS capabilities, including a range of accessible momentum transfer (Q) and a sufficiently high energy resolution at the elastic line. The new capabilities offered by BWAVES will expand the application of neutron scattering in ways not possible with existing neutron spectrometers.

Nanosecond structural dynamics of intrinsically disordered β-casein micelles by neutron spectroscopy

β-casein undergoes a reversible endothermic self-association, forming protein micelles of limited size. In its functional state, a single β-casein monomer is unfolded, which creates a high structural flexibility, which is supposed to play a major role in preventing the precipitation of calcium phosphate particles. We characterize the structural flexibility in terms of nanosecond molecular motions, depending on the temperature by quasielastic neutron scattering. Our major questions are: Does the self-association reduce the chain flexibility? How does the dynamic spectrum of disordered caseins differ from a compactly globular protein? How does the dynamic spectrum of β-casein in solution differ from that of a protein in hydrated powder states? We report on two relaxation processes on a nanosecond and a sub-nanosecond timescale for β-casein in solution. Both processes are analyzed by Brownian oscillator model, by which the spring constant can be defined in the isotropic parabolic potential. The slower process, which is analyzed by neutron spin echo, seems a characteristic feature of the unfolded structure. It requires bulk solvent and is not seen in hydrated protein powders. The faster process, which is analyzed by neutron backscattering, has a smaller amplitude and requires hydration water, which is also observed with folded proteins in the hydrated state. The self-association had no significant influence on internal relaxation, and thus, a β-casein protein monomer flexibility is preserved in the micelle. We derive spring constants of the faster and slower motions of β-caseins in solution and compared them with those of some proteins in various states (folded or hydrated powder).

Fed-batch production of deuterated protein in Escherichia coli for neutron scattering experimentation

Neutron scattering is a powerful technique for determining the structure and dynamics of biological materials in a variety of environmental conditions. A distinguishing property of the neutron is its sensitivity to detecting hydrogen and distinguishing it from its isotope deuterium. This enables unique types of experiments that take advantage of this differential sensitivity called isotopic contrast variation. Using this approach, the chemistry of the system is not changed, but the visibility of individual sample components can be tuned by varying the deuterium content of the system under investigation. Deuterated proteins are commonly produced in bacterial systems that are adapted to growth in D₂O minimal media. To maximize the yield of deuterium-labeled protein and efficiently utilize D₂O and occasionally the deuterated substrate, fed-batch processes are routinely used to maximize biomass production without compromising cell viability. A step-by-step procedure will be described along with a case study of the production of deuterated green fluorescent protein. Limitations of the process will also be discussed.

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.

SANS quantification of bound water in water-soluble polymers across multiple concentration regimes

Contrast-variation small-angle neutron scattering (CV-SANS) is a widely used technique for quantifying hydration water in soft matter systems, but it is predominantly applied in the dilute regime or for systems with a well-defined structure factor. Here, CV-SANS was used to quantify the number of hydration water molecules associating with three water-soluble polymers with different critical solution temperatures and types of water-solute interactions in dilute, semidilute, and concentrated solution through the exploration of novel methods of data fitting and analysis. Multiple SANS fitting workflows with varying levels of model assumptions were evaluated and compared to give insight into SANS model selection. These fitting pathways ranged from general, model-free algorithms to more standard form and structure factor fitting. In addition, Monte Carlo bootstrapping was evaluated as a method to estimate parameter uncertainty through simulation of technical replicates. The most robust fitting workflow for dilute solutions was found to be form factor fitting without CV-SANS (i.e. polymer in 100% D2O). For semidilute and concentrated solutions, while the model-free approach can be mathematically defined for CV-SANS data, the addition of a structure factor imposes physical constraints on the optimization problem, suggesting that the optimal fitting pathway should include appropriate form and structure factor models. The measured hydration numbers were consistent with the number of tightly bound water molecules associated with each monomer unit, and the concentration dependence of the hydration number was largely governed by the chemistry-specific interactions between water and polymer. Polymers with weaker water-polymer interactions (i.e. those with fewer hydration water molecules) were found to have more bound water at higher concentrations than those with stronger water-polymer interactions due to the increase in the number of forced water-polymer contacts in the concentrated system. This SANS-based method to count hydration water molecules can be applied to polymers in any concentration regime, which will lead to improved understanding of water-polymer interactions and their impact on materials design.

Decoupling between the translation and rotation of water in the proximity of a protein molecule

The interaction between water and biomacromolecules is of fundamental interest in biophysics, biochemistry and physical chemistry. By combining neutron scattering and molecular dynamics simulations on a perdeuterated protein at a series of hydration levels, we demonstrated that the translational motion of water is slowed down more significantly than its rotation, when water molecules approach the protein molecule. Further analysis of the simulation trajectories reveals that the observed decoupling results from the fact that the translational motion of water is more correlated over space and more retarded by the charged/polar residues and spatial confinement on the protein surface, than the rotation. Moreover, around the stable protein residues (with smaller atomic fluctuations), water exhibits more decoupled dynamics, indicating a connection between the observed translation-rotation decoupling in hydration water and the local stability of the protein molecule.

Multiscale Microstructure, Composition, and Stability of Surfactant/Polymer Foams

Inclusion of polymer additives is a known strategy to improve foam stability, but questions persist about the amount of polymer incorporated in the foam and the resulting structural changes that impact material performance. Here, we study these questions in sodium dodecyl sulfate (SDS)/hydroxypropyl methylcellulose (HPMC) foams using a combination of flow injection QTOF mass spectrometry and small-angle neutron scattering (SANS) measurements leveraging contrast matching. Mass spectrometry results demonstrate polymer incorporation and retention in the foam during drainage by measuring the HPMC-to-SDS ratio. The results confirm a ratio matching the parent solution and stability over the time of our measurements. The SANS measurements leverage precise contrast matching to reveal detailed descriptions of the micellar structure (size, shape, and aggregation number) along with the foam film thickness. The presence of HPMC leads to thicker films, correlating with increased foam stability over the first 15-20 min after foam production. Taken together, mass spectrometry and SANS present a structural and compositional picture of SDS/HPMC foams and an approach amenable to systematic study for foams, gathering mechanistic insights and providing formulation guidance for rational foam design.

Hydration/Dehydration Behavior of Hydroxyethyl Cellulose Ether in Aqueous Solution

Hydroxyethyl cellulose (HeC) maintains high water solubility over a wide temperature range even in a high temperature region where other nonionic chemically modified cellulose ethers, such as methyl cellulose (MC) and hydroxypropylmethyl cellulose (HpMC), demonstrate cloud points. In order to clarify the reason for the high solubility of HeC, the temperature dependence of the hydration number per glucopyranose unit, n_H, for the HeC samples was examined by using extremely high frequency dielectric spectrum measuring techniques up to 50 GHz over a temperature range from 10 to 70 °C. HeC samples with a molar substitution number (MS) per glucopyranose unit by hydroxyethyl groups ranging from 1.3 to 3.6 were examined in this study. All HeC samples dissolve into water over the examined temperature range and did not show their cloud points. The value of n_H for the HeC sample possessing the MS of 1.3 was 14 at 20 °C and decreased gently with increasing temperature and declined to 10 at 70 °C. The n_H values of the HeC samples are substantially larger than the minimum critical n_H value of ca. 5 necessary to be dissolved into water for cellulose ethers such as MC and HpMC, even in a high temperature range. Then, the HeC molecules possess water solubility over the wide temperature range. The temperature dependence of n_H for the HeC samples and triethyleneglycol, which is a model compound for substitution groups of HeC, is gentle and they are similar to each other. This observation strongly suggests that the hydration/dehydration behavior of the HeC samples was essentially controlled by that of their substitution groups.

Metal-Organic Frameworks for Enzyme Immobilization: Beyond Host Matrix Materials

Enzyme immobilization in metal-organic frameworks (MOFs) as a promising strategy is attracting the interest of scientists from different disciplines with the expansion of MOFs' development. Different from other traditional host materials, their unique strengths of high surface areas, large yet adjustable pore sizes, functionalizable pore walls, and diverse architectures make MOFs an ideal platform to investigate hosted enzymes, which is critical to the industrial and commercial process. In addition to the protective function of MOFs, the extensive roles of MOFs in the enzyme immobilization are being well-explored by making full use of their remarkable properties like well-defined structure, high porosity, and tunable functionality. Such development shifts the focus from the exploration of immobilization strategies toward functionalization. Meanwhile, this would undoubtedly contribute to a better understanding of enzymes in regards to the structural transformation after being hosted in a confinement environment, particularly to the orientation and conformation change as well as the interplay between enzyme and matrix MOFs. In this Outlook, we target a comprehensive review of the role diversities of the host matrix MOF based on the current enzyme immobilization research, along with proposing an outlook toward the future development of this field, including the representatives of potential techniques and methodologies being capable of studying the hosted enzymes.

Anomalous sub-diffusion of water in biosystems: From hydrated protein powders to concentrated protein solution to living cells

Water is essential to life and its translational motion in living systems mediates various biological processes, including transportation of function-required ingredients and facilitating the interaction between biomacromolecules. By combining neutron scattering and isotopic labeling, the present work characterizes translational motion of water on a biomolecular surface, in a range of systems: a hydrated protein powder, a concentrated protein solution, and in living Escherichia coli (E. coli) cells. Anomalous sub-diffusion of water is observed in all samples, which is alleviated upon increasing the water content. Complementary molecular dynamics simulations and coarse-grained numerical modeling demonstrated that the sub-diffusive behavior results from the heterogeneous distribution of microscopic translational mobility of interfacial water. Moreover, by comparing the experimental results measured on E. coli cells with those from a concentrated protein solution with the same amount of water, we show that water in the two samples has a similar average mobility, however the underlying distribution of motion is more heterogeneous in the living cell.

The molecular dynamics of bacterial spore and the role of calcium dipicolinate in core properties at the sub-nanosecond time-scale

Bacterial spores are among the most resistant forms of life on Earth. Their exceptional resistance properties rely on various strategies, among them the core singular structure, organization and hydration. By using elastic incoherent neutron scattering, we probed the dynamics of Bacillus subtilis spores to determine whether core macromolecular motions at the sub-nanosecond timescale could also contribute to their resistance to physical stresses. In addition, in order to better specify the role of the various spore components, we used different mutants lacking essential structure such as the coat (PS4150 mutant), or the calcium dipicolinic acid complex (CaDPA) located in the core (FB122 mutant). PS4150 allows to better probe the core’s dynamics, as proteins of the coat represent an important part of spore proteins, and FB122 gives information about the role of the large CaDPA depot for the mobility of core’s components. We show that core’s macromolecular mobility is not particularly constrained at the sub-nanosecond timescale in spite of its low water content as some dynamical characteristics as force constants are very close to those of vegetative bacteria such as Escherichia coli or to those of fully hydrated proteins. Although the force constants of the coatless mutant are similar to the wild-type’s ones, it has lower mean square displacements (MSDs) at high Q showing that core macromolecules are somewhat more constrained than the rest of spore components. However, no behavior reflecting the glassy state regularly evoked in the literature could be drawn from our data. As hydration and macromolecules’ mobility are highly correlated, the previous assumption, that core low water content might explain spores’ exceptional resistance properties seems unlikely. Thus, we confirm recent theories, suggesting that core water is mostly as free as bulk water and proteins/macromolecules are fully hydrated. The germination of spores leads to a much less stable system with a force constant of 0.1 N/m and MSDs ~2.5 times higher at low Q than in the dormant state. DPA has also an influence on core mobility with a slightly lower force constant for the DPA-less mutant than for the wild-type, and MSDs that are ~ 1.8 times higher on average than for the wild-type at low Q. At high Q, germinated and DPA-less spores were very similar to the wild-type ones, showing that DPA and core compact structure might influence large amplitude motions rather than local dynamics of macromolecules.

Differences between calcium rich and depleted alpha-lactalbumin investigated by molecular dynamics simulations and incoherent neutron scattering

We present a study comparing atomic motional amplitudes in calcium rich and depleted alpha-lactalbumin. The investigations were performed by elastic incoherent neutron scattering (EINS) and molecular dynamics (MD) simulations. As the variations were expected to be very small, three different hydration levels and timescales (instrumental resolutions) were measured. In addition, we used two models to extract the mean square displacements (MSDs) from the EINS data, one taking into account the motional heterogeneity of the MSD. At a timescale of several nanoseconds, small differences in the amplitudes between the calcium enriched and depleted alpha-lactalbumin are visible, whereas at lower timescales no changes can be concluded within the statistics. The results are compared to MD simulations at 280 and 300 K by extracting the MSDs of the trajectories in two separate ways: first by direct calculation, and second by a virtual neutron experiment using the same models as for the experimental data. We show that the simulated data give qualitatively similar results as the experimental data but quantitatively there are differences. Furthermore, the distribution of the MSDs in the simulations suggests that the inclusion of heterogeneity is reasonable for alpha-lactalbumin, but a bi-or trimodal approach may be sufficient.

Unraveling Dynamic Transitions in Time-Resolved Biomolecular Motions by A Dressed Diffusion Model

Recent experimental data reveal the complexity of diffusion dynamics beyond the scope of classical Brownian dynamics. The particles exhibit diverse diffusive motions from the anomalous toward classical diffusion over a wide range of temporal scales. Here a dressed diffusion model is developed to account for non-Brownian phenomena. By coupling the particle dynamics with a local field, the dressed diffusion model generalizes the Langevin equation through coupled damping kernels and generates the salient feature of time-dependent diffusion dynamics reported in the experimental measurements of biomolecules. The dressed diffusion model provides one quantitative aspect for future endeavors in this rapid-growing field.

Protein Backbone and Average Particle Dynamics in Reconstituted Discoidal and Spherical HDL Probed by Hydrogen Deuterium Exchange and Elastic Incoherent Neutron Scattering

Lipoproteins are supramolecular assemblies of proteins and lipids with dynamic characteristics critically linked to their biological functions as plasma lipid transporters and lipid exchangers. Among them, spherical high-density lipoproteins are the most abundant forms of high-density lipoprotein (HDL) in human plasma, active participants in reverse cholesterol transport, and associated with reduced development of atherosclerosis. Here, we employed elastic incoherent neutron scattering (EINS) and hydrogen-deuterium exchange mass spectrometry (HDX-MS) to determine the average particle dynamics and protein backbone local mobility of physiologically competent discoidal and spherical HDL particles reconstituted with human apolipoprotein A-I (apoA-I). Our EINS measurements indicated that discoidal HDL was more dynamic than spherical HDL at ambient temperatures, in agreement with their lipid-protein composition. Combining small-angle neutron scattering (SANS) with contrast variation and MS cross-linking, we showed earlier that the most likely organization of the three apolipoprotein A-I (apoA-I) chains in spherical HDL is a combination of a hairpin monomer and a helical antiparallel dimer. Here, we corroborated those findings with kinetic studies, employing hydrogen-deuterium exchange mass spectrometry (HDX-MS). Many overlapping apoA-I digested peptides exhibited bimodal HDX kinetics behavior, suggesting that apoA-I regions with the same amino acid composition located on different apoA-I chains had different conformations and/or interaction environments.

Molecular origins of bulk viscosity in liquid water

The shear viscous response of water is closely associated with changes in network connectivity on the sub ps timescale. The bulk viscous response is shown here to be associated with local density fluctuations and rotational motion around 1–3 ps.

Secondary structure drives self-assembly in weakly segregated globular protein–rod block copolymers

Imparting secondary structure to the polymer block can drive self-assembly in globular protein–helix block copolymers, increasing the effective segregation strength between blocks with weak or no repulsion.

Mechanism of Trehalose-Induced Protein Stabilization from Neutron Scattering and Modeling

The sugar molecule trehalose has been proven to be an excellent stabilizing cosolute for the preservation of biological materials. However, the stabilizing mechanism of trehalose has been much debated during the previous decades, and it is still not fully understood, partly because it has not been completely established how trehalose molecules structure around proteins. Here, we present a molecular model of a protein-water-trehalose system, based on neutron scattering results obtained from neutron diffraction, quasielastic neutron scattering, and different computer modeling techniques. The structural data clearly show how the proteins are preferentially hydrated, and analysis of the dynamical properties show that the protein residues are slowed down because of reduced dynamics of the protein hydration shell, rather than because of direct trehalose-protein interactions. These findings, thereby, strongly support previous models related to the preferential hydration model and contradict other models based on water replacement at the protein surface. Furthermore, the results are important for understanding the specific role of trehalose in biological stabilization and, more generally, for providing a likely mechanism of how cosolutes affect the dynamics of proteins.

Examination of Trifluoroethanol Interactions with Trp-Cage in Trifluoroethanol-Water at 298 K through Molecular Dynamics Simulations and Intermolecular Nuclear Overhauser Effects

Molecular dynamics simulations of the protein model Trp-cage in 42% trifluoroethanol (TFE)-water at 298 K have been carried out with the goal of exploring peptide hydrogen-solvent fluorine nuclear spin cross-relaxation. The TFE5 model of TFE developed in a previous work was used with the TIP5P-Ew model of water. System densities and component translational diffusion coefficients predicted by the simulations were within 20% of the experimental values. Consideration of the calculated relative amounts of TFE and water surrounding the hydrogens of Trp-cage indicated that the composition of the solvent mixture beyond ∼1.5 nm from the van der Waals surface of the peptide is close to the composition of the bulk solvent, but as observed by others, TFE accumulates preferentially near the peptide surface. In the simulations, both TFE and water molecules make contacts with the peptide surface; water molecules predominate in contacts with the peptide backbone atoms and TFE molecules generally preferentially interact with side chains. Translational diffusion of solvent molecules appears to be slowed near the surface of the peptide. Depending on the location in the structure, TFE molecules form complexes with the peptide that may persist for up to ∼7 ns. Many of the peptide spin-solvent fluorine cross-relaxation parameters (Σ_HF) for which experimental values are available are reasonably well-predicted from the simulations. However, the calculated Σ_HF values were too small for some hydrogens of the 6Trp indole ring and the amino acid hydrogens near this residue in the native structure, whereas Σ_HF values for hydrogens on the side chains of 1Asn, 4Ile, and 7Leu are too large. In 42% TFE-water, persistent conformations of Trp-cage are found, which differ from the conformation found in water by the orientation of the 3Tyr ring.

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.

Molecular motions, structure and hydration behaviour of glucose oligomers in aqueous solution

The presence of helical fragment configurations is confirmed for glucose oligomers in aqueous solution at temperatures lower than 40 °C.

Analysis of elastic incoherent neutron scattering data beyond the Gaussian approximation

This work addresses the use of the Gaussian approximation as a common tool to extract atomic motions in proteins from elastic incoherent neutron scattering and whether improvements in data analysis and additional information can be obtained when going beyond that. We measured alpha-lactalbumin with different levels of hydration on three neutron backscattering spectrometers, to be able to resolve a wide temporal and spatial range for dynamics. We demonstrate that the Gaussian approximation gives qualitatively similar results to models that include heterogeneity, if one respects a certain procedure to treat the intercept of the elastic intensities with the momentum transfer axis. However, the inclusion of motional heterogeneity provides better fits to the data. Our analysis suggests an approach of limited heterogeneity, where including only two kinds of motions appears sufficient to obtain more quantitative results for the mean square displacement. Finally, we note that traditional backscattering spectrometers pose a limit on the lowest accessible momentum transfer. We therefore suggest that complementary information about the spatial evolution of the elastic intensity close to zero momentum transfer can be obtained using other neutron methods, in particular, neutron spin-echo together with polarization analysis.

Neutron scattering in the biological sciences: progress and prospects

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

Overhauser Dynamic Nuclear Polarization for the Study of Hydration Dynamics, Explained

We outline the physical properties of hydration water that are captured by Overhauser Dynamic Nuclear Polarization (ODNP) relaxometry and explore the insights that ODNP yields about the water and the surface that this water is coupled to. As ODNP relies on the pairwise cross-relaxation between the electron spin of a spin probe and a proton nuclear spin of water, it captures the dynamics of single-particle diffusion of an ensemble of water molecules moving near the spin probe. ODNP principally utilizes the same physics as other nuclear magnetic resonance (NMR) relaxometry (i.e., relaxation measurement) techniques. However, in ODNP, electron paramagnetic resonance (EPR) excites the electron spins probes and their high net polarization acts as a signal amplifier. Furthermore, it renders ODNP parameters highly sensitive to water moving at rates commensurate with the EPR frequency of the spin probe (typically 10GHz). Also, ODNP selectively enhances the NMR signal contributions of water moving within close proximity to the spin label. As a result, ODNP can capture ps-ns movements of hydration waters with high sensitivity and locality, even in samples with protein concentrations as dilute as 10 µM. To date, the utility of the ODNP technique has been demonstrated for two major applications: the characterization of the spatial variation in the properties of the hydration layer of proteins or other surfaces displaying topological diversity, and the identification of structural properties emerging from highly disordered proteins and protein domains. The former has been shown to correlate well with the properties of hydration water predicted by MD simulations and has been shown capable of evaluating the hydrophilicity or hydrophobicity of a surface. The latter has been demonstrated for studies of an interhelical loop of proteorhodopsin, the partial structure of α-synuclein embedded at the lipid membrane surface, incipient structures adopted by tau proteins en route to fibrils, and the structure and hydration profile of a transmembrane peptide. This chapter focuses on offering a mechanistic understanding of the ODNP measurement and the molecular dynamics encoded in the ODNP parameters. In particular, it clarifies how the electron-nuclear dipolar coupling encodes information about the molecular dynamics in the nuclear spin self-relaxation and, more importantly, the electron-nuclear spin cross-relaxation rates. The clarification of the molecular dynamics underlying ODNP should assist in establishing a connection to theory and computer simulation that will offer far richer interpretations of ODNP results in future studies.

Nanosecond Tracer Diffusion as a Probe of the Solution Structure and Molecular Mobility of Protein Assemblies: The Case of Ovalbumin

Protein diffusion is not only an important process ensuring biological function but can also be used as a probe to obtain information on structural properties of protein assemblies in liquid solutions. Here, we explore the oligomerization state of ovalbumin at high protein concentrations by means of its short-time self-diffusion. We employ high-resolution incoherent quasielastic neutron scattering to access the self-diffusion on nanosecond timescales, on which interparticle contacts are not altered. Our results indicate that ovalbumin in aqueous (D₂O) solutions occurs in increasingly large assemblies of its monomeric subunits with rising protein concentration. It changes from nearly monomeric toward dimeric and ultimately larger than tetrameric complexes. Simultaneously, we access information on the internal molecular mobility of ovalbumin on the nanometer length scale and compare it with results obtained for bovine serum albumin, immunoglobulin, and β-lactoglobulin.

Gradual Crossover from Subdiffusion to Normal Diffusion: A Many-Body Effect in Protein Surface Water

Dynamics of hydration water is essential for the function of biomacromolecules. Previous studies have demonstrated that water molecules exhibit subdiffusion on the surface of biomacromolecules; yet the microscopic mechanism remains vague. Here, by performing neutron scattering, molecular dynamics simulations, and analytic modeling on hydrated perdeuterated protein powders, we found water molecules jump randomly between trapping sites on protein surfaces, whose waiting times obey a broad distribution, resulting in subdiffusion. Moreover, the subdiffusive exponent gradually increases with observation time towards normal diffusion due to a many-body volume-exclusion effect.