Structural features of cytochromes P450 and ligands that affect drug metabolism as revealed by X-ray crystallography and NMR

Comprehensive assessment of Cytochrome P450 reactions: A multiplex approach using real-time ESI-MS

BACKGROUND

The detailed analysis of Cytochrome P450 (CYP) catalyzed reactions is of great interest, since those are of importance for biotechnical applications, drug interaction studies and environmental research. Often cocktail approaches are carried out in order to monitor several CYP activities in a single experiment. Commonly in these approaches product formation is detected and IC50 values are determined.

METHODS

In the present work, the reactions of two different CYP isoforms were monitored using real-time electrospray ionization mass spectrometry. Multiplex experiments using the highly specific CYP2A6 with its corresponding substrate coumarin as well as the highly promiscuous CYP3A4 with testosterone were conducted. Product formation and substrate depletion were simultaneously monitored and compared to the single CYP experiments. The diffusion-controlled rate of reaction and conversion rates that are used as parameters to assess the enzymatic activity were calculated for all measurements conducted.

RESULTS

Differences in conversion rates and the theoretical rate of reaction that were observed for single CYP and multiplex experiments, respectively, reveal the complexity of the underlying mechanisms. Findings of this study imply that there might be distinct deviations between product formation and substrate degradation when mixtures are used.

CONCLUSIONS

Detailed results indicate that for a comprehensive assessment of these enzymatic reactions both product and substrate should be considered.

GENERAL SIGNIFICANCE

The direct hyphenation of enzymatic reactions to mass spectrometry allows for a comprehensive assessment of enzymatic behavior. Due to the benefits of this technique, the entire system which includes substrate, product and intermediates can be investigated. Thus, besides IC50 values further information regarding the enzymatic behavior offers the opportunity for a more detailed insight.

In silico methods for predicting drug-drug interactions with cytochrome P-450s, transporters and beyond

Drug-drug interactions (DDIs) are associated with severe adverse effects that may lead to the patient requiring alternative therapeutics and could ultimately lead to drug withdrawal from the market if they are severe. To prevent the occurrence of DDI in the clinic, experimental systems to evaluate drug interaction have been integrated into the various stages of the drug discovery and development process. A large body of knowledge about DDI has also accumulated through these studies and pharmacovigillence systems. Much of this work to date has focused on the drug metabolizing enzymes such as cytochrome P-450s as well as drug transporters, ion channels and occasionally other proteins. This combined knowledge provides a foundation for a hypothesis-driven in silico approach, using either cheminformatics or physiologically based pharmacokinetics (PK) modeling methods to assess DDI potential. Here we review recent advances in these approaches with emphasis on hypothesis-driven mechanistic models for important protein targets involved in PK-based DDI. Recent efforts with other informatics approaches to detect DDI are highlighted. Besides DDI, we also briefly introduce drug interactions with other substances, such as Traditional Chinese Medicines to illustrate how in silico modeling can be useful in this domain. We also summarize valuable data sources and web-based tools that are available for DDI prediction. We finally explore the challenges we see faced by in silico approaches for predicting DDI and propose future directions to make these computational models more reliable, accurate, and publically accessible.

Size- and time-dependent alteration in metabolic activities of human hepatic cytochrome P450 isozymes by gold nanoparticles via microsomal coincubations

Nano-sized particles are known to interfere with drug-metabolizing cytochrome P450 (CYP) enzymes, which can be anticipated to be a potential source of unintended adverse reactions, but the mechanisms underlying the inhibition are still not well understood. Herein we report a systematic investigation of the impacts of gold nanoparticles (AuNPs) on five major CYP isozymes under in vitro incubations of human liver microsomes (HLMs) with tannic acid (TA)-stabilized AuNPs in the size range of 5 to 100 nm. It is found that smaller AuNPs show more pronounced inhibitory effects on CYP2C9, CYP2C19, CYP2D6, and CYP3A4 in a dose-dependent manner, while 1A2 is the least susceptible to the AuNP inhibition. The size- and dose-dependent CYP-specific inhibition and the nonspecific drug-nanogold binding in the coincubation media can be significantly reduced by increasing the concentration ratio of microsomal proteins to AuNPs, probably via a noncompetitive mode. Remarkably, AuNPs are also found to exhibit a slow time-dependent inactivation of 2D6 and 3A4 in a β-nicotinamide adenine dinucleotide 2'-phosphate reduced tetrasodium salt hydrate (NADPH)-independent manner. During microsomal incubations, UV-vis spectroscopy, dynamic light scattering, and zeta-potential measurements were used to monitor the changes in particle properties under the miscellaneous AuNP/HLM/CYP dispersion system. An improved stability of AuNPs by mixing HLM with the gold nanocolloid reveals that the stabilization via AuNP-HLM interactions may occur on a faster time scale than the salt-induced nanoaggregation by incubation in phosphate buffer. The results suggest that the AuNP induced CYP inhibition can be partially attributed to its adhesion onto the enzymes to alter their structural conformations or onto the HLM membrane therefore impairing the integral membrane proteins. Additionally, AuNPs likely block the substrate pocket on the CYP surface, depending on both the particle characteristics and the structural diversity of the isozymes. These findings may represent additional mechanisms for the differential inhibitory effects arising from the coincubated AuNPs on the metabolic activities of the hepatic CYP isozymes.

