Dynamics of protein and its hydration water: neutron scattering studies on fully deuterated GFP

Backscattering silicon spectrometer (BASIS): sixteen years in advanced materials characterization

Quasielastic neutron scattering (QENS) is an experimental technique that can measure parameters of mobility, such as diffusion jump rate and jump length, as well as localized relaxations of chemical species (molecules, ions, and segments) at atomic and nanometer length scales. Due to the high penetrative power of neutrons and their sensitivity to neutron scattering cross-section of chemical species, QENS can effectively probe mobility inside most bulk materials. This review focuses on QENS experiments performed using a neutron backscattering silicon spectrometer (BASIS) to explore the dynamics in various materials and understand their structure-property relationship. BASIS is a time-of-flight near-backscattering inverted geometry spectrometer with very high energy resolution (approximately 0.0035 meV of full width at half maximum), allowing measurements of dynamics on nano to picosecond timescales. The science areas studied with BASIS are diverse, with a focus on soft matter topics, including traditional biological and polymer science experiments, as well as measurements of fluids ranging from simple hydrocarbons and aqueous solutions to relatively complex room-temperature ionic liquids and deep-eutectic solvents, either in the bulk state or confined. Additionally, hydrogen confined in various materials is routinely measured on BASIS. Other topics successfully investigated at BASIS include quantum fluids, spin glasses, and magnetism. BASIS has been in the user program since 2007 at the Spallation Neutron Source of the Oak Ridge National Laboratory, an Office of Science User Facility supported by the U.S. Department of Energy. Over the past sixteen years, BASIS has contributed to various scientific disciplines, exploring the structure and dynamics of many chemical species and their fabrication for practical applications. A comprehensive review of BASIS contributions and capabilities would be an asset to the materials science community, providing insights into employing the neutron backscattering technique for advanced materials characterization.

Decoupling of the onset of anharmonicity between a protein and its surface water around 200 K

The protein dynamical transition at ~200 K, where the biomolecule transforms from a harmonic, non-functional form to an anharmonic, functional state, has been thought to be slaved to the thermal activation of dynamics in its surface hydration water. Here, by selectively probing the dynamics of protein and hydration water using elastic neutron scattering and isotopic labeling, we found that the onset of anharmonicity in the two components around 200 K is decoupled. The one in protein is an intrinsic transition, whose characteristic temperature is independent of the instrumental resolution time, but varies with the biomolecular structure and the amount of hydration, while the one of water is merely a resolution effect.

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.

Direct Observation of the Mutual Coupling Effect in the Protein-Water-Glycerol Mixture by Combining Neutron Scattering and Selective Deuteration

Numerous studies have discussed the impact of cosolvents on the structure, dynamics, and stability of proteins in aqueous solutions. However, the dynamics of cosolvents in the protein-water-cosolvent ternary system is largely unexplored in experiments due to technical difficulty. Consequently, a comprehensive understanding of the interplay among proteins, water, and cosolvents is still lacking. Here, we employed selective deuteration and neutron scattering techniques to characterize the individual motions of each component in the protein/water/glycerol (GLY) mixture across various temperatures. The consistent dynamic onset temperatures and the correlation between the MSD of the protein and the viscosity of solvents revealed the mutual coupling effects among the three components. Furthermore, our experimental and simulation results showed that the hydrogen bond relaxation energy barrier in the ternary system is ∼43 kJ/mol, whereas in the protein-water binary system it is merely ∼35 kJ/mol. Therefore, we suggest that GLY can enhance hydrogen bond interactions in the ternary system through the mutual coupling effect, thereby serving as one of the protective mechanisms of protein preservation by GLY.

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.

Examination of Solvent Interactions with Trp-Cage in 1,1,1,3,3,3-Hexafluoro-2-propanol-water at 298 K through MD Simulations and Intermolecular Nuclear Overhauser Effects

