Molecular dynamics simulations of the docking of substituted N5-deazapterins to dihydrofolate reductase

Orientations of the deazapterin ring and the conformational preferences of groups appended to the deazapterin ring in a set of 8-substituted deazapterin cations docked into the dihydrofolate reductase (DHFR) binding site have been investigated using a methodology based on the simulated annealing technique within molecular dynamics (MD) simulations. Of five possible binding pockets for the 8-substituents, identified from a preliminary manual docking study, one has been definitively eliminated after an analysis of MD trajectories, while another remains uncertain. Using a new method based on standard thermodynamic cycles and a linear approximation of polar and non-polar free energy contributions from MD averages, binding affinities of the different ligands in each binding site have been correlated with experimental dissociation constants. The study has provided insights into structure-activity relationships for use in the design of modified inhibitors of DHFR.