A large-scale allosteric transition in cytochrome P450 3A4 revealed by luminescence resonance energy transfer (LRET)

Effector-induced allosteric transitions in cytochrome P450 3A4 (CYP3A4) were investigated by luminescence resonance energy transfer (LRET) between two SH-reactive probes attached to various pairs of distantly located cysteine residues, namely the double-cysteine mutants CYP3A4(C64/C468), CYP3A4(C377/C468) and CYP3A4(C64/C121). Successive equimolar labeling of these proteins with the phosphorescent probe erythrosine iodoacetamide (donor) and the near-infrared fluorophore DY-731 maleimide (acceptor) allowed us to establish donor/acceptor pairs sensitive to conformational motions. The interactions of all three double-labeled mutants with the allosteric activators α-naphthoflavone and testosterone resulted in an increase in the distance between the probes. A similar effect was elicited by cholesterol. These changes in distance vary from 1.3 to 8.5 Å, depending on the position of the donor/acceptor pair and the nature of the effector. In contrast, the changes in the interprobe distance caused by such substrates as bromocriptine or 1-pyrenebutanol were only marginal. Our results provide a decisive support to the paradigm of allosteric modulation of CYP3A4 and indicate that the conformational transition caused by allosteric effectors increases the spatial separation between the beta-domain of the enzyme (bearing residues Cys₆₄ and Cys₃₇₇) and the alpha-domain, where Cys₁₂₁ and Cys₄₆₈ are located.

Functional characterization of cytochromes P450 2B from the desert woodrat Neotoma lepida

Mammalian detoxification processes have been the focus of intense research, but little is known about how wild herbivores process plant secondary compounds, many of which have medicinal value or are drugs. cDNA sequences that code for three enzymes of the cytochrome P450 (CYP) 2B subfamily, here termed 2B35, 2B36, and 2B37 have been recently identified from a wild rodent, the desert woodrat (Malenke et al., 2012). Two variant clones of each enzyme were engineered to increase protein solubility and to facilitate purification, as reported for CYP2B enzymes from multiple species. When expressed in Escherichia coli each of the woodrat proteins gave the characteristic maximum at 450nm in a reduced carbon monoxide difference spectrum but generally expressed at lower levels than rat CYP2B1. Two enzymes, 2B36 and 2B37, showed dealkylation activity with the model substrates 7-ethoxy-4-(trifluoromethyl)coumarin and 7-benzyloxyresorufin, whereas 2B35 was inactive. Binding of the monoterpene (+)-α-pinene produced a Type I shift in the absorbance spectrum of each enzyme. Mutation of 2B37 at residues 114, 262, or 480, key residues governing ligand interactions with other CYP2B enzymes, did not significantly change expression levels or produce the expected functional changes. In summary, two catalytic and one ligand-binding assay are sufficient to distinguish among CYP2B35, 2B36, and 2B37. Differences in functional profiles between 2B36 and 2B37 are partially explained by changes in substrate recognition site residue 114, but not 480. The results advance our understanding of the mechanisms of detoxification in wild mammalian herbivores and highlight the complexity of this system.

Prediction of amino acid residues participated in substrate recognition by cytochrome P450 subfamilies with broad substrate specificity

Cytochromes P450 comprise a large superfamily and several of their isoforms play a crucial role in metabolism of xenobiotics, including drugs. Although these enzymes demonstrate broad and cross-substrate specificity, different cytochrome P450 subfamilies exhibit certain selectivity for some types of substrates. Analysis of amino acid residues of the active sites of six cytochrome subfamilies (CYP1А, CYP2А, CYP2С, CYP2D, CYP2E and CYP3А) enables to define subfamily-specific patterns that consist of four residues. These residues are located on the periphery of the active sites of these cytochromes. We suggest that they can form a primary binding site at the entrance to the active site, defining cytochrome substrate recognition.

