Hydrogen-bond dynamics at the bio-water interface in hydrated proteins: a molecular-dynamics study

Self-assembly of tamarind seed polysaccharide via enzymatic depolymerization and degalactosylation enhanced ice recrystallization inhibition activity

Polysaccharides are becoming potential candidates for developing food-grade cryoprotectants due to their extensive accessibility and health-promoting effects. However, unremarkable ice recrystallization inhibition (IRI) activity and high viscosity limit their practical applications in some systems. Our previous study found a galactoxyloglucan polysaccharide from tamarind seed (TSP) showing moderate IRI activity. Herein, the enhancement of the IRI performance of TSP via enzymatic depolymerization and degalactosylation-induced self-assembly was reported. TSP was depolymerized and subsequently removed ∼40 % Gal, which induced the formation of supramolecular rod-like fiber self-assembles and exhibited a severalfold enhancement of IRI. Ice shaping assay did not show obvious faceting of ice crystals, indicating that both depolymerized and self-assembled TSP showed very weak binding to ice. Molecular dynamics simulation confirmed the absence of molecular complementarity with ice. Further, it highlighted that degalactosylation did not cause significant changes in local hydration properties of TSP from the view of a single oligomer. The inconsistency between molecular simulation and macroscopic IRI effect proposed that the formation of unique supramolecular self-assemblies may be a key requirement for enhancing IRI activity. The findings of this study provided a new opportunity to enhance the applied potential of natural polysaccharides in food cryoprotection.

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.

Temperature-Dependent Dynamics at Protein-Solvent Interfaces

We perform differential scanning calorimetry, broadband dielectric spectroscopy (BDS), and nuclear magnetic resonance (NMR) studies to ascertain the molecular dynamics in mixtures of ethylene glycol with elastin or lysozyme over broad temperature ranges. To focus on the protein-solvent interface, we use mixtures with about equal numbers of amino acids and solvent molecules. The elastin and lysozyme mixtures show similar glass transition steps, which extend over a broad temperature range of 157-185K. The BDS and NMR studies yield fully consistent results for the fastest process P1, which is caused by the structural relaxation of ethylene glycol between the protein molecules and follows an Arrhenius law with an activation energy of Ea=0.63eV. It involves quasi-isotropic reorientation and is very similar in the elastin and lysozyme matrices but different from the alpha and beta relaxations of bulk ethylene glycol. Two slower BDS processes P2 and P3 have protein-dependent time scales, but exhibit a similar Arrhenius-like temperature dependence with an activation energy of Ea~0.81eV. However, P2 and P3 do not have a clear NMR signature. In particular, the NMR results for the lysozyme mixture reveal that the protein backbone does not show isotropic alpha-like motion on the P2 and P3 time scales but only restricted beta-like reorientation. The different activation energies of the P1 and P2/P3 processes do not support an intimate coupling of protein and ethylene glycol dynamics. The present results are compared with previous findings for mixtures of proteins with water or glycerol, implying qualitatively different dynamical couplings at various protein-solvent interfaces.

Ultrafast Solar-Vapor Harvesting Based on a Hierarchical Porous Hydrogel with Wettability Contrast and Tailored Water States

Optimizing the water bonding network in an evaporator is significant for efficient solar-driven vapor generation (SVG). Herein, we report a facile one-pot method to regulate the hydrated structure and wettability in a hierarchical porous hydrogel. An ovalbumin (OVA)-polyacrylamide hydrogel foam was fabricated in a cake-making fashion. Because of the enrichment of amphiphilic OVA at the interface, the hydrophobic walls of the air pores in the foam provide vaporization sites and help reduce parasitic heat loss, while the hydrophilic skeleton with the secondary pores effectively pumps capillary water. Notably, the proportion of intermediate water in the foam reaches 87.6% with the melting point as low as -10 °C. All these features contribute to an exceptional evaporation rate of 3.4-4.5 kg m^-2 h^-1 under 1 sun and robust SVG performances at high-humidity, weak sunlight, or cold weathers. The strategy of using amphiphilic molecules to optimize the hydrated structures both at the interface and in bulk promises the reasonable design of SVG materials with superior efficiency and weather adaptability.

