Structural and dynamical basis of broad substrate specificity, catalytic mechanism, and inhibition of cytochrome P450 3A4

Cytochrome P450 (CYP) 3A4 is responsible for the oxidative degradation of more than 50% of clinically used drugs. By means of molecular dynamics simulations with the newly developed force field parameters for the heme-thiolate group and its dioxygen adduct, we examine the differences in structural and dynamic properties between CYP3A4 in the resting form and its complexes with the substrate progesterone and the inhibitor metyrapone. The results indicate that the broad substrate specificity of CYP3A4 stems from the malleability of a loop (residues 211-218) that resides in the vicinity of the channel connecting the active site and bulk solvent. However, the high-amplitude motion of the flexible loop is found to be damped out upon binding of the inhibitor or the substrate in the active site. In the resting form of CYP3A4, a structural water molecule is bound to the sixth coordination position of the heme iron, stabilizing the octahedral coordination geometry. In addition to the direct coordination of metyrapone to the heme iron, the hydrogen bond interaction between the inhibitor carbonyl group and the side chain of Ser119 also contributes significantly to stabilizing the CYP3A4-metyrapone complex. On the other hand, progesterone is stabilized in the active site by the formation of two hydrogen bonds with Ser119 and Arg106, as well as by the van der Waals interactions with the heme and hydrophobic residues. The structural and dynamic features of the CYP3A4-progesterone complex indicate that the oxidative degradation of progesterone occurs through hydroxylation at the C16 position by the reactive oxygen coordinated to the heme iron.     

Phenotyping of CYP450 in human liver microsomes using the cocktail approach

The cocktail approach is an advantageous strategy used to monitor the activities of several cytochromes P450 (CYPs) in a single test to increase the throughput of in vitro phenotyping studies. In this study, a cocktail mixture was developed with eight CYP-specific probe substrates to simultaneously evaluate the activity of the most important CYPs, namely, CYP1A2, CYP2A6, CYP2B6, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and the CYP3A subfamily. After cocktail incubation in the presence of human liver microsomes (HLMs), the eight selected substrates and their specific metabolites were analyzed by ultra-high-pressure liquid chromatography and electrospray ionization quadrupole time-of-flight mass spectrometry. Qualitative and quantitative data were simultaneously acquired to produce an overview of the extended phase I biotransformation routes for each probe substrate in the HLMs and to generate phenotypic profiles of various HLMs. A comparison of the cocktail strategy with an individual substrate assay for each CYP produced similar results. Moreover, the cocktail was tested on HLMs with different allelic variants and/or in the presence of selective inhibitors. The results were in agreement with the genetic polymorphisms of the CYPs and the expected effect of the alterations. All of these experiments confirmed the reliability of this cocktail assay for phenotyping of the microsomal CYPs.     

Patterns in Protein Flexibility: A Comparison of NMR “Ensembles”, MD Trajectories,and Crystallographic B-Factors

Proteins are molecular machines requiring flexibility to function. Crystallographic B-factors and Molecular Dynamics (MD) simulations both provide insights into protein flexibility on an atomic scale. Nuclear Magnetic Resonance (NMR) lacks a universally accepted analog of the B-factor. However, a lack of convergence in atomic coordinates in an NMR-based structure calculation also suggests atomic mobility. This paper describes a pattern in the coordinate uncertainties of backbone heavy atoms in NMR-derived structural “ensembles” first noted in the development of FindCore2 (previously called Expanded FindCore: DA Snyder, J Grullon, YJ Huang, R Tejero, GT Montelione, Proteins: Structure, Function, and Bioinformatics 82 (S2), 219–230) and demonstrates that this pattern exists in coordinate variances across MD trajectories but not in crystallographic B-factors. This either suggests that MD trajectories and NMR “ensembles” capture motional behavior of peptide bond units not captured by B-factors or indicates a deficiency common to force fields used in both NMR and MD calculations.     

Pyrene.pyrene complexes at the active site of cytochrome P450 3A4: evidence for a multiple substrate binding site

