Membrane position of ibuprofen agrees with suggested access path entrance to cytochrome P450 2C9 active site

Nanodiscs in the studies of membrane-bound cytochrome P450 enzymes

Cytochromes P450 from eukaryotes and their native redox partners cytochrome P450 reductases both belong to the class of monotopic membrane proteins containing one transmembrane anchor. Incorporation into the lipid bilayer significantly affects their equilibrium and kinetic properties and plays an important role in their interactions. We describe here the detailed protocols developed in our group for the functional self-assembly of mammalian cytochromes P450 and cytochrome P450 reductases into Nanodiscs with controlled lipid composition. The resulting preparations are fully functional, homogeneous in size, composition and oligomerization state of the heme enzyme, and show an improved stability with respect to P420 formation. We provide a brief overview of applications of Nanodisc technology to the biophysical and biochemical mechanistic studies of cytochromes P450 involved in steroidogenesis, and of the most abundant xenobiotic-metabolizing human cytochrome P450 CYP3A4.

Cross-linking mass spectrometry and mutagenesis confirm the functional importance of surface interactions between CYP3A4 and holo/apo cytochrome b(5)

Cytochrome b(5) (cyt b(5)) is one of the key components in the microsomal cytochrome P450 monooxygenase system. Consensus has not been reached about the underlying mechanism of cyt b(5) modulation of CYP catalysis. Both cyt b(5) and apo b(5) are reported to stimulate the activity of several P450 isoforms. In this study, the surface interactions of both holo and apo b(5) with CYP3A4 were investigated and compared for the first time. Chemical cross-linking coupled with mass spectrometric analysis was used to identify the potential electrostatic interactions between the protein surfaces. Subsequently, the models of interaction of holo/apo b(5) with CYP3A4 were built using the identified interacting sites as constraints. Both cyt b(5) and apo b(5) were predicted to bind to the same groove on CYP3A4 with close contacts to the B-B' loop of CYP3A4, a substrate recognition site. Mutagenesis studies further confirmed that the interacting sites on CYP3A4 (Lys96, Lys127, and Lys421) are functionally important. Mutation of these residues reduced or abolished cyt b(5) binding affinity. The critical role of Arg446 on CYP3A4 in binding to cyt b(5) and/or cytochrome P450 reductase was also discovered. The results indicated that electrostatic interactions on the interface of the two proteins are functionally important. The results indicate that apo b(5) can dock with CYP3A4 in a manner analogous to that of holo b(5), so electron transfer from cyt b(5) is not required for its effects.

Conformational adaptation of human cytochrome P450 2B6 and rabbit cytochrome P450 2B4 revealed upon binding multiple amlodipine molecules

Structures of human cytochrome P450 2B6 and rabbit cytochrome P450 2B4 in complex with two molecules of the calcium channel blocker amlodipine have been determined by X-ray crystallography. The presence of two drug molecules suggests clear substrate access channels in each P450. According to a previously established nomenclature, amlodipine molecules were trapped in access pathway 2f in P450 2B6 and in pathway 2a or 2f in P450 2B4. These pathways overlap for part of the length and then diverge as they extend toward the protein surface. A previously described solvent channel was also found in each enzyme. The results indicate that key residues located on the surface and at the entrance of the substrate access channels in each of these P450s may play a crucial role in guiding substrate entry. In addition, the region of P450 2B6 and 2B4 involving helices B', F, F', and G' and part of helix G is substantially more open in the amlodipine complexes than in the corresponding 4-(4-chlorophenyl)imidazole complexes. The increased active site volume observed results from the major retraction of helices F, F', and B' and the β4 sheet region located close to the binding cavity to accommodate amlodipine. These structures demonstrate novel insight into distinct conformational states not observed with previous P450 2B structures and provide clear evidence of the substrate access channels in two drug-metabolizing P450s. In addition, the structures exhibit the versatility that can be exploited via in silico studies with other P450 2B6 ligands as large as raloxifene and itraconazole.

MOLEonline 2.0: interactive web-based analysis of biomacromolecular channels

Biomolecular channels play important roles in many biological systems, e.g. enzymes, ribosomes and ion channels. This article introduces a web-based interactive MOLEonline 2.0 application for the analysis of access/egress paths to interior molecular voids. MOLEonline 2.0 enables platform-independent, easy-to-use and interactive analyses of (bio)macromolecular channels, tunnels and pores. Results are presented in a clear manner, making their interpretation easy. For each channel, MOLEonline displays a 3D graphical representation of the channel, its profile accompanied by a list of lining residues and also its basic physicochemical properties. The users can tune advanced parameters when performing a channel search to direct the search according to their needs. The MOLEonline 2.0 application is freely available via the Internet at http://ncbr.muni.cz/mole or http://mole.upol.cz.

Investigation of indazole unbinding pathways in CYP2E1 by molecular dynamics simulations

Human microsomal cytochrome P450 2E1 (CYP2E1) can oxidize not only low molecular weight xenobiotic compounds such as ethanol, but also many endogenous fatty acids. The crystal structure of CYP2E1 in complex with indazole reveals that the active site is deeply buried into the protein center. Thus, the unbinding pathways and associated unbinding mechanisms remain elusive. In this study, random acceleration molecular dynamics simulations combined with steered molecular dynamics and potential of mean force calculations were performed to identify the possible unbinding pathways in CYP2E1. The results show that channel 2c and 2a are most likely the unbinding channels of CYP2E1. The former channel is located between helices G and I and the B-C loop, and the latter resides between the region formed by the F-G loop, the B-C loop and the β1 sheet. Phe298 and Phe478 act as the gate keeper during indazole unbinding along channel 2c and 2a, respectively. Previous site-directed mutagenesis experiments also supported these findings.

