Practical applications of time-averaged restrained molecular dynamics to ligand-receptor systems: FK506 bound to the Q50R,A95H,K98I triple mutant of FKBP-13.

The ability of time-averaged restrained molecular dynamics (TARMD) to escape local low-energy conformations and explore conformational space is compared with conventional simulated-annealing methods. Practical suggestions are offered for performing TARMD calculations with ligand-receptor systems, and are illustrated for the complex of the immunosuppressant FK506 bound to Q50R,A95H,K98I triple mutant FKBP-13. The structure of (13)C-labeled FK506 bound to triple-mutant FKBP-13 was determined using a set of 87 NOE distance restraints derived from HSQC-NOESY experiments. TARMD was found to be superior to conventional simulated-annealing methods, and produced structures that were conformationally similar to FK506 bound to wild-type FKBP-12. The individual and combined effects of varying the NOE restraint force constant, using an explicit model for the protein binding pocket, and starting the calculations from different ligand conformations were explored in detail.