MD simulations of the peptide Trp-cage dissolved in 28% hexafluoroisopropanol (HFIP)-water have been carried out at 298 K with the goal of exploring peptide hydrogen-solvent fluorine nuclear spin cross-relaxation. The work was motivated by the observation that most experimental fluoroalcohol-peptide cross-relaxation terms at 298 K are small, both positive and negative, and not always well predicted from simulations. The cross-relaxation terms for hydrogens of the caged tryptophan residue of Trp-cage are substantially negative, a result consistent with simulations. It was concluded that hexafluoroisopropanol interactions near this part of the peptide are particularly long-lived. While both HFIP and water are present in all regions of the simulation box, the composition of the solvent mixture is not homogeneous throughout the system. HFIP generally accumulates near the peptide surface, while water molecules are preferentially found in regions that are more than 1.5 nm from the surface of the peptide. However, some water remains in higher-than-expected amounts in the solvent layer surrounding 6Trp, 9Asp, Ser13, and Ser14 residues in the helical region of Trp-cage. As observed in simulations of this system at 278 K, HFIP molecules aggregate into clusters that continually form and re-form. Translational diffusion of both HFIP and water appears to be slowed near the surface of the peptide with reduction in diffusion near the 6Trp residue 2- to 3-fold larger than calculated for solvent interactions with other regions of Trp-cage.

Creatinase: Using Increased Entropy to Improve the Activity and Thermostability

Improving protein thermostability in mutagenesis-based enzyme engineering was often achieved by enhancing interresidue interactions via mutation to increase the enthalpy penalty of unfolding. However, this approach may trade off the functional activity due to the loss of structural flexibility of the biomolecule. Here, by performing X-ray crystallography, enzymatic kinetic experiments, neutron scattering, and thermodynamical measurements, we compared the structures, catalytic behaviors, dynamics, and thermostability between a wild-type creatinase and its four-point mutant. We found that the mutant is an entropy-driven thermostable protein with higher structural flexibility, i.e., higher conformational entropy, in the folded state compared to the wild type. The increased conformational entropy of the mutant in the folded state can reduce the entropy gain during unfolding and thus renders it greater thermostability. Moreover, the increased structural flexibility, particularly around the catalytic site, can broaden the mutant's working temperature range and considerably improve its activity at ambient conditions, which is crucial for its application in diagnosing kidney diseases. Complementary all-atom molecular dynamics simulations indicated that the four mutations replaced several of the strong interresidue interactions (electrostatic interactions and hydrogen bonds) with weak hydrophobic interactions. These substitutions not only release the structural flexibility to promote the thermostability and enzymatic activity of the protein but they also preserve the protein structure from collapsing. Our findings may pave a route for the entropy-driven strategy to design proteins with high thermostability and activity.

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.

Vibrational Analysis of Hydration-Layer Water around Ubiquitin, Unpeeled Layer by Layer: Molecular-Dynamics Perceptions

Classical molecular-dynamics simulations have been performed to examine the interplay between ubiquitin and its hydration-water sub-layers, chiefly from a vibrational-mode and IR viewpoint-where we analyse individual sub-layers characteristics. The vibrational Density of States (VDOS) revealed that the first solvation sub-shell indicates a confined character therein. For layers of increasing distance from the surface, the adoption of greater bulk-like spectral behaviour was evident, suggesting that vibrational harmonisation to bulk occurs within 6-7 Å of the surface.

Tools shaping drug discovery and development

