Disentangling ligand migration and heme pocket relaxation in cytochrome P450cam

Ligand diffusion in proteins via enhanced sampling in molecular dynamics

Computational simulations in biophysics describe the dynamics and functions of biological macromolecules at the atomic level. Among motions particularly important for life are the transport processes in heterogeneous media. The process of ligand diffusion inside proteins is an example of a complex rare event that can be modeled using molecular dynamics simulations. The study of physical interactions between a ligand and its biological target is of paramount importance for the design of novel drugs and enzymes. Unfortunately, the process of ligand diffusion is difficult to study experimentally. The need for identifying the ligand egress pathways and understanding how ligands migrate through protein tunnels has spurred the development of several methodological approaches to this problem. The complex topology of protein channels and the transient nature of the ligand passage pose difficulties in the modeling of the ligand entry/escape pathways by canonical molecular dynamics simulations. In this review, we report a methodology involving a reconstruction of the ligand diffusion reaction coordinates and the free-energy profiles along these reaction coordinates using enhanced sampling of conformational space. We illustrate the above methods on several ligand-protein systems, including cytochromes and G-protein-coupled receptors. The methods are general and may be adopted to other transport processes in living matter.

Enzyme transient state kinetics in crystal and solution from the perspective of a time-resolved crystallographer

With recent technological advances at synchrotrons [Graber et al., J. Synchrotron Radiat. 18, 658-670 (2011)], it is feasible to rapidly collect time-resolved crystallographic data at multiple temperature settings [Schmidt et al., Acta Crystallogr. D 69, 2534-2542 (2013)], from which barriers of activation can be extracted. With the advent of fourth generation X-ray sources, new opportunities emerge to investigate structure and dynamics of biological macromolecules in real time [M. Schmidt, Adv. Condens. Matter Phys. 2013, 1-10] in crystals and potentially from single molecules in random orientation in solution [Poon et al., Adv. Condens. Matter Phys. 2013, 750371]. Kinetic data from time-resolved experiments on short time-scales must be interpreted in terms of chemical kinetics [Steinfeld et al., Chemical Kinetics and Dynamics, 2nd ed. (Prentience Hall, 1985)] and tied to existing time-resolved experiments on longer time-scales [Schmidt et al., Acta Crystallogr. D 69, 2534-2542 (2013); Jung et al., Nat. Chem. 5, 212-220 (2013)]. With this article, we will review and outline steps that are required to routinely determine the energetics of reactions in biomolecules in crystal and solution with newest X-ray sources. In eight sections, we aim to describe concepts and experimental details that may help to inspire new approaches to collect and interpret these data.

Cryoradiolysis and cryospectroscopy for studies of heme-oxygen intermediates in cytochromes p450

Cryogenic radiolytic reduction is one of the most straightforward and convenient methods of generation and stabilization of reactive iron-oxygen intermediates for mechanistic studies in chemistry and biochemistry. The method is based on one-electron reduction of the precursor complex in frozen solution via exposure to the ionizing radiation at cryogenic temperatures. Such approach allows for accumulation of the fleeting reactive complexes which otherwise could not be generated at sufficient amount for structural and mechanistic studies. Application of this method allowed for characterizing of peroxo-ferric and hydroperoxo-ferric intermediates, which are common for the oxygen activation mechanism in cytochromes P450, heme oxygenases, and nitric oxide synthases, as well as for the peroxide metabolism by peroxidases and catalases.

Protein dynamics in cytochrome P450 molecular recognition and substrate specificity using 2D IR vibrational echo spectroscopy

Cytochrome (cyt) P450s hydroxylate a variety of substrates that can differ widely in their chemical structure. The importance of these enzymes in drug metabolism and other biological processes has motivated the study of the factors that enable their activity on diverse classes of molecules. Protein dynamics have been implicated in cyt P450 substrate specificity. Here, 2D IR vibrational echo spectroscopy is employed to measure the dynamics of cyt P450(cam) from Pseudomonas putida on fast time scales using CO bound at the active site as a vibrational probe. The substrate-free enzyme and the enzyme bound to both its natural substrate, camphor, and a series of related substrates are investigated to explicate the role of dynamics in molecular recognition in cyt P450(cam) and to delineate how the motions may contribute to hydroxylation specificity. In substrate-free cyt P450(cam), three conformational states are populated, and the structural fluctuations within a conformational state are relatively slow. Substrate binding selectively stabilizes one conformational state, and the dynamics become faster. Correlations in the observed dynamics with the specificity of hydroxylation of the substrates, the binding affinity, and the substrates' molecular volume suggest that motions on the hundreds of picosecond time scale contribute to the variation in activity of cyt P450(cam) toward different substrates.

Interactions of cytochrome P450s with their ligands

