Hydration dynamics of protein molecules in aqueous solution: Unity among diversity

Dielectric dispersion and NMRD experiments have revealed that a significant fraction of water molecules in the hydration shell of various proteins do not exhibit any slowing down of dynamics. This is usually attributed to the presence of the hydrophobic residues (HBR) on the surface, although HBRs alone cannot account for the large amplitude of the fast component. Solvation dynamics experiments and also computer simulation studies, on the other hand, repeatedly observed the presence of a non-negligible slow component. Here we show, by considering three well-known proteins (lysozyme, myoglobin and adelynate kinase), that the fast component arises partly from the response of those water molecules that are hydrogen bonded with the backbone oxygen (BBO) atoms. These are structurally and energetically less stable than those with the side chain oxygen (SCO) atoms. In addition, the electrostatic interaction energy distribution (EIED) of individual water molecules (hydrogen bonded to SCO) with side chain oxygen atoms shows a surprising two peak character with the lower energy peak almost coincident with the energy distribution of water hydrogen bonded to backbone oxygen atoms (BBO). This two peak contribution appears to be quite general as we find it for lysozyme, myoglobin and adenylate kinase (ADK). The sharp peak of EIED at small energy (at less than 2 k(B)T) for the BBO atoms, together with the first peak of EIED of SCO and the HBRs on the protein surface, explain why a large fraction (similar to 80%) of water in the protein hydration layer remains almost as mobile as bulk water Significant slowness arises only from the hydrogen bonds that populate the second peak of EIED at larger energy (at about 4 k(B)T). Thus, if we consider hydrogen bond interaction alone, only 15-20% of water molecules in the protein hydration layer can exhibit slow dynamics, resulting in an average relaxation time of about 5-10 ps. The latter estimate assumes a time constant of 20-100 ps for the slow component. Interestingly, relaxation of water molecules hydrogen bonded to back bone oxygen exhibit an initial component faster than the bulk, suggesting that hydrogen bonding of these water molecules remains frustrated. This explanation of the heterogeneous and non-exponential dynamics of water in the hydration layer is quantitatively consistent with all the available experimental results, and provides unification among diverse features.

Stabilized SPH-based simulations of impact dynamics using acceleration-corrected artificial viscosity

Artificial viscosity in SPH-based computations of impact dynamics is a numerical artifice that helps stabilize spurious oscillations near the shock fronts and requires certain user-defined parameters. Improper choice of these parameters may lead to spurious entropy generation within the discretized system and make it over-dissipative. This is of particular concern in impact mechanics problems wherein the transient structural response may depend sensitively on the transfer of momentum and kinetic energy due to impact. In order to address this difficulty, an acceleration correction algorithm was proposed in Shaw and Reid (''Heuristic acceleration correction algorithm for use in SPH computations in impact mechanics'', Comput. Methods Appl. Mech. Engrg., 198, 3962-3974) and further rationalized in Shaw et al. (An Optimally Corrected Form of Acceleration Correction Algorithm within SPH-based Simulations of Solid Mechanics, submitted to Comput. Methods Appl. Mech. Engrg). It was shown that the acceleration correction algorithm removes spurious high frequency oscillations in the computed response whilst retaining the stabilizing characteristics of the artificial viscosity in the presence of shocks and layers with sharp gradients. In this paper, we aim at gathering further insights into the acceleration correction algorithm by further exploring its application to problems related to impact dynamics. The numerical evidence in this work thus establishes that, together with the acceleration correction algorithm, SPH can be used as an accurate and efficient tool in dynamic, inelastic structural mechanics. (C) 2011 Elsevier Ltd. All rights reserved.

Anomalous structural and dynamical phase transitions of soft colloidal binary mixtures

We report microscopic structural and dynamical measurements on binary mixtures of homopolymers and polymer grafted nanoparticles at high densities in good solvent. We find strong and unexpected anomalies in the structure and dynamics of these binary mixtures, including appearance of spontaneous orientational alignment, as a function of added homopolymers of different molecular weights. Our experiments point to the possibility of exploiting the phase space in density and homopolymer size, of such hybrid systems, to create new materials with novel structural and physical properties.

Diffusion of pentane isomers in faujasite-type zeolites : NMR and molecular dynamics study