Structural and Dynamical Properties of Liquids in Confinements: A Review of Molecular Dynamics Simulation Studies

Molecular dynamics (MD) simulations are a powerful tool for detailed studies of altered properties of liquids in confinement, in particular, of changed structures and dynamics. They allow, on one hand, for perfect control and systematic variation of the geometries and interactions inherent in confinement situations and, on the other hand, for type-selective and position-resolved analyses of a huge variety of structural and dynamical parameters. Here, we review MD simulation studies on various types of liquids and confinements. The main focus is confined aqueous systems, but also ionic liquids and polymer and silica melts are discussed. Results for confinements featuring different interactions, sizes, shapes, and rigidity will be presented. Special attention will be given to situations in which the confined liquid and the confining matrix consist of the same type of particles and, hence, disparate liquid-matrix interactions are absent. Findings for the magnitude and the range of wall effects on molecular positions and orientations and on molecular dynamics, including vibrational motion and structural relaxation, are reviewed. Moreover, their dependence on the parameters of the confinement and their relevance to theoretical approaches to the glass transition are addressed.

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.

Experimental demonstration of the novel "van-Hove integral method (vHI)" for measuring diffusive dynamics by elastic neutron scattering

Quasi-elastic neutron scattering (QENS)-based on the seminal work of Nobel Laureate Brockhouse-has been one of the major methods for studying pico-second to nano-second diffusive dynamics over the past 70 years. This is regarded as an "inelastic" method for dynamics. In contrast, we recently proposed a new neutron-scattering method for dynamics, which uses the elastic line of the scattering to access system dynamics directly in the time domain (Benedetto and Kearley in Sci Rep 9:11284, 2019). This new method has been denoted "vHI" that stands for "van Hove Integral". The reason is that, under certain conditions, the measured elastic intensity corresponds to the running-time integral of the intermediate scattering function, [Formula: see text], up to a time that is inversely proportional to the energy band-width incident on the sample. As a result, [Formula: see text] is accessed from the time derivative of the measured vHI profile. vHI has been supported by numerical and Monte-Carlo simulations, but has been difficult to validate experimentally due to the lack of a suitable instrument. Here we show that vHI works in practice, which we achieved by using a simple modification to the standard QENS backscattering spectrometer methodology. Basically, we varied the neutron-energy band-widths incident at the sample via a step-wise variation of the frequency of the monochromator Doppler-drive. This provides a measurement of the vHI profile at the detectors. The same instrument and sample were also used in standard QENS mode for comparison. The intermediate scattering functions, [Formula: see text], obtained by the two methods-vHI and QENS-are strikingly similar providing a direct experimental validation of the vHI method. Perhaps surprisingly, the counting statistics of the two methods are comparable even though the instrument used was expressly designed for QENS. This shows that the methodology modification adopted here can be used in practice to access vHI profiles at many of the backscattering spectrometers worldwide. We also show that partial integrations of the measured QENS spectrum cannot provide the vHI profile, which clarifies a common misconception. At the same time, we show a novel approach which does access [Formula: see text] from QENS spectra.

Why Should Metformin Not Be Given in Advanced Kidney Disease? Potential Leads from Computer Simulations

Metformin is considered as the go-to drug in the treatment of diabetes. However, it is either prescribed in lower doses or not prescribed at all to patients with kidney problems. To find a potential explanation for this practice, we employed atomistic-level computer simulations to simulate the transport of metformin through multidrug and toxin extrusion 1 (MATE1), a protein known to play a key role in the expulsion of metformin into urine. Herein, we examine the hydrogen bonding between MATE1 and one or more metformin molecules. The simulation results indicate that metformin continuously forms and breaks off hydrogen bonds with MATE1 residues. However, the mean hydrogen bond lifetimes increase for an order of magnitude when three metformin molecules are inserted instead of one. This new insight into the metformin transport process may provide the molecular foundation behind the clinical practice of not prescribing metformin to kidney disease patients.

Mechanism of deep eutectic solvents enhancing catalytic function of cytochrome P450 enzymes in biosynthesis and organic synthesis