Cytochrome P450s (CYPs) are heme-containing monooxygenases that contribute to an enormous range of enzymatic function including biosynthetic and detoxification roles. This review summarizes recent studies concerning interactions of CYPs with ligands including substrates, inhibitors, and diatomic heme-ligating molecules. These studies highlight the complexity in the relationship between the heme spin state and active site occupancy, the roles of water in directing protein-ligand and ligand-heme interactions, and the details of interactions between heme and gaseous diatomic CYP ligands. Both kinetic and thermodynamic aspects of ligand binding are considered.

Spectroscopic features of cytochrome P450 reaction intermediates

Cytochromes P450 constitute a broad class of heme monooxygenase enzymes with more than 11,500 isozymes which have been identified in organisms from all biological kingdoms [1]. These enzymes are responsible for catalyzing dozens chemical oxidative transformations such as hydroxylation, epoxidation, N-demethylation, etc., with very broad range of substrates [2,3]. Historically these enzymes received their name from 'pigment 450' due to the unusual position of the Soret band in UV-vis absorption spectra of the reduced CO-saturated state [4,5]. Despite detailed biochemical characterization of many isozymes, as well as later discoveries of other 'P450-like heme enzymes' such as nitric oxide synthase and chloroperoxidase, the phenomenological term 'cytochrome P450' is still commonly used as indicating an essential spectroscopic feature of the functionally active protein which is now known to be due to the presence of a thiolate ligand to the heme iron [6]. Heme proteins with an imidazole ligand such as myoglobin and hemoglobin as well as an inactive form of P450 are characterized by Soret maxima at 420nm [7]. This historical perspective highlights the importance of spectroscopic methods for biochemical studies in general, and especially for heme enzymes, where the presence of the heme iron and porphyrin macrocycle provides rich variety of specific spectroscopic markers available for monitoring chemical transformations and transitions between active intermediates of catalytic cycle.

Reaction of the Co(II)-substrate radical pair catalytic intermediate in coenzyme B12-dependent ethanolamine ammonia-lyase in frozen aqueous solution from 190 to 217 K

The decay kinetics of the aminoethanol-generated Co(II)-substrate radical pair catalytic intermediate in ethanolamine ammonia-lyase from Salmonella typhimurium have been measured on timescales of <10(5) s in frozen aqueous solution from 190 to 217 K. X-band continuous-wave electron paramagnetic resonance (EPR) spectroscopy of the disordered samples has been used to continuously monitor the full radical pair EPR spectrum during progress of the decay after temperature step reaction initiation. The decay to a diamagnetic state is complete and no paramagnetic intermediate states are detected. The decay exhibits three kinetic regimes in the measured temperature range, as follows. i), Low temperature range, 190 < or = T < or = 207 K: the decay is biexponential with constant fast (0.57 +/- 0.04) and slow (0.43 +/- 0.04) phase amplitudes. ii), Transition temperature range, 207 < T < 214 K: the amplitude of the slow phase decreases to zero with a compensatory rise in the fast phase amplitude, with increasing temperature. iii), High temperature range, T > or = 214 K: the decay is monoexponential. The observed first-order rate constants for the monoexponential (k(obs,m)) and the fast phase of the biexponential decay (k(obs,f)) adhere to the same linear relation on an lnk versus T(-1) (Arrhenius) plot. Thus, k(obs,m) and k(obs,f) correspond to the same apparent Arrhenius prefactor and activation energy (logA(app,f) (s(-1)) = 13.0, E(a,app,f) = 15.0 kcal/mol), and therefore, a common decay mechanism. We propose that k(obs,m) and k(obs,f) represent the native, forward reaction of the substrate through the radical rearrangement step. The slow phase rate constant (k(obs,s)) for 190 < or = T < or = 207 K obeys a different linear Arrhenius relation (logA(app,s) (s(-1)) = 13.9, E(a,app,s) = 16.6 kcal/mol). In the transition temperature range, k(obs,s) displays a super-Arrhenius increase with increasing temperature. The change in E(a,app,s) with temperature and the narrow range over which it occurs suggest an origin in a liquid/glass or dynamical transition. A discontinuity in the activation barrier for the chemical reaction is not expected in the transition temperature range. Therefore, the transition arises from a change in the properties of the protein. We propose that a protein dynamical contribution to the reaction, which is present above the transition temperature, is lost below the transition temperature, owing to an increase in the activation energy barrier for protein motions that are coupled to the reaction. For both the fast and slow phases of the low temperature decay, the dynamical transition in protein motions that are obligatorily coupled to the reaction of the Co(II)-substrate radical pair lies below 190 K.

CO migration pathways in cytochrome P450cam studied by molecular dynamics simulations

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.