3 Nsec molecular dynamics simulation of the protein ubiquitin and comparison with X-ray crystal and solution NMR structures

Mainly due to computational limitations, past protein molecular dynamics simulations have rarely been extended to 300 psec; we are not aware of any published results beyond 350 psec. The present work compares a 3000 psec simulation of the protein ubiquitin with the available x-ray crystallographic and solution NMR structures. Aside from experimental structure availability, ubiquitin was studied because of its relatively small size (76 amino acids) and lack of disulfide bridges. An implicit solvent model was used except for explicit treatment of waters of crystallization. We found that the simulated average structure retains most of the character of the starting x-ray crystal structure. In two highly surface accessible regions, the simulation was not in agreement with the x-ray structure. In addition, there are six backbone-backbone hydrogen bonds that are in conflict between the solution NMR and x-ray crystallographic structures; two are bonds that the NMR does not locate, and four are ones that the two methods disagree upon the donor. Concerning these six backbone-backbone hydrogen bonds, the present simulation agrees with the solution NMR structure in five out-of-the six cases, in that if a hydrogen bond is present in the x-ray structure and not in the NMR structure, the bond breaks within 700 psec. Of the two hydrogen bonds that are found in the NMR structure and not in the x-ray structure, one forms at 1400 psec and the other forms rarely. The present results suggest that relatively long molecular dynamics simulations, that use protein x-ray crystal coordinates for the starting structure and a computationally efficient solvent representation, may be used to gain an understanding of conformational and dynamic differences between the solid-crystal and dilute-solution states.