Spectroscopic, scattering, and imaging methods play an important role in advancing the study of pharmaceutical and biopharmaceutical therapies. The tools more familiar to scientists within industry and beyond, such as nuclear magnetic resonance and fluorescence spectroscopy, serve two functions: as simple high-throughput techniques for identification and purity analysis, and as potential tools for measuring dynamics and structures of complex biological systems, from proteins and nucleic acids to membranes and nanoparticle delivery systems. With the expansion of commercial small-angle x-ray scattering instruments into the laboratory setting and the accessibility of industrial researchers to small-angle neutron scattering facilities, scattering methods are now used more frequently in the industrial research setting, and probe-less time-resolved small-angle scattering experiments are now able to be conducted to truly probe the mechanism of reactions and the location of individual components in complex model or biological systems. The availability of atomic force microscopes in the past several decades enables measurements that are, in some ways, complementary to the spectroscopic techniques, and wholly orthogonal in others, such as those related to nanomechanics. As therapies have advanced from small molecules to protein biologics and now messenger RNA vaccines, the depth of biophysical knowledge must continue to serve in drug discovery and development to ensure quality of the drug, and the characterization toolbox must be opened up to adapt traditional spectroscopic methods and adopt new techniques for unraveling the complexities of the new modalities. The overview of the biophysical methods in this review is meant to showcase the uses of multiple techniques for different modalities and present recent applications for tackling particularly challenging situations in drug development that can be solved with the aid of fluorescence spectroscopy, nuclear magnetic resonance spectroscopy, atomic force microscopy, and small-angle scattering.

Structure and diffusive dynamics of aspartate α-decarboxylase (ADC) liganded with D-serine in aqueous solution

Incoherent neutron spectroscopy, in combination with dynamic light scattering, was used to investigate the effect of ligand binding on the center-of-mass self-diffusion and internal diffusive dynamics of Escherichia coli aspartate α-decarboxylase (ADC). The X-ray crystal structure of ADC in complex with the d-serine inhibitor was also determined, and molecular dynamics simulations were used to further probe the structural rearrangements that occur as a result of ligand binding. These experiments reveal that d-serine forms hydrogen bonds with some of the active site residues, that higher order oligomers of the ADC tetramer exist on ns–ms time-scales, and also show that ligand binding both affects the ADC internal diffusive dynamics and appears to further increase the size of the higher order oligomers.

Neutron spectroscopy, dynamic light scattering, X-ray diffraction, and MD-simulations were used to investigate the effect of ligand binding on the structure and diffusive dynamics of Escherichia coli aspartate alpha-decarboxylase.

The Impact of Surface Charges of Carboxylated Cellulose Nanofibrils on the Water Motions in Hydrated Films

Cellulose nanofibrils (CNFs) with carboxylated surface ligands are a class of materials with tunable surface functionality, good mechanical properties, and bio-/environmental friendliness. They have been used in many applications as scaffold, reinforcing, or functional materials, where the interaction between adsorbed moisture and the CNF could lead to different properties and structures and become critical to the performance of the materials. In this work, we exploited multiple experimental methods to study the water movement in hydrated films made of carboxylated CNFs prepared by TEMPO oxidation with two different surface charges of 600 and 1550 μmol·g^-1. A combination of quartz crystal microbalance with dissipation (QCM-D) and small-angle X-ray scattering (SAXS) shows that both the surface charge of a single fibril and the films' network structure contribute to the moisture uptake. The films with 1550 μmol·g^-1 surface charges take up twice the amount of moisture per unit mass, leading to the formation of nanostructures with an average radius of gyration of 2.1 nm. Via the nondestructive quasi-elastic neutron scattering (QENS), a faster motion is explained as a localized movement of water molecules inside confined spheres, and a slow diffusive motion is found with the diffusion coefficient close to bulk water at room temperature via a random jump diffusion model and regardless of the surface charge in films made from CNFs.

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.

Low-temperature librations and dynamical transition in proteins at differing hydration levels

Hydration of water affects the dynamics and in turn the activity of biomacromolecules. We investigated the dependence of the librational oscillations and the dynamical transition on the hydrating conditions of two globular proteins with different structure and size, namely β-lactoglobulin (βLG) and human serum albumin (HSA), by spin-label electron paramagnetic resonance (EPR) in the temperature range of 120–270 K. The proteins were spin-labeled with 5-maleimide spin-label on free cysteins and prepared in the lyophilized state, at low (h = 0.12) and full (h = 2) hydration levels in buffer. The angular amplitudes of librations are small and almost temperature independent for both lyophilized proteins. Therefore, in these samples, the librational dynamics is restricted and the dynamical transition is absent. In the small and compact beta-structured βLG, the angular librational amplitudes increase with temperature and hydrating conditions, whereas hydration-independent librational oscillations whose amplitudes rise with temperature are recorded in the large and flexible alpha-structured HSA. Both βLG and HSA at low and fully hydration levels undergo the dynamical transition at about 230 K. The overall results indicate that protein librational dynamics is activated at the low hydration level h = 0.12 and highlight biophysical properties that are common to other biosamples at cryogenic temperatures.

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.

