Temperature dependence of protein dynamics simulated with three different water models

The effect of variation of the water model on the temperature dependence of protein and hydration water dynamics is examined by performing molecular dynamics simulations of myoglobin with the TIP3P, TIP4P, and TIP5P water models and the CHARMM protein force field at temperatures between 20 and 300 K. The atomic mean-square displacements, solvent reorientational relaxation times, pair angular correlations between surface water molecules, and time-averaged structures of the protein are all found to be similar, and the protein dynamical transition is described almost indistinguishably for the three water potentials. The results provide evidence that for some purposes changing the water model in protein simulations without a loss of accuracy may be possible

Liquid-like and solid-like motions in proteins

Recent analyses of molecular dynamics simulations of hydrated C-phycocyanin suggest that the internal single-particle dynamics of this protein can be decomposed into four almost decoupled motion types: (1) diffusion of residues ("beads") in an effective harmonic potential, (2) corresponding vibrations in a local potential well, (3) purely rotational rigid side-chain diffusion, and (4) residue deformations. Each residue bead is represented by the corresponding C[α] carbon atom on the main chain. The effective harmonic residue potential can be imagined as the envelope of many local wells which are separated by small energy barriers. The residue friction matrix is assumed diagonal and the individual friction constants can be related to the density of the surrounding atoms. In this article we show that our model can be applied to lysozyme in solution as well, the only difference being that the side-chain deformations are more important and seem to be strongly correlated with the side-chain rotations. Comparing the simulated coherent scattering function of C-phycocyanin to a neutron spin-echo spectrum we show that our model can also describe collective motions in proteins at the residue level

A molecular view of melting in anhydrous phospholipidic membranes

A high-flux backscattering spectrometer and a time-of-flight disk chopper spectrometer are used to probe the molecular mobility of model freeze-dried phospholipid liposomes at a range of temperatures surrounding the main melting transition. Using specific deuteration, quasielastic neutron scattering provides evidence that, in contrast to the hydrocarbon chains, the headgroups of the phospholipid molecules do not exhibit a sharp melting transition. The onset of motion in the tails is located at temperatures far below the calorimetric transition. Long-range motion is achieved through the onset of whole-lipid translation at the melting temperature. Atomistic simulations are performed on a multibilayer model at conditions corresponding to the scattering experiments. The model provides a good description of the dynamics of the system, with predictions of the scattering functions that agree with experimental results. The analysis of both experimental data and results of simulations supports a picture of a gradual melting of the heterogeneous hydrophobic domain, with part of the chains spanning increasingly larger volumes and part of them remaining effectively immobile until the thermodynamic phase transition occurs

Nanoscale Imaging of Buried Structures via Scanning Near-Field

16. Materials and methods are available as supportin

Slow Solvation Dynamics at the Active Site of an Enzyme: Implications for Catalysis

Solvation dynamics at the active site of an enzyme, glutaminyl-tRNA synthetase (GlnRS), was studied using a fluorescence probe, acrylodan, site-specifically attached at cysteine residue C229, near the active site. The picosecond time-dependent fluorescence Stokes shift indicates slow solvation dynamics at the active site of the enzyme, in the absence of any substrate. The solvation dynamics becomes still slower when the substrate (glutamine or tRNAGln) binds to the enzyme. A mutant Y211H-GlnRS was constructed in which the glutamine binding site is disrupted. The mutant Y211H-GlnRS labeled at C229 with acrylodan exhibited significantly different solvent relaxation, thus demonstrating that the slow dynamics is indeed associated with the active site. Implications for catalysis and specificity have been discussed