Diffusion of pentane isomers in zeolites NaX has been investigated using pulsed field gradient nuclear magnetic resonance (PFG-NMR) and molecular dynamics (MD) techniques respectively. Temperature and concentration dependence of diffusivities have been studied. The diffusivities obtained from NMR are roughly an order of magnitude smaller than those obtained from MD. The dependence of diffusivity on loading at high temperatures exhibits a type I behavior according to the classification of Karger and Pfeifer 1]. NMR diffusivities of the isomers exhibit the order D(n-pentane) > D(isopentane) > D(neopentane). The results from MD suggest that the diffusivities of the isomers follow the order D(n-pentane) < D(isopentane) < D(neopentane). The activation energies from NMR show E-a(n-pentane) < E-a(isopentane) < E-a(neopentane) whereas those from MD suggest the order E-a(n-pentane) > (isopentane) > E-a(neopentane). The latter follows the predictions of levitation effect whereas those of NMR appears to be due to the presence of defects in the zeolite crystals. The differences between diffusivities estimated by NMR and MD are attributed to the longer time and length scales sampled by the NMR technique, as compared to MD. (C) 2012 Elsevier Inc. All rights reserved.

Gap-Dependent Quasiparticle Dynamics and Coherent Acoustic Phonons in CaFe2As2 across Spin Density Wave Phase Transition

We report ultrafast quasiparticle (QP) dynamics and coherent acoustic phonons in undoped CaFe2As2 iron pnictide single crystals exhibiting spin-density wave (SDW) and concurrent structural phase transition at temperature T-SDW similar to 165K using femtosecond time-resolved pump-probe spectroscopy. The contributions in transient differential reflectivity arising from exponentially decaying QP relaxation and oscillatory coherent acoustic phonon mode show large variations in the vicinity of T-SDW. From the temperature-dependence of the QP recombination dynamics in the SDW phase, we evaluate a BCS-like temperature dependent charge gap with its zero-temperature value of similar to(1.6 perpendicular to 0.2)k(B)T(SDW), whereas, much above T-SDW, an electron-phonon coupling constant of similar to 0.13 has been estimated from the linear temperature-dependence of the QP relaxation time. The long-wavelength coherent acoustic phonons with typical time-period of similar to 100 ps have been analyzed in the light of propagating strain pulse model providing important results for the optical constants, sounds velocity and the elastic modulus of the crystal in the whole temperature range of 3 to 300 K.

Solvation dynamics of tryptophan in water-dimethyl sulfoxide binary mixture: In search of molecular origin of composition dependent multiple anomalies

Experimental and simulation studies have uncovered at least two anomalous concentration regimes in water-dimethyl sulfoxide (DMSO) binary mixture whose precise origin has remained a subject of debate. In order to facilitate time domain experimental investigation of the dynamics of such binary mixtures, we explore strength or extent of influence of these anomalies in dipolar solvation dynamics by carrying out long molecular dynamics simulations over a wide range of DMSO concentration. The solvation time correlation function so calculated indeed displays strong composition dependent anomalies, reflected in pronounced non-exponential kinetics and non-monotonous composition dependence of the average solvation time constant. In particular, we find remarkable slow-down in the solvation dynamics around 10%-20% and 35%-50% mole percentage. We investigate microscopic origin of these two anomalies. The population distribution analyses of different structural morphology elucidate that these two slowing down are reflections of intriguing structural transformations in water-DMSO mixture. The structural transformations themselves can be explained in terms of a change in the relative coordination number of DMSO and water molecules, from 1DMSO:2H(2)O to 1H(2)O:1DMSO and 1H(2)O:2DMSO complex formation. Thus, while the emergence of first slow down (at 15% DMSO mole percentage) is due to the percolation among DMSO molecules supported by the water molecules (whose percolating network remains largely unaffected), the 2nd anomaly (centered on 40%-50%) is due to the formation of the network structure where the unit of 1DMSO:1H(2)O and 2DMSO:1H(2)O dominates to give rise to rich dynamical features. Through an analysis of partial solvation dynamics an interesting negative cross-correlation between water and DMSO is observed that makes an important contribution to relaxation at intermediate to longer times.

Relation Between the Diffusivity, Viscosity, and Ionic Radius of LiCl in Water, Methanol, and Ethylene Glycol: A Molecular Dynamics Simulation

