Molecular dynamics simulation of the neuroglobin crystal: Comparison with the simulation in solution

Neuroglobin (Ngb) is a monomeric protein that, despite the small sequence similarity with other globins, displays the typical globin fold. In the absence of exogenous ligands, the ferric and the ferrous forms of Ngb are both hexacoordinated to the distal and proximal histidines. In the ferrous form, oxygen, nitric oxide or carbon monoxide can displace the distal histidine, yielding a reversible adduct. Crystallographic data show that the binding of an exogenous ligand is associated to structural changes involving heme sliding and a topological reorganization of the internal cavities. Molecular dynamics (MD) simulations in solution show that the heme oscillates between two positions, much as the ones observed in the crystal structure, although the occupancy is different. The simulations also suggest that ligand binding in solution can affect the flexibility and conformation of residues connecting the C and D helices, referred to as the CD corner, which is coupled to the configuration adopted by the distal histidine. In this study, we report the results of 30 ns MD simulations of CO-bound Ngb in the crystal. Our goal was to compare the protein dynamical behavior in the crystal with the results supplied by the previous MD simulation of CO-bound Ngb in solution and the x-ray experimental data. The results show that the different environments (crystal or solution) affect the dynamics of the heme group and of the CD corner. © 2008 by the Biophysical Society.