Indigo is an insoluble blue dye, which generates serious pollution in its production process. Increasing focus has come to the biosynthesis of indigo that are more environment-preserved and high-efficient. Hence, this study was designed to explore the specific role of various deep eutectic solvents (DESs) on cytochromeP45-BM-3 catalyzing indole to produce indigo. DESs were synthesized by heating and stirring. The structure of the solvent was analyzed by nuclear magnetic resonance (NMR) and fourier transform infrared spectrum (FT-IR), and the relationship between the viscosity, density and refractive index of the solvent, and the water content of the solvent was examined. Circular dichroism spectrometer was used to detect the tertiary structure of the enzyme protein. The effect of solvent type, concentration, pH, temperature, and water content on the catalytic activity and stability of P450 BM-3 was measured using an ultraviolet spectrophotometer. A new solvent biphasic system was established using DESs and buffers, and indigo was prepared using recombinant E. coli-biocatalyzed indole. DESs were low-melting eutectics formed by molecules interaction of components through hydrogen bonding. The physical properties of DESs such as density, viscosity, and refractive index varied with water content and temperature of the solvent. The pH, water content, and temperature of DESs were positively correlated with the catalytic activity of P450 BM-3. To sum up, DESs can improve the catalytic activity and thermal stability of P450 BM-3. Indigo can be efficiently prepared using the DESs-buffer biphasic system.

Mechanisms of action of ionic liquids on living cells: the state of the art

Ionic liquids (ILs) are a relatively new class of organic electrolytes composed of an organic cation and either an organic or inorganic anion, whose melting temperature falls around room-temperature. In the last 20 years, the toxicity of ILs towards cells and micro-organisms has been heavily investigated with the main aim to assess the risks associated with their potential use in (industrial) applications, and to develop strategies to design greener ILs. Toxicity, however, is synonym with affinity, and this has stimulated, in turn, a series of biophysical and chemical-physical investigations as well as few biochemical studies focused on the mechanisms of action (MoAs) of ILs, key step in the development of applications in bio-nanomedicine and bio-nanotechnology. This review has the intent to present an overview of the state of the art of the MoAs of ILs, which have been the focus of a limited number of studies but still sufficient enough to provide a first glimpse on the subject. The overall picture that emerges is quite intriguing and shows that ILs interact with cells in a variety of different mechanisms, including alteration of lipid distribution and cell membrane viscoelasticity, disruption of cell and nuclear membranes, mitochondrial permeabilization and dysfunction, generation of reactive oxygen species, chloroplast damage (in plants), alteration of transmembrane and cytoplasmatic proteins/enzyme functions, alteration of signaling pathways, and DNA fragmentation. Together with our earlier review work on the biophysics and chemical-physics of IL-cell membrane interactions (Biophys. Rev. 9:309, 2017), we hope that the present review, focused instead on the biochemical aspects, will stimulate a series of new investigations and discoveries in the still new and interdisciplinary field of "ILs, biomolecules, and cells."

From protein and its hydration water dynamics to controlling mechano-elasticity of cellular lipid membranes and cell migration via ionic liquids

In this invited Commentary, as requested, I will walk the reader through my research path starting from my first works on proteins and their hydration water dynamics to my most recent activity on the use of ionic liquids (ILs) as molecular handles to control and manipulate cell membrane mechano-elasticity and cell migration. In doing so I will comment on my research activity on polymers, proteins, natural bioprotectants, phospholipid bilayers, amyloids and cells, which I have carried out by combining several different experimental and computational approaches including neutron scattering, atomic force microscopy, classical molecular dynamics and ab initio calculations, used in tandem with several biological assays and a palette of complementary techniques ranging from calorimetry to static and dynamic light scattering. In parallel to this biophysical/chemical-physical core activity, a smaller portion of my interest and effort has been-I may say always-dedicated to the development of a new neutron scattering method/spectroscopy for dynamics based on "elastic" scattering. I will comment on this instrumental side of my research as well and show its relationship with the biophysical core of my activity. The overall picture that emerges is, from my personal prospective, of a coherent 13-year research path based on curiosity and a problem-solving approach, in which the fundamental importance of inter- and trans-disciplinary approaches and collaborations is emerging on the way, forecasting a prosper and intriguing future for physics in biology and in nanomedicine and bionanotechnology applications.

A Quantitative Comparison of the Counting Significance of van Hove Integral Spectroscopy and Quasielastic Neutron Scattering