A molecular dynamics (MD) investigation of LiCl in water, methanol, and ethylene glycol (EG) at 298 K is reported. Several; structural and dynamical properties of the ions as well as the solvent such as self-diffusivity, radial distribution functions, void and neck distributions, velocity autocorrelation functions, and mean residence times of solvent in the first solvation shell have been computed. The results show that the reciprocal relationship between the self-diffusivity of the ions and the viscosity is valid in almost all solvents with the exception of water. From an analysis of radial distribution functions and coordination numbers the nature of hydrogen bonding within the solvent and its influence on the void and neck distribution becomes evident. It is seen that the solvent solvent interaction is important in EG while solute solvent interactions dominate in water and methanol. From Voronoi tessellation, it is seen that the voids and necks within methanol are larger as compared to those within water or EG. On the basis of the void and neck distributions obtained from MD simulations and literature experimental data of limiting ion conductivity for various ions of different sizes we show that there is a relation between the void and neck radius on e one hand and dependence of conductivity on the ionic radius on the other. It is shown that the presence of large diameter voids and necks in methanol is responsible for maximum in limiting ion conductivity (lambda(0)) of TMA(+), while in water in EG, the maximum is seen for Rb+. In the case of monovalent anions, maximum in lambda(0) as a function ionic radius is seen for Br- in water EG but for the larger ClO4- ion in methanol. The relation between the void and neck distribution and the variation in lambda(0) with ionic radius arises via the Levitation effect which is discussed. These studies show the importance of the solvent structure and the associated void structure.

Layerwise decomposition of water dynamics in reverse micelles: A simulation study of two-dimensional infrared spectrum

We present computer simulation study of two-dimensional infrared spectroscopy (2D-IR) of water confined in reverse micelles (RMs) of various sizes. The present study is motivated by the need to understand the altered dynamics of confined water by performing layerwise decomposition of water, with an aim to quantify the relative contributions of different layers water molecules to the calculated 2D-IR spectrum. The 0-1 transition spectra clearly show substantial elongation, due to in-homogeneous broadening and incomplete spectral diffusion, along the diagonal in the surface water layer of different sized RMs. Fitting of the frequency fluctuation correlation functions reveal that the motion of the surface water molecules is sub-diffusive and indicate the constrained nature of their dynamics. This is further supported by two peak nature of the angular analogue of van Hove correlation function. With increasing system size, the water molecules become more diffusive in nature and spectral diffusion almost completes in the central layer of the larger size RMs. Comparisons between experiments and simulations establish the correspondence between the spectral decomposition available in experiments with the spatial decomposition available in simulations. Simulations also allow a quantitative exploration of the relative role of water, sodium ions, and sulfonate head groups in vibrational dephasing. Interestingly, the negative cross correlation between force on oxygen and hydrogen of O-H bond in bulk water significantly decreases in the surface layer of each RM. This negative cross correlation gradually increases in the central water pool with increasing RMs size and this is found to be partly responsible for the faster relaxation rate of water in the central pool. (C) 2013 AIP Publishing LLC.

Positive temperature coefficient and structural relaxations in selectively localized MWNTs in PE/PEO blends

The dispersion state of multiwall carbon nanotubes (MWNTs) in melt mixed polyethylene/polyethylene oxide (PE/PEO) blends has been assessed by both surface and volume electrical conductivity measurements and the structural relaxations have been assessed by broadband dielectric spectroscopy. The selective localization of MWNTs in the blends was controlled by the flow characteristics of the components, which led to their localization in the energetically less favored phase (PE). The electrical conductivity and positive temperature co-efficient (PTC) measurements were carried out on hot pressed samples. The neat blends exhibited only a negative temperature coefficient (NTC) effect while the blends with MWNTs exhibited both a PTC and a NTC at the melting temperatures of PE and PEO respectively. These phenomenal changes were corroborated with the different crystalline morphology in the blends. It was deduced that during compression molding, the more viscous PEO phase spreads less in contrast to the less viscous PE phase. This has further resulted in a gradient in morphology as well as the distribution state of the MWNTs in the samples and was supported by scanning electron and scanning acoustic microscopy (SAM) studies and contact angle measurements. SAM from different depths of the samples revealed a gradient in the microstructure in the PE/PEO blends which is contingent upon the flow characteristics of the components. Interestingly, the surface and volume electrical conductivity was different due to the different dispersion state of the MWNTs at the surface and bulk. The observed surface and volume electrical conductivity measurements were corroborated with the evolved morphology during processing. The structural relaxations in both PE and PEO were discerned from broadband dielectric spectroscopy. The segmental dynamics below and above the melting temperature of PEO were significantly different in the presence of MWNTs.

mulPBA: an efficient multiple protein structure alignment method based on a structural alphabet