Cytochrome P450 monooxygenases (CYPs) metabolize nearly all drugs and toxins. Recently, it has become clear that CYPs exhibit both homotropic and heterotropic allosteric kinetics for many substrates. However, the mechanism of cooperative kinetics has not been established for any specific human CYP/substrate combination. Suggested mechanisms include binding of multiple substrates within distinct, static, subsites of a single large active site or binding of multiple substrates within a single fluid active site. CYP3A4 hydroxylates pyrene with positive cooperativity. Therefore, experiments were designed to exploit the fluorescence properties of pyrene, which diagnostically distinguish between pyrene.pyrene complexes versus spatially separated pyrene substrates. Pyrene complexes (excimers) yield an emission spectrum clearly distinct from pyrene monomers. In lipid-free aqueous/glycerol solutions of CYP3A4, addition of pyrene affords a concentration-dependent low-spin to high-spin conversion of the CYP3A4 heme prosthetic group, indicating occupancy of the active site by pyrene. Under the same conditions, in the presence of CYP3A4 but not other heme proteins, the excimer/monomer ratio (E/M) of pyrene was decreased in emission spectra, compared to pyrene alone. However, excitation spectra indicate a CYP3A4-dependent increase in the wavelength shift for the excimer excitation spectrum versus the monomer excitation spectrum, as well as changes in the excimer excitation peak shape and vibronic structure. These changes are reversed by the CYP3A4 substrate testosterone. Together, the results demonstrate that pyrene.pyrene ground-state complexes occupy the CYP3A4 active site, and they provide the first spectroscopic evidence for substrate complexes within a single fluid active site. Functional implications include the possibility that turnover rate, regioselectivity, and stereoselectivity of the reaction are determined by the substrate.substrate complex rather than individual substrates.     

Characterization of the CYP3A4 active site by homology modeling

Human microsomal cytochrome P450s participate in drug metabolism and detoxification. Among them, CYP3A4 is the most important isoform for drug-drug interactions. To gain a better understanding of the active site, a homology model of CYP3A4 was constructed based on the crystallographic coordinates of mammalian CYP2C5. The putative active site is much larger than that of CYP2C5 and is divided into three parts (i.e. a proximal and two distal sites from the heme). Most residues reported to be important for ligand-binding are located in the active site of the model. Moreover, some inhibitors (paclitaxel etc.) docked into the model have complementary shapes to the pocket. Pharmacophore docking of 14 substrates was also performed using Ph4Dock of MOE. Calculated interaction energies showed a moderate correlation with the logarithm of apparent K(m) values. These results suggest that this model is reliable enough to be used in the design of compounds for removing undesirable CYP3A4 inhibition.     

Fourier transform infrared spectroscopy as a tool to study structural properties of cytochromes P450 (CYPs)

Cytochrome P450 proteins (CYPs) are a big class of heme proteins which are involved in various metabolic processes of living organisms. CYPs are the terminal catalytically active components of monooxygenase systems where the substrate binds and is hydroxylated. In order to be functionally competent, the protein structures of CYPs possess specific properties that must be explored in order to understand structure–function relationships and mechanistic aspects. Fourier transform infrared spectroscopy (FTIR) is one tool that is used to study these structural properties. The application of FTIR spectroscopy to the secondary structures of CYP proteins, protein unfolding, protein–protein interactions and the structure and dynamics of the CYP heme pocket is reviewed. A comparison with other thiolate heme proteins (nitric oxide synthase and chloroperoxidase) is also included. Figure The protein secondary structure, protein unfolding, redox-partner protein–protein interaction, structural changes induced by the reduction of the heme iron, and the structure and dynamics of the active site of cytochromes P450 (CYP) can be studied using Fourier transform infrared spectroscopy (FTIR). FTIR spectroscopy is a good approach for gaining a deeper insight into structure–function relationships in CYPs. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Molecular Docking of 3-Methylindole-containing Drugs Binding into CYP3A4

Drugs SPD-304(6,7-dimethyl-3-{[methyl-(2-{methyl-[1-(3-trifluoromethyl-phenyl)-1H-indol-3-ylmethyl]- amino}-ethyl)-amino]-methyl}-chromen-4-one) and zafirlukast contain a common structural element of 3-substituted indole moiety which closely relates to a dehydrogenated reaction catalyzed by cytochrome P450s(CYPs). It was reported that the dehydrogenation can produce a reactive electrophilic intermediate which cause toxicities and inactivate CYPs. Drug L-745,870(3-{[4-(4-chlorophenyl)piperazin-1-yl]-methyl}-1H-pyrrolo- 2,3-β-pyridine) might have similar effect since it contains the same structural element. We used molecular docking approach combined with molecular dynamics(MD) simulation to model three-dimensional(3D) complex structures of SPD-304, zafirlukast and L-745,870 into CYP3A4, respectively. The results show that these three drugs can stably bind into the active site and the 3-methylene carbons of the drugs keep a reasonable reactive distance from the heme iron. The complex structure of SPD-304-CYP3A4 is in agreement with experimental data. For zafirlukast, the calculation results indicate that 3-methylene carbon might be the dehydrogenation reaction site. Docking model of L-745,870-CYP3A4 shows a potential possibility of L-745,870 dehydrogenated by CYP3A4 at 3-methylene carbon which is in agreement with experiment in vivo. In addition, residues in the phenylalanine cluster as well as S119 and R212 play a critical role in the ligands binding based on our calculations. The docking models could provide some clues to understand the metabolic mechanism of the drugs by CYP3A4.     