Probing Small-Angle Molecular Motions with EPR Spectroscopy: Dynamical Transition and Molecular Packing in Disordered Solids

Disordered molecular solids present a rather broad class of substances of different origin—amorphous polymers, materials for photonics and optoelectronics, amorphous pharmaceutics, simple molecular glass formers, and others. Frozen biological media in many respects also may be referred to this class. Theoretical description of dynamics and structure of disordered solids still does not exist, and only some phenomenological models can be developed to explain results of particular experiments. Among different experimental approaches, electron paramagnetic resonance (EPR) applied to spin probes and labels also can deliver useful information. EPR allows probing small-angle orientational molecular motions (molecular librations), which intrinsically are inherent to all molecular solids. EPR is employed in its conventional continuous wave (CW) and pulsed—electron spin echo (ESE)—versions. CW EPR spectra are sensitive to dynamical librations of molecules while ESE probes stochastic molecular librations. In this review, different manifestations of small-angle motions in EPR of spin probes and labels are discussed. It is shown that CW-EPR-detected dynamical librations provide information on dynamical transition in these media, similar to that explored with neutron scattering, and ESE-detected stochastic librations allow elucidating some features of nanoscale molecular packing. The possible EPR applications are analyzed for gel-phase lipid bilayers, for biological membranes interacting with proteins, peptides and cryoprotectants, for supercooled ionic liquids (ILs) and supercooled deep eutectic solvents (DESs), for globular proteins and intrinsically disordered proteins (IDPs), and for some other molecular solids.

Seeking Solvation: Exploring the Role of Protein Hydration in Silk Gelation

The mechanism by which arthropods (e.g., spiders and many insects) can produce silk fibres from an aqueous protein (fibroin) solution has remained elusive, despite much scientific investigation. In this work, we used several techniques to explore the role of a hydration shell bound to the fibroin in native silk feedstock (NSF) from Bombyx mori silkworms. Small angle X-ray and dynamic light scattering (SAXS and DLS) revealed a coil size (radius of gyration or hydrodynamic radius) around 12 nm, providing considerable scope for hydration. Aggregation in dilute aqueous solution was observed above 65 °C, matching the gelation temperature of more concentrated solutions and suggesting that the strength of interaction with the solvent (i.e., water) was the dominant factor. Infrared (IR) spectroscopy indicated decreasing hydration as the temperature was raised, with similar changes in hydration following gelation by freezing or heating. It was found that the solubility of fibroin in water or aqueous salt solutions could be described well by a relatively simple thermodynamic model for the stability of the protein hydration shell, which suggests that the affected water is enthalpically favoured but entropically penalised, due to its reduced (vibrational or translational) dynamics. Moreover, while the majority of this investigation used fibroin from B. mori, comparisons with published work on silk proteins from other silkworms and spiders, globular proteins and peptide model systems suggest that our findings may be of much wider significance.

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⁻¹.

Notes on Fitting and Analysis Frameworks for QENS Spectra of (Soft) Colloid Suspensions

With continuously improving signal-to-noise ratios, a statistically sound analysis of quasi-elastic neutron scattering (QENS) spectra requires to fit increasingly complex models which poses several challenges. Simultaneous fits of the spectra for all recorded values of the momentum transfer become a standard approach. Spectrometers at spallation sources can have a complicated non-Gaussian resolution function which has to be described most accurately. At the same time, to speed up the fitting, an analytical convolution with this resolution function is of interest. Here, we discuss basic concepts to efficient approaches for fits of QENS spectra based on standard MATLAB and Python fit algorithms. We illustrate the fits with example data from IN16B, BASIS, and BATS.