Structural and thermodynamic basis of (+)-α-pinene binding to human cytochrome P450 2B6

Despite recent advances in atomic-level understanding of drug and inhibitor interactions with human cytochromes P450, the decades-old questions of chemical and structural determinants of hydrocarbon binding are still unanswered. (+)-α-Pinene is a monoterpene hydrocarbon that is widely distributed in the environment and a potent P450 2B inhibitor. Therefore, a combined biophysical and structural analysis of human P450 2B6 interactions with (+)-α-pinene was undertaken to elucidate the basis of the very high affinity binding. Binding of (+)-α-pinene to the P450 active site was demonstrated by a Type I spectral shift. Thermodynamics of ligand binding were explored using isothermal titration calorimetry and compared to those of P450 2A6, which is much less flexible than 2B6 based on comparison of multiple X-ray crystal structures. Consistent with expectation, entropy is the major driving force for hydrocarbon binding to P450 2A6, as evidenced by the calorimetric results. However, formation of the 2B6-(+)-α-pinene complex has a significant enthalpic component. A 2.0 Å resolution crystal structure of this enzyme-ligand complex reveals that the highly plastic 2B6 utilizes previously unrecognized rearrangements of protein motifs. The results indicate that the specific components of enthalpic contribution to ligand binding are closely tied to the degree of enzyme flexibility.

A structural snapshot of CYP2B4 in complex with paroxetine provides insights into ligand binding and clusters of conformational states

An X-ray crystal structure of CYP2B4 in complex with the drug paroxetine [(3S,4R)-3-[(2H-1,3-benzodioxol-5-yloxy)methyl]-4-(4-fluorophenyl)piperidine] was solved at 2.14 Å resolution. The structure revealed a conformation intermediate to that of the recently solved complex with amlodipine and that of the more compact complex with 4-(4-chlorophenyl)imidazole in terms of the placement of the F-G cassette. Moreover, comparison of the new structure with 15 previously solved structures of CYP2B4 revealed some new insights into the determinants of active-site size and shape. The 2B4-paroxetine structure is nearly superimposable on a previously solved closed structure in a ligand-free state. Despite the overall conformational similarity among multiple closed structures, the active-site cavity volume of the paroxetine complex is enlarged. Further analysis of the accessible space and binding pocket near the heme reveals a new subchamber that resulted from the movement of secondary structural elements and rearrangements of active-site side chains. Overall, the results from the comparison of all 16 structures of CYP2B4 demonstrate a cluster of protein conformations that were observed in the presence or absence of various ligands.

Structural insights into aldosterone synthase substrate specificity and targeted inhibition

Aldosterone is a major mineralocorticoid hormone that plays a key role in the regulation of electrolyte balance and blood pressure. Excess aldosterone levels can arise from dysregulation of the renin-angiotensin-aldosterone system and are implicated in the pathogenesis of hypertension and heart failure. Aldosterone synthase (cytochrome P450 11B2, CYP11B2) is the sole enzyme responsible for the production of aldosterone in humans. Blocking of aldosterone synthesis by mediating aldosterone synthase activity is thus a recently emerging pharmacological therapy for hypertension, yet a lack of structural information has limited this approach. Here, we present the crystal structures of human aldosterone synthase in complex with a substrate deoxycorticosterone and an inhibitor fadrozole. The structures reveal a hydrophobic cavity with specific features associated with corticosteroid recognition. The substrate binding mode, along with biochemical data, explains the high 11β-hydroxylase activity of aldosterone synthase toward both gluco- and mineralocorticoid formation. The low processivity of aldosterone synthase with a high extent of intermediates release might be one of the mechanisms of controlled aldosterone production from deoxycorticosterone. Although the active site pocket is lined by identical residues between CYP11B isoforms, most of the divergent residues that confer additional 18-oxidase activity of aldosterone synthase are located in the I-helix (vicinity of the O(2) activation path) and loops around the H-helix (affecting an egress channel closure required for retaining intermediates in the active site). This intrinsic flexibility is also reflected in isoform-selective inhibitor binding. Fadrozole binds to aldosterone synthase in the R-configuration, using part of the active site cavity pointing toward the egress channel. The structural organization of aldosterone synthase provides critical insights into the molecular mechanism of catalysis and enables rational design of more specific antihypertensive agents.

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.

Plasticity of CYP2B enzymes: structural and solution biophysical methods