We have recently proposed a new method to access system dynamics via neutron scattering based on measuring the elastic scattered intensity: By varying the energy band-width that impinges on the sample (also known as instrumental energy resolution), the purely elastic-scattering from this variation is the running time-integral of the intermediate scattering function (I(t)) [Benedetto and Kearley, Sci. Rep. 9, 11284, 2019]. In this correspondence we denote our method "vHI", which stands for "van Hove Integral". The method is now widely accepted as "valid" and here we focus on the efficiency of the vHI method compared with the standard quasi-elastic neutron scattering (QENS) method. We use a numerical Monte-Carlo simulation of an instrument that is equally capable of measuring QENS and vHI under identical conditions. For an "experiment" in which the same number of neutrons enter the instrument, we present comparisons between QENS and vHI at three levels of data-reduction. Firstly, at the raw-data level vHI achieves 100 times more neutrons at the detector than QENS. Secondly, vHI has a factor of 2 less statistical error, which would translate to an overall gain of 4 for vHI in counting-time. Lastly, we compare the distortions caused in obtaining the final I(t) via time-Fourier transform (QENS) and polynomial time-derivative (vHI). Here, the statistical error is 10 times smaller for vHI. This last comparison is the most important result where the 10 times smaller residual for vHI gives a net gain in counting time of 100 better than QENS to obtain the same underlying dynamics of the system under study.

Quantification of spatio-temporal scales of dynamical heterogeneity of water near lipid membranes above supercooling

A hydrated 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid membrane is investigated using an all atom molecular dynamics simulation at 308 K to determine the physical sources of universal slow relaxations of hydration layers and length-scale of the spatially heterogeneous dynamics. Continuously residing interface water (IW) molecules hydrogen bonded to different moieties of lipid heads in the membrane are identified. The non-Gaussian parameters of all classes of IW molecules show a cross-over from cage vibration to translational diffusion. A significant non-Gaussianity is observed for the IW molecules exhibiting large length correlations in translational van Hove functions. Two time-scales for the ballistic motions and hopping transitions are obtained from the self intermediate scattering functions of the IW molecules with an additional long relaxation, which disappears for bulk water. The long relaxation time-scales for the IW molecules obtained from the self intermediate scattering functions are in good accordance with the hydrogen bond relaxation time-scales irrespective of the nature of the chemical confinement and the confinement lifetime. Employing a block analysis approach, the length-scale of dynamical heterogeneities is captured from a transition from non-Gaussianity to Gaussianity in van Hove correlation functions of the IW molecules. The heterogeneity length-scale is comparable to the wave-length of the small and weak undulations of the membrane calculated by Fourier transforms of lipid tilts. This opens up a new avenue towards a possible correlation between heterogeneity length-scale and membrane curvature more significant for rippled membranes. Thus, our analyses provide a measure towards the spatio-temporal scale of dynamical heterogeneity of confined water near membranes.

Strong influence of weak hydrogen bonding on actinide-phosphonate complexation: accurate predictions from DFT followed by experimental validation

Among the varied classes of weak hydrogen bond, the CHO type is one of immense interest as it governs the finer structures of biological and chemical molecules, hence determining their functionalities. In the present work, this weak hydrogen bond has been shown to strongly influence the complexation behaviour of uranyl nitrate [UO2(NO3)2] with diamyl-H-phosphonate (DAHP) and its branched isomer disecamyl-H-phosphonate (DsAHP). The structures of the bare ligands and complexes have been optimized by density functional theory (DFT) calculations. Surprisingly, despite having the same chemical composition the branched UO2(NO3)2·2DsAHP complex shows a remarkably higher stability (by ∼14 kcal mol-1) compared to the UO2(NO3)2·2DAHP complex. Careful inspection of the optimized structures reveals the existence of multiple CHO hydrogen-bonding interactions between the nitrate oxygens or U[double bond, length as m-dash]O oxygens and the α-hydrogens in the alkyl chains of the ligands. Comparatively stronger such bonds are found in the UO2(NO3)2·2DsAHP complex. The binding free energies associated with the complexes are computed and favoured superior binding energetics for the more stable UO2(NO3)2·2DsAHP complex. Calculations involving diisoamyl-H-phosphonate (DiAHP) and its complexes have also been performed. Theoretical predictions are experimentally tested by carrying out the extraction of U(vi) from nitric acid media using these ligands. DAHP, DsAHP and DiAHP are synthesised, characterised by NMR and evaluated for their physicochemical properties viz. viscosity, density and aqueous solubility. It was experimentally discovered that indeed DsAHP complexation with uranyl nitrate is more favoured. H-phosphonates are generically classified as acidic extractants owing to the formation of an enol tautomer at lower acidities, hence complexing the metal ion by proton exchange. Our experiments interestingly reveal a neutral ligand characteristic for DsAHP alone which is generically an acidic extractant. Furthermore, the enol tautomer of H-phosphonates that governs their extraction profiles at low acidities is also explored by DFT and the anomalous pH dependent complexation trend of DsAHP could be successfully explained. The extractions of Pu(iv) and Th(iv) have also been carried out in addition to U(vi). Solvent extraction behaviour of Am(iii) was also studied with all three ligands; the positive binding energies computed for the Am(iii) complexation corroborate with our experimental results on the poor extraction of Am(iii).