Experimental mapping of short-wavelength phonons in proteins

Phonons are quasi-particles, observed as lattice vibrations in periodic materials, that often dampen in the presence of structural perturbations. Nevertheless, phonon-like collective excitations exist in highly complex systems, such as proteins, although the origin of such collective motions has remained elusive. Here we present a picture of temperature and hydration dependence of collective excitations in green fluorescent protein (GFP) obtained by inelastic neutron scattering. Our results provide evidence that such excitations can be used as a measure of flexibility/softness and are possibly associated with the protein’s activity. Moreover, we show that the hydration water in GFP interferes with the phonon propagation pathway, enhancing the structural rigidity and stability of GFP.

•

Quantum phenomena in biology have long fascinated people around the world

•

This work presents a direct experimental observation of phonons, the quantum vibrations in a protein

•

The collective excitations or phonons in proteins were detected by utilizing inelastic neutron scattering technique at Oak Ridge National Laboratory

•

Our results illustrate the flexibility-activity relationship in proteins by mapping the temperature and hydration dependence of these collective excitations

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.

Hydration of methemoglobin studied by in silico modeling and dielectric spectroscopy

The hemoglobin concentration of 35 g/dl of human red blood cells is close to the solubility threshold. Using microwave dielectric spectroscopy, we have assessed the amount of water associated with hydration shells of methemoglobin as a function of its concentration in the presence or absence of ions. We estimated water-hemoglobin interactions to interpret the obtained data. Within the concentration range of 5-10 g/dl of methemoglobin, ions play an important role in defining the free-to-bound water ratio competing with hemoglobin to recruit water molecules for the hydration shell. At higher concentrations, hemoglobin is a major contributor to the recruitment of water to its hydration shell. Furthermore, the amount of bound water does not change as the hemoglobin concentration is increased from 15 to 30 g/dl, remaining at the level of ∼20% of the total intracellular water pool. The theoretical evaluation of the ratio of free and bound water for the hemoglobin concentration in the absence of ions corresponds with the experimental results and shows that the methemoglobin molecule binds about 1400 water molecules. These observations suggest that within the concentration range close to the physiological one, hemoglobin molecules are so close to each other that their hydration shells interact. In this case, the orientation of the hemoglobin molecules is most likely not stochastic, but rather supports partial neutralization of positive and negative charges at the protein surface. Furthermore, deformation of the red blood cell shape results in the rearrangement of these structures.

Examination of Interactions of Hexafluoro-2-propanol with Trp-Cage in Hexafluoro-2-propanol-Water by MD Simulations and Intermolecular Nuclear Overhauser Effects

All-atom molecular dynamic simulations of the peptide Trp-cage in 30% hexafluoro-2-propanol- water (V/V) at 278 K have been carried out with the goal of exploring peptide hydrogen-solvent fluorine nuclear spin cross relaxation. Force field parameters for HFIP reported by Fioroni et al. along with the fluorine parameters of the TFE5 model reported by this lab were used. Water was represented by the TIP5P-Ew model. Peptide modeling used the AMBER99SB-ILDN force field. Translational diffusion coefficients of solution components at 278 K were predicted to within 35% of experimental values using these parameter sets. The simulations indicate that the solvent mixture is not homogeneous, with HFIP molecules clustered into aggregates as large as 53 fluoroalcohol molecules. The solvent environment of surface atoms of Trp-cage fluctuates between being HFIP-rich and more water-rich about every 10 ns. In accord with previous studies by other groups, the average concentration of HFIP near the surface of the peptide is significantly enhanced over the concentration of HFIP in the bulk solvent. In the simulations, ∼7% of the initial contacts between HFIP molecules and Trp-cage develop into peptide-fluoroalcohol interactions that persist for times as long as 8 ns. Most of the available experimental nuclear spin cross-relaxation rates (Σ_HF) for hydrogens of the Trp-cage in 30% HFIP-water are reproduced from the MD trajectories to within uncertainties of the experimental data and the simulations. However, a few calculated Σ_HF values for hydrogens of the Trp-cage do not agree with experiment. These tend to be situations where long-lived peptide-HFIP interactions are predicted. The disagreements between observed and calculated Σ_HF in these instances signal defects in the modeling parameters and procedures that are presently unrecognized.