In the past three years, major advances in understanding cytochrome P450 2B (CYP2B) structure-function relationships have been made through determination of multiple ligand-bound and one ligand-free X-ray crystal structure of CYP2B4 and one ligand-bound X-ray crystal structure of CYP2B6. These structures have provided insight into the features that provide the high degree of plasticity of the enzymes. A combination of a phenylalanine cluster that allows for concerted movement of helices F through G and a conserved set of electrostatic interactions involving Arg(262) facilitates movement of this region to accommodate binding of ligands of various sizes without perturbing most of the P450 fold. Integrating solution based techniques such as NMR or deuterium exchange mass spectrometry (DXMS) with computational methods including molecular docking has provided further insight into enzyme behavior upon ligand binding. In addition, extended molecular dynamics simulations have provided a link between an open and a closed conformation of ligand-free CYP2B4 found in crystal structures. Other studies revealed the utility of rational engineering in improving stability of P450s to facilitate structural studies. The solution and computational results combined with the X-ray crystal structures yield a comprehensive picture of how these enzymes adopt different conformations to bind various ligands.

Predictive models for cytochrome p450 isozymes based on quantitative high throughput screening data

The human cytochrome P450 (CYP450) isozymes are the most important enzymes in the body to metabolize many endogenous and exogenous substances including environmental toxins and therapeutic drugs. Any unnecessary interactions between a small molecule and CYP450 isozymes may raise a potential to disarm the integrity of the protection. Accurately predicting the potential interactions between a small molecule and CYP450 isozymes is highly desirable for assessing the metabolic stability and toxicity of the molecule. The National Institutes of Health Chemical Genomics Center (NCGC) has screened a collection of over 17,000 compounds against the five major isozymes of CYP450 (1A2, 2C9, 2C19, 2D6, and 3A4) in a quantitative high throughput screening (qHTS) format. In this study, we developed support vector classification (SVC) models for these five isozymes using a set of customized generic atom types. The CYP450 data sets were randomly split into equal-sized training and test sets. The optimized SVC models exhibited high predictive power against the test sets for all five CYP450 isozymes with accuracies of 0.93, 0.89, 0.89, 0.85, and 0.87 for 1A2, 2C9, 2C19, 2D6, and 3A4, respectively, as measured by the area under the receiver operating characteristic (ROC) curves. The important atom types and features extracted from the five models are consistent with the structural preferences for different CYP450 substrates reported in the literature. We also identified novel features with significant discerning power to separate CYP450 actives from inactives. These models can be useful in prioritizing compounds in a drug discovery pipeline or recognizing the toxic potential of environmental chemicals.

Structures of cytochrome P450 2B6 bound to 4-benzylpyridine and 4-(4-nitrobenzyl)pyridine: insight into inhibitor binding and rearrangement of active site side chains

The biochemical, biophysical, and structural analysis of the cytochrome P450 2B subfamily of enzymes has provided a wealth of information regarding conformational plasticity and substrate recognition. The recent X-ray crystal structure of the drug-metabolizing P450 2B6 in complex with 4-(4-chlorophenyl)imidazole (4-CPI) yielded the first atomic view of this human enzyme. However, knowledge of the structural basis of P450 2B6 specificity and inhibition has remained limited. In this study, structures of P450 2B6 were determined in complex with the potent inhibitors 4-benzylpyridine (4-BP) and 4-(4-nitrobenzyl)pyridine (4-NBP). Comparison of the present structures with the previous P450 2B6-4-CPI complex showed that reorientation of side chains of the active site residue Phe206 on the F-helix and Phe297 on the I-helix was necessary to accommodate the inhibitors. However, P450 2B6 does not require any major side chain rearrangement to bind 4-NBP compared with 4-BP, and the enzyme provides no hydrogen-bonding partners for the polar nitro group of 4-NBP within the hydrophobic active site. In addition, on the basis of these new structures, substitution of residue 172 with histidine as observed in the single nucleotide polymorphism Q172H and in P450 2B4 may contribute to a hydrogen bonding network connecting the E- and I-helices, thereby stabilizing active site residues on the I-helix. These results provide insight into the role of active site side chains upon inhibitor binding and indicate that the recognition of the benzylpyridines in the closed conformation structure of P450 2B6 is based solely on hydrophobicity, size, and shape.

Structure and function of cytochromes P450 2B: from mechanism-based inactivators to X-ray crystal structures and back

