Distal pocket residues affect picosecond ligand recombination in myoglobin. An experimental and molecular dynamics...

Time courses for intramolecular NO and O2 recombination to native and three position 29 mutants of sperm whale myoglobins were measured after laser photolysis on picosecond and nanosecond time scales. The rates for the first phase of NO recombination were 1.8, 2.5, 29, and > or = 100 ns-1 for Ala29, Val29, Leu29 (native), and Phe29 myoglobin, respectively, at room temperature. This order is not correlated with the overall association rate constants for NO binding which were all in the range 20-50 x 10(6) M-1 s-1 and is the opposite of that observed for the rate constants for the overall thermal dissociation of NO which were 5.0, 2.8, 0.98, and 0.21 x 10(-4) s-1 for Ala29, Val29, Leu29 (native), and Phe29 myoglobin, respectively, at 20 degrees C. This inverse correlation suggests that photo- and thermally dissociated ligand molecules experience similar kinetic and equilibrium barriers to rebinding. The larger side chains of Leu29 and Phe29 inhibit rapid movement of the ligand away from the iron atom facilitating geminate recombination. The smaller side chains of Val29 and Ala29 increase the space available to the ligand, decreasing the rate of geminate recombination and enhancing complete dissociation. Diffusion of NO in the distal pocket of myoglobin was simulated using a variant of the molecular dynamics program CHARMM that includes the locally enhanced sampling protocol (Elber, R., and Karplus, M. (1991) J. Am. Chem. Soc. 112, 9161-9175; Roitberg, A., and Elber, R. (1991) J. Chem. Phys. 95, 9277-9287) and the x-ray structures of Carver et al. (Carver, T. E., Brantley, R. E., Jr., Singleton, E. W., Arduini, R. M., Quillin, M. L., Phillips, G. N., Jr., and Olson, J. S. (1992) J. Biol. Chem. 267, 14443-14450). Both accelerated (5,000 K) and room temperature ligands were used, and comparisons were made between simulations with a complete hydration shell surrounding the protein and those with only eight water molecules near the distal histidine. Photodissociated ligands initially move away from the heme plane, past Leu29, and toward Leu32, Phe33, Ile107, and Ile111. These theoretical results confirm that a complete description of picosecond ligand recombination must include the dynamics of ligand movement in the distal portion of the heme pocket.