Water-Polymer Coupling Induces a Dynamical Transition in Microgels

The long debated protein dynamical transition was recently found also in nonbiological macromolecules, such as poly- N-isopropylacrylamide (PNIPAM) microgels. Here, by using atomistic molecular dynamics simulations, we report a description of the molecular origin of the dynamical transition in these systems. We show that PNIPAM and water dynamics below the dynamical transition temperature T _d are dominated by methyl group rotations and hydrogen bonding, respectively. By comparing with bulk water, we unambiguously identify PNIPAM-water hydrogen bonding as mainly responsible for the occurrence of the transition. The observed phenomenology thus crucially depends on the water-macromolecule coupling, being relevant to a wide class of hydrated systems, independently from the biological function.

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.

Correlating interfacial water dynamics with protein-protein interaction in complex of GDF-5 and BMPRI receptors

GDF-5 mediated signal transduction regulating chondrogenesis and skeletogenesis involves three different type-I receptors viz. Act-RI, BMPRIA and BMPRIB. BMPRIA and BMPRIB generally shows temporal and spatial co-expression but some spatially different expression pattern has also been observed. BMPRIA receptor is the key receptor implicated in BMP signalling during osteogenesis and is expressed in osteoblasts during the course of bone formation. However, BMPRIB appears to be primarily expressed in mesenchymal pre-cartilage condensations and also found in differentiated osteoblast and chondrocytes. The extracellular pH affects bone cell function and it is experimentally known that mineralization of bone is affected by shift of pH in cultured osteoblast. Here we report the effect of pH on dynamics of water present at the interface of GDF-5:BMPRIA and GDF-5:BMPRIB and binding interaction energy of these complexes. Water dynamics at different pH was analysed using residence time and hydrogen bond relaxation kinetics. pH influences the interaction energy between GDF-5 and BMPRIA and BMPRIB receptors indicating the electrostatic environment modulating the activity of two receptors. This pH dependence of interaction energy is further supported by similar behaviour of hydrogen bond existence of buried water molecules at the interface. In contrast to this the slow and fast exchanging water molecules do not show similar pH dependence of hydrogen bonding relaxation kinetics. Hence; we conclude that only buried water molecule at the interface influences the protein-protein interaction and the electrostatic environment of the extracellular fluid might decide the specificity of the two receptors.

Low-Temperature Decoupling of Water and Protein Dynamics Measured by Neutron Scattering

Water plays a major role in biosystems, greatly contributing to determine their structure, stability, and function. It is well known, for instance, that proteins require a minimum amount of water to be fully functional. Despite many years of intensive research, however, the detailed nature of protein-hydration water interactions is still partly unknown. The widely accepted "protein dynamical transition" scenario is based on perfect coupling between the dynamics of proteins and that of their hydration water, which has never been probed in depth experimentally. I present here high-resolution elastic neutron scattering measurements of the atomistic dynamics of lysozyme in water. The results show for the first time that the dynamics of proteins and of their hydration water are actually decoupled at low temperatures. This important result challenges the "protein dynamical transition" scenario and requires a new model to link protein dynamics to the dynamics of its hydration water.