Quasielastic neutron scattering studies on couplings of protein and water dynamics in hydrated elastin

Performing quasielastic neutron scattering measurements and analyzing both elastic and quasielasic contributions, we study protein and water dynamics of hydrated elastin. At low temperatures, hydration-independent methyl group rotation dominates the findings. It is characterized by a Gaussian distribution of activation energies centered at about E_m = 0.17 eV. At ∼195 K, coupled protein-water motion sets in. The hydration water shows diffusive motion, which is described by a Gaussian distribution of activation energies with E_m = 0.57 eV. This Arrhenius behavior of water diffusion is consistent with previous results for water reorientation, but at variance with a fragile-to-strong crossover at ∼225 K. The hydration-related elastin backbone motion is localized and can be attributed to the cage rattling motion. We speculate that its onset at ∼195 K is related to a secondary glass transition, which occurs when a β relaxation of the protein has a correlation time of τ_β ∼ 100 s. Moreover, we show that its temperature-dependent amplitude has a crossover at the regular glass transition T_g = 320 K of hydrated elastin, where the α relaxation of the protein obeys τ_α ∼ 100 s. By contrast, we do not observe a protein dynamical transition when water dynamics enters the experimental time window at ∼240 K.

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.

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.

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.

X-ray crystallographic studies on the hydrogen isotope effects of green fluorescent protein at sub-ångström resolutions

Hydrogen atoms are critical to the nature and properties of proteins, and thus deuteration has the potential to influence protein function. In fact, it has been reported that some deuterated proteins show different physical and chemical properties to their protiated counterparts. Consequently, it is important to investigate protonation states around the active site when using deuterated proteins. Here, hydrogen isotope effects on the S65T/F99S/M153T/V163A variant of green fluorescent protein (GFP), in which the deprotonated B form is dominant at pH 8.5, were investigated. The pH/pD dependence of the absorption and fluorescence spectra indicates that the protonation state of the chromophore is the same in protiated GFP in H₂O and protiated GFP in D₂O at pH/pD 8.5, while the pK_a of the chromophore became higher in D₂O. Indeed, X-ray crystallographic analyses at sub-ångström resolution revealed no apparent changes in the protonation state of the chromophore between the two samples. However, detailed comparisons of the hydrogen OMIT maps revealed that the protonation state of His148 in the vicinity of the chromophore differed between the two samples. This indicates that protonation states around the active site should be carefully adjusted to be the same as those of the protiated protein when neutron crystallographic analyses of proteins are performed.

Role of hydration water in the onset of protein structural dynamics

Proteins are the molecular workhorses in a living organism. Their 3D structures are animated by a multitude of equilibrium fluctuations and specific out-of-equilibrium motions that are required for proteins to be biologically active. When studied as a function of temperature, functionally relevant dynamics are observed at and above the so-called protein dynamical transition (~240 K) in hydrated, but not in dry proteins. In this review we present and discuss the main experimental and computational results that provided evidence for the dynamical transition, with a focus on the role of hydration water dynamics in sustaining functional protein dynamics. The coupling and mutual influence of hydration water dynamics and protein dynamics are discussed and the hypotheses illustrated that have been put forward to explain the physical origin of their onsets.

Water in protein hydration and ligand recognition