Selenolate complexes of CYP101 and the heme-bound hHO-1/H25A proximal cavity mutant

Thiolate and selenolate complexes of CYP101 (P450cam) and the H25A proximal cavity mutant of heme-bound human heme oxygenase-1 (hHO-1) have been examined by UV-vis spectroscopy. Both thiolate and selenolate ligands bound to the heme distal side in CYP101 and gave rise to characteristic hyperporphyrin spectra. Thiolate ligands also bound to the proximal side of the heme in the cavity created by the H25A mutation in hHO-1, giving a Soret absorption similar to that of the H25C hHO-1 mutant. Selenolate ligands also bound to this cavity mutant under anaerobic conditions but reduced the heme iron to the ferrous state, as shown by the formation of a ferrous CO complex. Under aerobic conditions, the selenolate ligand but not the thiolate ligand was rapidly oxidized. These results indicate that selenocysteine-coordinated heme proteins will not be stable species in the absence of a redox potential stabilizing effect.     

The crystal structures of 4-methoxybenzoate bound CYP199A2 and CYP199A4: structural changes on substrate binding and the identification of an anion binding site

The crystal structures of the 4-methoxybenzoate bound forms of cytochrome P450 enzymes CYP199A2 and CYP199A4 from the Rhodopseudomonas palustris strains CGA009 and HaA2 have been solved. The structures of these two enzymes, which share 86% sequence identity, are very similar though some differences are found on the proximal surface. In these structures the enzymes have a closed conformation, in contrast to the substrate-free form of CYP199A2 where an obvious substrate access channel is observed. The switch from an open to a closed conformation arises from pronounced residue side-chain movements and alterations of ion pair and hydrogen bonding interactions at the entrance of the access channel. A chloride ion bound just inside the protein surface caps the entrance to the active site and protects the substrate and the heme from the external solvent. In both structures the substrate is held in place via hydrophobic and hydrogen bond interactions. The methoxy group is located over the heme iron, accounting for the high activity and selectivity of these enzymes for oxidative demethylation of the substrate. Mutagenesis studies on CYP199A4 highlight the involvement of hydrophobic (Phe185) and hydrophilic (Arg92, Ser95 and Arg243) amino acid residues in the binding of para-substituted benzoates by these enzymes.     

Insights on cytochrome p450 enzymes and inhibitors obtained through QSAR studies

The cytochrome P450 (CYP) superfamily of heme enzymes play an important role in the metabolism of a large number of endogenous and exogenous compounds, including most of the drugs currently on the market. Inhibitors of CYP enzymes have important roles in the treatment of several disease conditions such as numerous cancers and fungal infections in addition to their critical role in drug-drug interactions. Structure activity relationships (SAR), and three-dimensional quantitative structure activity relationships (3D-QSAR) represent important tools in understanding the interactions of the inhibitors with the active sites of the CYP enzymes. A comprehensive account of the QSAR studies on the major human CYPs 1A1, 1A2, 1B1, 2A6, 2B6, 2C9, 2C19, 2D6, 2E1, 3A4 and a few other CYPs are detailed in this review which will provide us with an insight into the individual/common characteristics of the active sites of these enzymes and the enzyme-inhibitor interactions.     

Coexpression of CPR from Various Origins Enhances Biotransformation Activity of Human CYPs in S. pombe

Cytochrome P450 enzymes (CYPs or P450s) are the most important enzymes involved in the phase I metabolism of drugs (and other xenobiotics) in humans, and the corresponding drug metabolites are needed as reference substances for their structural confirmation and for pharmacological or toxicological characterization. We have previously shown that biotechnological synthesis of such metabolites is feasible by whole-cell biotransformation with human CYPs recombinantly expressed in the fission yeast Schizosaccharomyces pombe. It was the aim of this study to compare the activity of seven human microsomal CYPs (CYP2C9, CYP2D6, CYP3A4, CYP3A5, CYP3A7, CYP17, and CYP21) upon coexpression with NADPH-cytochrome P450 oxidoreductases (CPRs) from various origins, namely, human CPR (hCPR) and its homologues from fission yeast (ccr1) and the bishop’s weed Ammi majus (AmCPR), respectively. For this purpose, 28 recombinant strains were needed, with five of them having been constructed previously and 23 strains being newly constructed. Bioconversion experiments showed that coexpression of a CPR does not only influence the reaction rate but, in some cases, also exerts an influence on the metabolite pattern. For CYP3A enzymes, coexpression of hCPR yielded the best results, while for another two, hCPR was equally helpful as ccr1 (both CYP17 and CYP21) or AmCPR (CYP17 only), respectively. Interestingly, CYP2D6 displayed its highest activity when coexpressed with ccr1 and CYP2C9 with AmCPR. These results corroborate the view of CPR as a well-suited bio-brick in synthetic biology for the construction of artificial enzyme complexes.     

