The mechanistic basis of substrate-recognition by archtypal metalloenzyme cytochrome P450cam lacks clarity due to reports of contrasting crystallographic poses in its substrate-free form. Here we employ atomic resolution computer simulation to quantitatively dissect the conformational diversities of cytochrome P450cam and reconcile them with its substrate-binding mechanism. By initiating unbiased molecular dynamics simulation trajectories from multiple crystal structures of substrate-free P450cam, we first explore the statistical probabilities of each crystallographic pose of making dynamically transition to the other crystallographic pose and vice versa. Our approach allows us to build a comprehensive conformational ensemble of substrate-free P450cam using an aggregated 54 µs simulation trajectories. A Markov state model constructed using the conformational ensemble identifies multiple key macrostates in a broad and heterogenous free energy landscape and in particular, reveals a dynamic equilibrium between a pair of conformations in which the substrate-recognition sites are closed and open respectively. However, the prediction of a significantly high population of closed conformation, coupled with higher flux and faster rate from open → closed transition than its reverse process, dictates that the net conformational equilibrium would be swayed in favour of closed conformation. The investigation infers that, while a potential substrate of P450cam would in principle sample the diverse array of conformations of substrate-free protein, it would mostly encounter a closed conformation and hence would follow an induced-fit based recognition process. Taken together, this work establishes how a quantitative elucidation of conformational heterogeneity in protein would provide crucial insights in the mechanism of potential substrate-recognition process.
IntroductionCytochrome P450s are key hemoprotein monooxygenase family which catalyzes a diverse range of biochemical processes, including metabolism of almost all drugs, lipid and steroid biosynthesis and decomposition of pollutants. 1 They also play a pivotal role in the detoxification of xenobiotics, cellular metabolism, homeostasis and therefore, are of major clinical significance. 2 These superfamily are especially notable for catalysing the hydroxylation process of unactivated hydrophobic molecules, 3 making them attractive models for investigation of determinants of biological molecular recognition process. 4-6 Among various Cytochrome P450 family, CYP101A1, from the Pseudomonas Putida has served as the archtypal model for thermodynamic and kinetic investigation of substrate recognition in its heme centre. 7,8 Popularly known as P450cam, this specific subfamily of cytochrome P450 catalyzes the hydroxylation process of D-camphor to 5-exohydroxycamphor. 9,10 P450cam is known for its high substrate specificity in comparison to the specificities of many members of the P450 superfamily and this has been attributed to its relatively small, rigid active site that showed few ...