This review describes selected basics of water in biomolecular recognition. We focus on a qualitative understanding of the most important physical aspects, how these change in magnitude between bulk water and protein environment, and how the roles that water plays for proteins arise from them. These roles include mechanical support, thermal coupling, dielectric screening, mass and charge transport, and the competition with a ligand for the occupation of a binding site. The presence or absence of water has ramifications that range from the thermodynamic binding signature of a single ligand up to cellular survival. The large inhomogeneity in water density, polarity and mobility around a solute is hard to assess in experiment. This is a source of many difficulties in the solvation of protein models and computational studies that attempt to elucidate or predict ligand recognition. The influence of water in a protein binding site on the experimental enthalpic and entropic signature of ligand binding is still a point of much debate. The strong water‐water interaction in enthalpic terms is counteracted by a water molecule's high mobility in entropic terms. The complete arrest of a water molecule's mobility sets a limit on the entropic contribution of a water displacement process, while the solvent environment sets limits on ligand reactivity.

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.

Dynamics from elastic neutron-scattering via direct measurement of the running time-integral of the van Hove distribution function

We present a new neutron-scattering approach to access the van Hove distribution function directly in the time domain, I(t), which reflects the system dynamics. Currently, I(t) is essentially determined from neutron energy-exchange. Our method consists of the straightforward measurement of the running time-integral of I(t), by computing the portion of scattered neutrons corresponding to species at rest within a time t, (conceptually elastic scattering). Previous attempts failed to recognise this connection. Starting from a theoretical standpoint, a practical realisation is assessed via numerical methods and an instrument simulation.

Water dynamics in MCF-7 breast cancer cells: a neutron scattering descriptive study

Water mobility in cancer cells could be a powerful parameter to predict the progression or remission of tumors. In the present descriptive work, new insight into this concept was achieved by combining neutron scattering and thermal analyses. The results provide the first step to untangle the role played by water dynamics in breast cancer cells (MCF-7) after treatment with a chemotherapy drug. By thermal analyses, the cells were probed as micrometric reservoirs of bulk-like and confined water populations. Under this perspective we showed that the drug clearly alters the properties of the confined water. We have independently validated this idea by accessing the cellular water dynamics using inelastic neutron scattering. Finally, analysis of the quasi-elastic neutron scattering data allows us to hypothesize that, in this particular cell line, diffusion increases in the intracellular water in response to the action of the drug on the nanosecond timescale.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

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.

Dynamics Properties of Photosynthetic Microorganisms Probed by Incoherent Neutron Scattering

Studies on the dynamical properties of photosynthetic membranes of land plants and purple bacteria have been previously performed by neutron spectroscopy, revealing a tight coupling between specific photochemical reactions and macromolecular dynamics. Here, we probed the intrinsic dynamics of biotechnologically useful mutants of the green alga Chlamydomonas reinhardtii by incoherent neutron scattering coupled with prompt chlorophyll fluorescence experiments. We brought to light that single amino acid replacements in the plastoquinone (PQ)-binding niche of the photosystem II D1 protein impair electron transport (ET) efficiency between quinones and confer increased flexibility to the host membranes, expanding to the entire cells. Hence, a more flexible environment in the PQ-binding niche has been associated to a less efficient ET. A similar function/dynamics relationship was also demonstrated in Rhodobacter sphaeroides reaction centers having inhibited ET, indicating that flexibility at the quinones region plays a crucial role in evolutionarily distant organisms. Instead, a different functional/dynamical correlation was observed in algal mutants hosting a single amino acid replacement residing in a D1 domain far from the PQ-binding niche. Noteworthy, this mutant displayed the highest degree of flexibility, and besides having a nativelike ET efficiency in physiological conditions, it acquired novel, to our knowledge, phenotypic traits enabling it to preserve a high maximal quantum yield of photosystem II photochemistry in extreme habitats. Overall, in the nanosecond timescale, the degree of the observed flexibility is related to the mutation site; in the picosecond timescale, we highlighted the presence of a more pronounced dynamic heterogeneity in all mutants compared to the native cells, which could be related to a marked chemically heterogeneous environment.

Dynamics of a family of cyan fluorescent proteins probed by incoherent neutron scattering