How Mutations Affecting the Ligand-receptor Interactions:a Combined MD and QM/MM Calculation on CYP2E1 and Its Two Mutants

Cytochrome P450(CYP) 2E1 is a dual function monoxygenase with a crucial role in the metabolism of 6% of drugs on the market at present. The enzyme is of tremendous interest for its association with alcohol consumption, diabetes, obesity and fasting. Despite the abundant experimental mutagenesis data, the molecular origin and the structural motifs for the enzymatic activity deficiencies have not been rationalized at the atomic level. In this regard, we have investigated the effects of mutation on the structural and energetic characteristics upon single point mutations in CYP2E1, N219D and S366C. The molecular dynamics(MD) simulation combined with quantum mechanics/molecular mechanics(QM/MM) and noncovalent interaction(NCI) analysis was carried out on CYP2E1 and its two mutants. The results highlight the critical role of Phe207, which is responsible for both structural flexibility and energetic variation, shortening the gap between the theory and the experimentally observed results of enzymatic activity decrease. The underlying molecular mechanism of the enzymatic activity deficiencies for mutants may be attributed to the changes of spatial position of Phe207 in the two mutants. This work provides particular explanations to how mutations affect ligand-receptor interactions based on combined MD and QM/MM calculations. Furthermore, the mutational effects on the activity of CYP2E1 obtained in the present study are beneficial to both the experimental and the computational works of CYPs and may allow researchers to achieve desirable changes in enzymatic activity.     

Controlling the Substrate Specificity of an Enzyme through Structural Flexibility by Varying the Salt-Bridge Density

Many enzymes, particularly in one single family, with highly conserved structures and folds exhibit rather distinct substrate specificities. The underlying mechanism remains elusive, the resolution of which is of great importance for biochemistry, biophysics, and bioengineering. Here, we performed a neutron scattering experiment and molecular dynamics (MD) simulations on two structurally similar CYP450 proteins; CYP101 primarily catalyzes one type of ligands, then CYP2C9 can catalyze a large range of substrates. We demonstrated that it is the high density of salt bridges in CYP101 that reduces its structural flexibility, which controls the ligand access channel and the fluctuation of the catalytic pocket, thus restricting its selection on substrates. Moreover, we performed MD simulations on 146 different kinds of CYP450 proteins, spanning distinct biological categories including Fungi, Archaea, Bacteria, Protista, Animalia, and Plantae, and found the above mechanism generally valid. We demonstrated that, by fine changes of chemistry (salt-bridge density), the CYP450 superfamily can vary the structural flexibility of its member proteins among different biological categories, and thus differentiate their substrate specificities to meet the specific biological needs. As this mechanism is well-controllable and easy to be implemented, we expect it to be generally applicable in future enzymatic engineering to develop proteins of desired substrate specificities.     

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

Cytochrome P450 2C9 (CYP2C9) is a membrane-anchored human microsomal protein involved in the drug metabolism in liver. CYP2C9 consists of an N-terminal transmembrane anchor and a catalytic cytoplasmic domain. While the structure of the catalytic domain is well-known from X-ray experiments, the complete structure and its incorporation into the membrane remains unsolved. We constructed an atomistic model of complete CYP2C9 in a dioleoylphosphatidylcholine membrane and evolved it by molecular dynamics simulations in explicit water on a 100+ ns time-scale. The model agrees well with known experimental data about membrane positioning of cytochromes P450. The entry to the substrate access channel is proposed to be facing the membrane interior while the exit of the product egress channel is situated above the interface pointing toward the water phase. The positions of openings of the substrate access and product egress channels correspond to free energy minima of CYP2C9 substrate ibuprofen and its metabolite in the membrane, respectively.     

Using a homology model of cytochrome P450 2D6 to predict substrate site of metabolism