The increasing number of available protein structures requires efficient tools for multiple structure comparison. Indeed, multiple structural alignments are essential for the analysis of function, evolution and architecture of protein structures. For this purpose, we proposed a new web server called multiple Protein Block Alignment (mulPBA). This server implements a method based on a structural alphabet to describe the backbone conformation of a protein chain in terms of dihedral angles. This sequence-like' representation enables the use of powerful sequence alignment methods for primary structure comparison, followed by an iterative refinement of the structural superposition. This approach yields alignments superior to most of the rigid-body alignment methods and highly comparable with the flexible structure comparison approaches. We implement this method in a web server designed to do multiple structure superimpositions from a set of structures given by the user. Outputs are given as both sequence alignment and superposed 3D structures visualized directly by static images generated by PyMol or through a Jmol applet allowing dynamic interaction. Multiple global quality measures are given. Relatedness between structures is indicated by a distance dendogram. Superimposed structures in PDB format can be also downloaded, and the results are quickly obtained. mulPBA server can be accessed at www.dsimb.inserm.fr/dsimb_tools/mulpba/.

Molecular Dynamics Perspective on the Protein Thermal Stability: A Case Study Using SAICAR Synthetase

The enzyme SAICAR synthetase ligates aspartate with CAIR (5'-phosphoribosyl-4-carboxy-5-aminoimidazole) forming SAICAR (5-amino-4-imidazole-N-succinocarboxamide ribonucleotide) in the presence of ATP. In continuation with our previous study on the thermostability of this enzyme in hyper-/thermophiles based on the structural aspects, here, we present the dynamic aspects that differentiate the mesophilic (E. coli, E. chaffeensis), thermophilic (G. kaustophilus), and hyperthermophilic (M. jannaschii, P. horikoshii) SAICAR synthetases by carrying out a total of 11 simulations. The five functional dimers from the above organisms were simulated using molecular dynamics for a period of 50 ns each at 300 K, 363 K, and an additional simulation at 333 K for the thermophilic protein. The basic features like root-mean-square deviations, root-mean-square fluctuations, surface accessibility, and radius of gyration revealed the instability of mesophiles at 363 K. Mean square displacements establish the reduced flexibility of hyper-/thermophiles at all temperatures. At the simulations time scale considered here, the long-distance networks are considerably affected in mesophilic structures at 363 K. In mesophiles, a comparatively higher number of short-lived (having less percent existence time) C alpha, hydrogen bonds, hydrophobic interactions are formed, and long-lived (with higher percentage existence time) contacts are lost. The number of time-averaged salt-bridges is at least 2-fold higher in hyperthermophiles at 363 K. The change in surface accessibility of salt-bridges at 363 K from 300 K is nearly doubled in mesophilic protein compared to proteins from other temperature classes.

Cooperativity and Structural Relaxations in PVDF/PMMA Blends in the Presence of MWNTs: An Assessment through SAXS and Dielectric Spectroscopy

Intermolecular cooperativity and structural relaxations in PVDF/PMMA blends were studied in this work with respect to different surface modified (amine, similar to NH2; carboxyl acid, similar to COOH and pristine) multiwalled nanotubes (MWNTs) at 1 wt % near blend's T-g and in the vicinity of demixing using dielectric spectroscopy, SAXS, DSC, and WAXD. Intermolecular cooperativity at T-g and configurational entropy was addressed in the framework of cooperative rearranging region (CRR) at T-g. Because of specific interactions between PVDF and NH2-MWNTs, the local composition fluctuates at its average value resulting in a broad T-g. The scale of cooperativity (xi(CRR)) and the number of segments in the cooperative volume (N-CRR) is comparatively smaller in the blends with NH2-MWNTs. This clearly suggests that the number of segments cooperatively relaxing is reduced in the blends due to specific interactions leading to more heterogeneity. The configurational entropy at T-g, as derived from Vogel-Fulcher and Adam-Gibbs analysis, was reduced in the blends in presence of MWNTs manifesting in entropic penalty of the chains. The crystallite size and the amorphous miscibility was evaluated using SAXS and was observed to be strongly contingent on the surface functional groups on MWNTs. Three distinct relaxations-alpha(c) due to relaxations in the crystalline phase of PVDF, alpha(m) indicating the amorphous miscibility in PVDF/PMMA blends, and alpha beta concerning the segmental dynamics of PMMA-were observed in the blends in the temperature range T-g < T < T-c. The dynamics as well as the nature of relaxations were observed to be dependent the surface functionality on the MWNTs. The dielectric permittivity was also enhanced in presence of MWNTs, especially with NH2-MWNTs, with minimal losses. The influence of the MWNTs on the spherulite size and crystalline morphology of the blends was also confirmed by POM and SEM.

Sensitivity of Water Dynamics to Biologically Significant Surfaces of Monomeric Insulin: Role of Topology and Electrostatic Interactions

In addition to the biologically active monomer of the protein insulin circulating in human blood, the molecule also exists in dimeric and hexameric forms that are used as storage. The insulin monomer contains two distinct surfaces, namely, the dimer forming surface (DFS) and the hexamer forming surface (HFS), that are specifically designed to facilitate the formation of the dimer and the hexamer, respectively. In order to characterize the structural and dynamical behavior of interfacial water molecules near these two surfaces (DFS and HFS), we performed atomistic molecular dynamics simulations of insulin with explicit water. Dynamical characterization reveals that the structural relaxation of the hydrogen bonds formed between the residues of DFS and the interfacial water molecules is faster than those formed between water and that of the HFS. Furthermore, the residence times of water molecules in the protein hydration layer for both the DFS and HFS are found to be significantly higher than those for some of the other proteins studied so far, such as HP-36 and lysozyme. In particular, we find that more structured water molecules, with higher residence times (similar to 300-500 ps), are present near HFS than those near DFS. A significant slowing down is observed in the decay of associated rotational auto time correlation functions of O-H bond vector of water in the vicinity of HFS. The surface topography and the arrangement of amino acid residues work together to organize the water molecules in the hydration layer in order to provide them with a preferred orientation. HFS having a large polar solvent accessible surface area and a convex extensive nonpolar region, drives the surrounding water molecules to acquire predominantly an outward H-atoms directed, clathrate-like structure. In contrast, near the DFS, the surrounding water molecules acquire an inward H-atoms directed orientation owing to the flat curvature of hydrophobic surface and the interrupted hydrophilic residual alignment. We have followed escape trajectory of several such quasi-bound water molecules from both the surfaces that reveal the significant differences between the two hydration layers.

Non-equilibrium segmental dynamics driven by multiwall carbon nanotubes in PS/PVME blends

The present paper discusses the effect of multiwall carbon nanotubes (MWNTs) on the structural relaxation and the intermolecular cooperativity in dynamically asymmetric blends of PS/PVME (polystyrene/poly(vinyl methyl ether)). The temperature regime where chain connectivity effects dominate the thermodynamic concentration fluctuation (T/T-g > 0.75, T-g is the glass transition temperature of the blends) was studied using dielectric spectroscopy (DS). Interestingly, in the blends with MWNTs a bimodal distribution of relaxation was obtained in the loss modulus spectra. This plausibly is due to different environments experienced by the faster component (PVME) in the presence of MWNTs. The segmental dynamics of PVME was observed to be significantly slowed down in the presence of MWNTs and an Arrhenius-type behavior, weakly dependent on temperature, is observed at higher frequencies. This non-equilibrium dynamics of PVME is presumed to be originating from interphase regions near the surface of MWNTs. The length scale of the cooperative rearranging region (xi CRR) at T-g, assessed by calorimetric measurements, was observed to be higher in the case of blends with MWNTs. An enhanced molecular level miscibility driven by MWNTs in the blends corroborates with the larger xi CRR and comparatively more number of segments in CRR (in contrast to neat blends) around T-g. The configurational entropy and length scale of the cooperative volume was mapped as a function of temperature in the temperature regime, Tg < T < T-g + 60 K. The blends phase separated by spinodal decomposition which further led to an interconnected PVME network in PS. This further led to materials with very high electrical conductivity upon demixing.

Anomalous power law decay in solvation dynamics of DNA: a mode coupling theory analysis of ion contribution

Several time dependent fluorescence Stokes shift (TDFSS) experiments have reported a slow power law decay in the hydration dynamics of a DNA molecule. Such a power law has neither been observed in computer simulations nor in some other TDFSS experiments. Here we observe that a slow decay may originate from collective ion contribution because in experiments DNA is immersed in a buffer solution, and also from groove bound water and lastly from DNA dynamics itself. In this work we first express the solvation time correlation function in terms of dynamic structure factors of the solution. We use mode coupling theory to calculate analytically the time dependence of collective ionic contribution. A power law decay in seen to originate from an interplay between long-range probe-ion direct correlation function and ion-ion dynamic structure factor. Although the power law decay is reminiscent of Debye-Falkenhagen effect, yet solvation dynamics is dominated by ion atmosphere relaxation times at longer length scales (small wave number) than in electrolyte friction. We further discuss why this power law may not originate from water motions which have been computed by molecular dynamics simulations. Finally, we propose several experiments to check the prediction of the present theoretical work.