Abstact

The solution structure of Desulfovibrio vulgaris Hildenborough (DvH) ferrocytochrome c553 has been determined by nuclear magnetic resonance spectroscopy and combined simulated annealing/high temperature restrained molecular dynamics calculations. This three-stage protocol consists of an initial determination of overall fold from randomised co-ordinates, followed by a 20 picosecond exploratory stage, during which the non-bonded terms are simplified to facilitate as broad a sampling of conformational space as possible, and a 26 picosecond refinement stage, using the full AMBER force field. This latter stage systematically improved the energetic and convergence characteristics of the ensemble, while still satisfying the experimental restraints. Forty structures have been obtained from a total of 875 distance constraints for this protein of 79 amino acid residues. The root-mean-square deviation over all residues with respect to the mean is 0.70(+/- 0.12)A for the backbone (N, C alpha and C') atoms. Two conformations of the turn motif at the solvent/heme cleft interface have been identified, both fulfilling the experimental data and having equally viable energetic characteristics. The stability of the ensemble and the dynamic characteristics have been further investigated by subjecting ten of the structures to constraint-free molecular dynamics calculations (130 picoseconds) in vacuo. The structures were found to be stable to within 1.5 A of the initial backbone conformation. Comparison with the dynamic behaviour of the restrained molecular dynamics calculations has been used to identify regions of inherent flexibility in the molecule.