CYP2D6 is an important enzyme that is involved in first pass metabolism and is responsible for metabolizing ~25% of currently marketed drugs. A homology model of CYP2D6 was built using X-ray structures of ligand-bound CYP2C5 complexes as templates. This homology model was used in docking studies to rationalize and predict the site of metabolism of known CYP2D6 substrates. While the homology model was generally found to be in good agreement with the recently solved apo (ligand-free) X-ray structure of CYP2D6, significant differences between the structures were observed in the B′ and F–G helical region. These structural differences are similar to those observed between ligand-free and ligand-bound structures of other CYPs and suggest that these conformational changes result from induced-fit adaptations upon ligand binding. By docking to the homology model using Glide, it was possible to identify the correct site of metabolism for a set of 16 CYP2D6 substrates 85% of the time when the 5 top scoring poses were examined. On the other hand, docking to the apo CYP2D6 X-ray structure led to a loss in accuracy in predicting the sites of metabolism for many of the CYP2D6 substrates considered in this study. These results demonstrate the importance of describing substrate-induced conformational changes that occur upon binding. The best results were obtained using Glide SP with van der Waals scaling set to 0.8 for both the receptor and ligand atoms. A discussion of putative binding modes that explain the distribution of metabolic sites for substrates, as well as a relationship between the number of metabolic sites and substrate size, are also presented. In addition, analysis of these binding modes enabled us to rationalize the typical hydroxylation and O-demethylation reactions catalyzed by CYP2D6 as well as the less common N-dealkylation.     

Allosteric effects on substrate dissociation from cytochrome P450 3A4 in nanodiscs observed by ensemble and single-molecule fluorescence spectroscopy

Cytochrome P450 (CYP) 3A4 is a major human drug-metabolizing enzyme and displays pharmacologically relevant allosteric kinetics caused by multiple substrate and/or effector binding. Here, in the first single-molecule (SM) fluorescence studies of CYPs, we use total internal reflection fluorescence microscopy to measure residence times of the fluorescent dye Nile Red in CYP3A4 incorporated in surface-immobilized lipid Nanodiscs, with and without the effector alpha-naphthoflavone. We find direct evidence that CYP3A4 effectors can decrease substrate off-rates, providing a possible mechanism for effector-mediated enhancement of substrate metabolism. These interesting results highlight the potential of SM methods in studies of CYP allosteric mechanisms.     

Investigation of the Substrate Range of CYP199A4: Modification of the Partition between Hydroxylation and Desaturation Activities by Substrate and Protein Engineering

The cytochrome P450 enzyme CYP199A4, from Rhodopseudomonas palustris HaA2, can efficiently demethylate 4‐methoxybenzoic acid. It is also capable of oxidising a range of other related substrates. By investigating substrates with different substituents and ring systems we have been able to show that the carboxylate group and the nature of the ring system and the substituent are all important for optimal substrate binding and activity. The structures of the veratric acid, 2‐naphthoic acid and indole‐6‐carboxylic acid substrate‐bound CYP199A4 complexes reveal the substrate binding modes and the side‐chain conformational changes of the active site residues to accommodate these larger substrates. They also provide a rationale for the selectivity of product oxidation. The oxidation of alkyl substituted benzoic acids by CYP199A4 is more complex, with desaturation reactions competing with hydroxylation activity. The structure of 4‐ethylbenzoic acid‐bound CYP199A4 revealed that the substrate is held in a similar position to 4‐methoxybenzoic acid, and that the C^β C? H bonds of the ethyl group are closer to the heme iron than those of the C^α (3.5 vs. 4.8 Å). This observation, when coupled to the relative energies of the reaction intermediates, indicates that the positioning of the alkyl group relative to the heme iron may be critical in determining the amount of desaturation that is observed. By mutating a single residue in the active site of CYP199A4 (Phe185) we were able to convert the enzyme into a 4‐ethylbenzoic acid desaturase.     

Contact angle hysteresis at the nanoscale: a molecular dynamics simulation study

In this paper, the contact angle hysteresis (CAH) of nanodroplets on both rigid and flexible substrates with different wettabilities was studied using molecular dynamics (MD) simulations. The critical shear stress (CSS) that determines the motion of the contact line (CL) was investigated. A theoretical correlation between CAH and CSS was proposed. Both CAH and CSS reflect the energy dissipation at the CL of the droplet in response to the exerted force. MD results of CAH are qualitatively consistent with the theoretical model. Simulation results also show that, for the same liquid–solid interactions, CAH on the flexible substrate is larger than that on the rigid substrate. These findings aim to enhance our understanding of the mechanism of the CAH at the nanoscale.