This article reviews work from the author dating back to 1978 and focuses on the structural basis of cytochrome P450 (P450) function using available contemporary techniques. Early studies used mechanism-based inactivators that bound to the protein moiety of hepatic P450s to try to localize the active site. Subsequent studies used cDNA cloning, heterologous expression, site-directed mutagenesis, and homology modeling based on multiple bacterial P450 X-ray crystal structures to predict the active sites of CYP2B enzymes with considerable accuracy. Breakthroughs in engineering and expression of mammalian P450s enabled us to determine our first X-ray crystal structure of ligand-free rabbit CYP2B4. To date, we have solved 11 CYP2B4 and three human CYP2B6 structures, which represent four significantly different conformations. The plasticity of CYP2B4 has been confirmed by deuterium exchange mass spectrometry and is substantiated by molecular dynamics simulations. In addition to major movement of secondary structure elements, more subtle reorientation of active site side chains, especially Phe206, Phe297, and Glu301, contributes to the ability of CYP2B enzymes to bind various ligands. Isothermal titration calorimetry has proven to be a useful tool for studying the thermodynamics of ligand binding to CYP2B4 and CYP2B6, and NMR has enabled study of ligand binding orientation in solution as an adjunct to X-ray crystallography. A major challenge remains to harness the power of the various approaches to facilitate prediction of CYP2B specificity and inhibition.

Analysis of Cytochrome P450 Conserved Sequence Motifs between Helices E and H: Prediction of Critical Motifs and Residues in Enzyme Functions

Rational approaches have been extensively used to investigate the role of active site residues in cytochrome P450 (CYP) functions. However, recent studies using random mutagenesis suggest an important role for non-active site residues in CYP functions. Meta-analysis of the random mutants showed that 75% of the functionally important non-active site residues are present in 20% of the entire protein between helices E and H (E-H) and conserved sequence motif (CSM) between 7 and 11. The CSM approach was developed recently to investigate the functional role of non-active site residues in CYP2B4. Furthermore, we identified and analyzed the CSM in multiple CYP families and subfamilies in the E-H region. Results from CSM analysis showed that CSM 7, 8, 10, and 11 are conserved in CYP1, CYP2, and CYP3 families, while CSM 9 is conserved only in CYP2 family. Analysis of different CYP2 subfamilies showed that CYP2B and CYP2C have similar characteristics in the CSM, while the characteristics of CYP2A and CYP2D subfamilies are different. Finally, we analyzed CSM 7, 8, 10, and 11, which are common in all the CYP families/subfamilies analyzed, in fifteen important drug-metabolizing CYPs. The results showed that while CSM 8 is most conserved among these CYPs, CSM 7, 9, and 10 have significant variations. We suggest that CSM8 has a common role in all the CYPs that have been analyzed, while CSM 7, 10, and 11 may have relatively specific role within the subfamily. We further suggest that these CSM play important role in opening and closing of the substrate access/egress channel by modulating the flexible/plastic region of the protein. Thus, site-directed mutagenesis of these CSM can be used to study structure-function and dynamic/plasticity-function relationships and to design CYP biocatalysts.

Comparison of in vitro metabolism of ticlopidine by human cytochrome P450 2B6 and rabbit cytochrome P450 2B4

A recent X-ray crystal structure of a rabbit cytochrome P450 2B4 (CYP2B4)-ticlopidine complex indicated that the compound could be modeled with either the thiophene or chlorophenyl group oriented toward the heme prosthetic group. Subsequent NMR relaxation and molecular docking studies suggested that orientation with the chlorophenyl ring closer to the heme was the preferred one. To evaluate the predictive value of these findings, the oxidation of ticlopidine by reconstituted CYP2B4 was studied and compared with CYP2B6, in which the thiophene portion of the molecule likely orients toward the heme. In vitro incubation of ticlopidine with both enzymes yielded the same set of metabolites: 7-hydroxyticlopidine (M1), 2-oxoticlopidine (M2), 5-(2-chlorobenzyl)thieno[3,2-c]pyridin-5-ium metabolite (M3), 5-(2-chlorobenzyl)thieno[3,2-c]pyridin-5-ium metabolite (M4), ticlopidine N-oxide (M5), and ticlopidine S-oxide dimer, a dimerization product of ticlopidine S-oxide (M6). The rates of metabolite formation deviated markedly from linearity with time, consistent with the known inactivation of CYP2B6 by ticlopidine. Fitting to a first-order equation yielded similar rate constants (k(obs)) for both enzymes. However, the amplitude (R(max)) of M1 and M6 formation was 4 to 5 times higher for CYP2B6 than CYP2B4, indicating a greater residence time of ticlopidine with its thiophene ring closer to heme in CYP2B6. In contrast, CYP2B4 formed M4 and M5 in more abundance than CYP2B6, indicating an alternate orientation. Overall, the results suggest that the preferential orientation of ticlopidine in the active site of CYP2B4 predicted by X-ray crystallography and NMR studies is unproductive and that ticlopidine likely reorients within CYP2B4 to a more productive mode.