Previous laser flash photolysis investigations between 100 and 300 K have shown that the kinetics of CO rebinding with cytochrome P450 cam (camphor) consist of up to four different processes revealing a complex internal dynamics after ligand dissociation. In the present work, molecular dynamics simulations were undertaken on the ternary complex P450 cam (cam)(CO) to explore the CO migration pathways, monitor the internal cavities of the protein, and localize the CO docking sites. One trajectory of 1 nsec with the protein in a water box and 36 trajectories of 1 nsec in the vacuum were calculated. In each trajectory, the protein contained only one CO ligand on which no constraints were applied. The simulations were performed at 200, 300, and 320 K. The results indicate the presence of seven CO docking sites, mainly hydrophobic, located in the same moiety of the protein. Two of them coincide with xenon binding sites identified by crystallography. The protein matrix exhibits eight persistent internal cavities, four of which corresponding to the ligand docking sites. In addition, it was observed that water molecules entering the protein were mainly attracted into the polar pockets, far away from the CO docking sites. Finally, the identified CO migration pathways provide a consistent interpretation of the experimental rebinding kinetics. (Srajer et al. 2001;Bourgeois et al. 2003;Schotte et al. 2003;Hummer et al. 2004), cryocrystallography (Brunori et al. 2000Chu et al. 2000;Ostermann et al. 2000;Schlichting 2000), spectroscopy (Engler et al. 2000;Lamb et al. 2002;Nienhaus et al. 2002;Kriegl et al. 2003), kinetic competition experiments with xenon (Scott and Gibson 1997;Scott et al. 2001; Nienhaus et al. 2003a,b;Tetreau et al. 2004, and molecular dynamics (MD) simulations (Elber and Karplus 1990;Carlson et al. 1996;Meller and Elber 1998;Amara et al. 2001;Mouawad et al. 2005). Article published online ahead of print. Article and publication date are at