Cyan fluorescent proteins (CFPs) are variants of green fluorescent proteins in which the central tyrosine of the chromophore has been replaced by a tryptophan. The increased bulk of the chromophore within a compact protein and the change in the positioning of atoms capable of hydrogen bonding have made it difficult to optimize their fluorescence properties, which took approximately 15 years between the availability of the first useable CFP, enhanced cyan fluorescent protein (ECFP), and that of a variant with almost perfect fluorescence efficiency, mTurquoise2. To understand the molecular bases of the progressive improvement in between these two CFPs, we have studied by incoherent neutron scattering the dynamics of five different variants exhibiting progressively increased fluorescence efficiency along the evolution pathway. Our results correlate well with the analysis of the previously determined X-ray crystallographic structures, which show an increase in flexibility between ECFP and the second variant, Cerulean, which is then hindered in the three later variants, SCFP3A (Super Cyan Fluorescent Protein 3A), mTurquoise and mTurquoise2. This confirms that increasing the rigidity of the direct environment of the fluorescent chromophore is not the sole parameter leading to brighter fluorescent proteins and that increased flexibility in some cases may be helpful.

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.

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.

Dynamical transition in molecular glasses and proteins observed by spin relaxation of nitroxide spin probes and labels

In glassy substances and biological media, dynamical transitions are observed in neutron scattering that manifests itself as deviations of the translational mean-squared displacement, 〈x²〉, of hydrogen atoms from harmonic dynamics. In biological media, the deviation occurs at two temperature intervals, at ∼100-150 K and at ∼170-230 K, and it is attributed to the motion of methyl groups in the former case and to the transition from harmonic to anharmonic or diffusive motions in the latter case. In this work, electron spin echo (ESE) spectroscopy-a pulsed version of electron paramagnetic resonance-is applied to study the spin relaxation of nitroxide spin probes and labels introduced in molecular glass former o-terphenyl and in protein lysozyme. The anisotropic contribution to the rate of the two-pulse ESE decay, ΔW, is induced by spin relaxation appearing because of restricted orientational stochastic molecular motion; it is proportional to 〈α²〉τ_c, where 〈α²〉 is the mean-squared angle of reorientation of the nitroxide molecule around the equilibrium position and τ_c is the correlation time of reorientation. The ESE time window allows us to study motions with τ_c < 10^-7 s. For glassy o-terphenyl, the 〈α²〉τ_c temperature dependence shows a transition near 240 K, which is in agreement with the literature data on 〈x²〉. For spin probes of essentially different size, the obtained data were found to be close, which evidences that motion is cooperative, involving a nanocluster of several neighboring molecules. For the dry lysozyme, the 〈α²〉τ_c values below 260 K were found to linearly depend on the temperature in the same way as it was observed in neutron scattering for 〈x²〉. As spin relaxation is influenced only by stochastic motion, the harmonic motions seen in ESE must be overdamped. In the hydrated lysozyme, ESE data show transitions near 130 K for all nitroxides, near 160 K for the probe located in the hydration layer, and near 180 K for the label in the protein interior. For this system, the two latter transitions are not observed in neutron scattering. The ESE-detected transitions are suggested to be related with water dynamics in the nearest hydration shell: with water glass transition near 130 K and with the onset of overall water molecular reorientations near 180 K; the disagreement with neutron scattering is ascribed to the larger time window for ESE-detected motions.

Protein crystal lattices are dynamic assemblies: the role of conformational entropy in the protein condensed phase

Until recently, the occurrence of conformational entropy in protein crystal contacts was considered to be a very unlikely event. A study based on the most accurately refined protein structures demonstrated that side-chain conformational entropy and static disorder might be common in protein crystal lattices. The present investigation uses structures refined using ensemble refinement to show that although paradoxical, conformational entropy is likely to be the major factor in the emergence and integrity of the protein condensed phase. This study reveals that the role of shape entropy and local entropic forces expands beyond the onset of crystallization. For the first time, the complete pattern of intermolecular interactions by protein atoms in crystal lattices is presented, which shows that van der Waals interactions dominate in crystal formation.