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

Crystal structure of a DNA decamer showing a novel pseudo four-way helix-helix junction.

The crystal structure of the decanucleotide d(CGCAATTGCG)2 has been solved by a combination of molecular replacement and heavy-atom procedures and has been refined to an R factor of 20.2% at 2.7 A. It is not a fully base-paired duplex but has a central core of eight Watson-Crick base pairs flanked by unpaired terminal guanosines and cytosines. These participate in hydrogen-bonding arrangements with adjacent decamer duplexes in the crystal lattice. The unpaired guanosines are bound in the G+C regions of duplex minor grooves. The cytosines have relatively high mobility, even though they are constrained to be in one region where they are involved in base-paired triplets with G.C base pairs. The 5'-AATT sequence in the duplex region has a narrow minor groove, providing further confirmation of the sequence-dependent nature of groove width

The Fe-CO Bond Energy in Myoglobin: A QM/MM Study of the Effect of Tertiary Structure

The Fe-CO bond dissociation energy (BDE) in myoglobin (Mb) has been calculated with B3LYP quantum mechanics/molecular mechanics methods for 22 different Mb conformations, generated from molecular dynamics simulations. Our average BDE of 8.1 kcal/mol agrees well with experiment and shows that Mb weakens the Fe-CO bond by 5.8 kcal/mol; the calculations provide detailed atomistic insight into the origin of this effect. BDEs for Mb conformations with the R carbonmonoxy tertiary structure are on average 2.6 kcal/mol larger than those with the T deoxy tertiary structure, suggesting two functionally distinct allosteric states. This allostery is partly explained by the reduction in distal cavity steric crowding as Mb moves from its T to R tertiary structure