A novel type of allosteric regulation: functional cooperativity in monomeric proteins

Cooperative functional properties and allosteric regulation in cytochromes P450 play an important role in xenobiotic metabolism and define one of the main mechanisms of drug-drug interactions. Recent experimental results suggest that ability to bind simultaneously two or more small organic molecules can be the essential feature of cytochrome P450 fold, and often results in rich and complex pattern of allosteric behavior. Manifestations of non-Michaelis kinetics include homotropic and heterotropic activation and inhibition effects depending on the stoichiometric ratios of substrate and effector, changes in the regio- and stereospecificity of catalytic transformations, and often give rise to the clinically important drug-drug interactions. In addition, functional response of P450 systems is modulated by the presence of specific and non-specific effector molecules, metal ions, membrane incorporation, formation of homo- and hetero-oligomers, and interactions with the protein redox partners. In this article we briefly overview the main factors contributing to the allosteric effects in cytochromes P450 with the main focus on the sources of cooperative behavior in xenobiotic metabolizing monomeric heme enzymes with their conformational flexibility and extremely broad substrate specificity. The novel mechanism of functional cooperativity in P450 enzymes does not require substantial binding cooperativity, rather it implies the presence of one or more binding sites with higher affinity than the single catalytically active site in the vicinity of the heme iron.

Convergence of Free Energy Profile of Coumarin in Lipid Bilayer

Atomistic molecular dynamics (MD) simulations of druglike molecules embedded in lipid bilayers are of considerable interest as models for drug penetration and positioning in biological membranes. Here we analyze partitioning of coumarin in dioleoylphosphatidylcholine (DOPC) bilayer, based on both multiple, unbiased 3 μs MD simulations (total length) and free energy profiles along the bilayer normal calculated by biased MD simulations (∼7 μs in total). The convergences in time of free energy profiles calculated by both umbrella sampling and z-constraint techniques are thoroughly analyzed. Two sets of starting structures are also considered, one from unbiased MD simulation and the other from "pulling" coumarin along the bilayer normal. The structures obtained by pulling simulation contain water defects on the lipid bilayer surface, while those acquired from unbiased simulation have no membrane defects. The free energy profiles converge more rapidly when starting frames from unbiased simulations are used. In addition, z-constraint simulation leads to more rapid convergence than umbrella sampling, due to quicker relaxation of membrane defects. Furthermore, we show that the choice of RESP, PRODRG, or Mulliken charges considerably affects the resulting free energy profile of our model drug along the bilayer normal. We recommend using z-constraint biased MD simulations based on starting geometries acquired from unbiased MD simulations for efficient calculation of convergent free energy profiles of druglike molecules along bilayer normals. The calculation of free energy profile should start with an unbiased simulation, though the polar molecules might need a slow pulling afterward. Results obtained with the recommended simulation protocol agree well with available experimental data for two coumarin derivatives.

Crystal structure of human cytochrome P450 2D6 with prinomastat bound

Human cytochrome P450 2D6 contributes to the metabolism of >15% of drugs used in clinical practice. This study determined the structure of P450 2D6 complexed with a substrate and potent inhibitor, prinomastat, to 2.85 Å resolution by x-ray crystallography. Prinomastat binding is well defined by electron density maps with its pyridyl nitrogen bound to the heme iron. The structure of ligand-bound P450 2D6 differs significantly from the ligand-free structure reported for the P450 2D6 Met-374 variant (Protein Data Bank code 2F9Q). Superposition of the structures reveals significant differences for β sheet 1, helices A, F, F', G", G, and H as well as the helix B-C loop. The structure of the ligand complex exhibits a closed active site cavity that conforms closely to the shape of prinomastat. The closure of the open cavity seen for the 2F9Q structure reflects a change in the direction and pitch of helix F and introduction of a turn at Gly-218, which is followed by a well defined helix F' that was not observed in the 2F9Q structure. These differences reflect considerable structural flexibility that is likely to contribute to the catalytic versatility of P450 2D6, and this new structure provides an alternative model for in silico studies of substrate interactions with P450 2D6.

Structural differences between soluble and membrane bound cytochrome P450s

The superfamily of cytochrome P450s forms a large class of heme monooxygenases with more than 13,000 enzymes represented in organisms from all biological kingdoms. Despite impressive variability in sizes, sequences, location, and function, all cytochrome P450s from various organisms have very similar tertiary structures within the same fold. Here we show that systematic comparison of all available X-ray structures of cytochrome P450s reveals the presence of two distinct structural classes of cytochrome P450s. For all membrane bound enzymes, except the CYP51 family, the beta-domain and the A-propionate heme side chain are shifted towards the proximal side of the heme plane, which may result in an increase of the volume of the substrate binding pocket and an opening of a potential channel for the substrate access and/or product escape directly into the membrane. This structural feature is also observed in several soluble cytochrome P450s, such as CYP108, CYP151, and CYP158A2, which catalyze transformations of bulky substrates. Alternatively, both beta-domains and the A-propionate side chains in the soluble isozymes extend towards the distal site of the heme. This difference between the structures of soluble and membrane bound cytochrome P450s can be rationalized through the presence of several amino acid inserts in the latter class which are involved in direct interactions with the membrane, namely the F'- and G'-helices. Molecular dynamics using the most abundant human cytochrome P450, CYP3A4, incorporated into a model POPC bilayer reveals the facile conservation of a substrate access channel, directed into the membrane between the B-C loop and the beta domain, and the closure of the peripheral substrate access channel directed through the B-C loop. This is in contrast to the case when the same simulation is run in buffer, where no such channel closing occurs. Taken together, these results reveal a key structural difference between membrane bound and soluble cytochrome P450s with important functional implications induced by the lipid bilayer.