Binding free energies for nicotine analogs inhibiting cytochrome P450 2A6 by a combined use of molecular dynamics simulations and QM/MM-PBSA calculations

Molecular dynamics (MD) simulations and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations have been performed to explore the dynamic behaviors of cytochrome P450 2A6 (CYP2A6) binding with nicotine analogs (that are typical inhibitors) and to calculate their binding free energies in combination with Poisson–Boltzmann surface area (PBSA) calculations. The combined MD simulations and QM/MM-PBSA calculations reveal that the most important structural parameters affecting the CYP2A6-inhibitor binding affinity are two crucial internuclear distances, that is, the distance between the heme iron atom of CYP2A6 and the coordinating atom of the inhibitor, and the hydrogen-bonding distance between the N297 side chain of CYP2A6 and the pyridine nitrogen of the inhibitor. The combined MD simulations and QM/MM-PBSA calculations have led to dynamic CYP2A6-inhibitor binding structures that are consistent with the observed dynamic behaviors and structural features of CYP2A6-inhibitor binding, and led to the binding free energies that are in good agreement with the experimentally-derived binding free energies. The agreement between the calculated binding free energies and the experimentally-derived binding free energies suggests that the combined MD and QM/MM-PBSA approach may be used as a valuable tool to accurately predict the CYP2A6-inhibitor binding affinities in future computational design of new, potent and selective CYP2A6 inhibitors.     

Targeting of Splice Variants of Human Cytochrome P450 2C8 (CYP2C8) to Mitochondria and Their Role in Arachidonic Acid Metabolism and Respiratory Dysfunction

In this study, we found that the full-length CYP2C8 (WT CYP2C8) and N-terminal truncated splice variant 3 (∼44-kDa mass) are localized in mitochondria in addition to the endoplasmic reticulum. Analysis of human livers showed that the mitochondrial levels of these two forms varied markedly. Molecular modeling based on the x-ray crystal structure coordinates of CYP2D6 and CYP2C8 showed that despite lacking the N-terminal 102 residues variant 3 possessed nearly complete substrate binding and heme binding pockets. Stable expression of cDNAs in HepG2 cells showed that the WT protein is mostly targeted to the endoplasmic reticulum and at low levels to mitochondria, whereas variant 3 is primarily targeted to mitochondria and at low levels to the endoplasmic reticulum. Enzyme reconstitution experiments showed that both microsomal and mitochondrial WT CYP2C8 efficiently catalyzed paclitaxel 6-hydroxylation. However, mitochondrial variant 3 was unable to catalyze this reaction possibly because of its inability to stabilize the large 854-Da substrate. Conversely, mitochondrial variant 3 catalyzed the metabolism of arachidonic acid into 8,9-, 11,12-, and 14,15-epoxyeicosatrienoic acids and 20-hydroxyeicosatetraenoic acid when reconstituted with adrenodoxin and adrenodoxin reductase. HepG2 cells stably expressing variant 3 generated higher levels of reactive oxygen species and showed a higher level of mitochondrial respiratory dysfunction. This study suggests that mitochondrially targeted variant 3 CYP2C8 may contribute to oxidative stress in various tissues.     

Evaluation of comparative cytochrome P450 2B4 model by photoaffinity labeling

A homology model of rabbit CYP 2B4 was constructed on the basis of the crystallographic structure of truncated mammalian CYP 2C5/3 and bacterial soluble CYPs. To validate the CYP 2B4 homology model photoaffinity labeling was employed. Three probes (I-III) containing a photo-labile azido-group and an amino-group on opposite ends of the molecule were designed for photoaffinity labeling of the CYP 2B4 in increasing distance from the heme iron. Spectroscopic data proved probes I (the shortest) and II (a middle sized) to be coordinated with the heme iron via their amino-groups in the enzyme active center while the probe III (the longest) was not bound in this way. This binding orientation of probes I and II is in accordance with the model predicting ion-pairing of the negatively charged side chain of CYP 2B4 Asp 105 and a positively charged nitrogen located in an appropriate position in structures of probes I and II, only. The lack of heme binding of the probe III is clear from its docking into the CYP 2B4 model since no Asp 105 ion-pairing is possible. The target of photoactivated probe II, Arg 197, in a distance of about 16.5 A from the heme iron, exactly matches the position of that amino acid residue, predicted from the CYP 2B4 homology model. Moreover, using this technique, a substrate access channel has been identified. To assess the predicted substrate-binding pocket, an interaction of a specific CYP 2B4 substrate, diamantane, was examined. In "silico" docking revealed strong binding of diamantane in an orientation allowing experimentally observed C4-hydroxylation. Our homology model of CYP 2B4 is thus consistent with experimental metabolic and photoaffinity labeling data.     

Theoretical study of the ligand-CYP2B4 complexes: effect of structure on binding free energies and heme spin state

The molecular origins of temperature-dependent ligand-binding affinities and ligand-induced heme spin state conversion have been investigated using free energy analysis and DFT calculations for substrates and inhibitors of cytochrome P450 2B4 (CYP2B4), employing models of CYP2B4 based on CYP2C5(3LVdH)/CYP2C9 crystal structures, and the results compared with experiment. DFT calculations indicate that large heme-ligand interactions (ca. -15 kcal/mol) are required for inducing a high to low spin heme transition, which is correlated with large molecular electrostatic potentials (approximately -45 kcal/mol) at the ligand heteroatom. While type II ligands often contain oxygen and nitrogen heteroatoms that ligate heme iron, DFT results indicate that BP and MF heme complexes, with weak substrate-heme interactions (ca. -2 kcal/mol), and modest MEPS minima (>-35 kcal/mol) are high spin. In contrast, heme complexes of the CYP2B4 inhibitor, 4PI, the product of benzphetamine metabolism, DMBP, and water are low spin, have substantial heme-ligand interaction energies (<-15 kcal/mol) and deep MEPS minima (<-45 kcal/mol) near their heteroatoms. MMPBSA analysis of MD trajectories were made to estimate binding free energies of these ligands at the heme binding site of CYP2B4. In order to initially assess the realism of this approach, the binding free energy of 4PI inhibitor was computed and found to be a reasonable agreement with experiment: -7.7 kcal/mol [-7.2 kcal/mol (experiment)]. BP was determined to be a good substrate [-6.3 kcal/mol (with heme-ligand water), -7.3 kcal/mol (without ligand water)/-5.8 kcal/mol (experiment)], whereas the binding of MF was negligible, with only marginal binding binding free energy of -1.7 kcal/mol with 2-MF bound [-3.8 kcal/mol (experiment)], both with and without retained heme-ligand water. Analysis of the free energy components reveal that hydrophobic/nonpolar contributions account for approximately 90% of the total binding free energy of these substrates and are the source of their differential and temperature-dependent CYP2B4 binding. The results indicate the underlying origins of the experimentally observed differential binding affinities of BP and MF, and indicate the plausibility of the use of models derived from moderate sequence identity templates in conjunction with approximate free energy methods in the estimation of ligand-P450 binding affinities.     

Molecular dynamics simulations of RNA kissing-loop motifs reveal structural dynamics and formation of cation-binding pockets

Explicit solvent molecular dynamics (MD) simulations were carried out for three RNA kissing–loop complexes. The theoretical structure of two base pairs (2 bp) complex of H3 stem–loop of Moloney murine leukemia virus agrees with the NMR structure with modest violations of few NMR restraints comparable to violations present in the NMR structure. In contrast to the NMR structure, however, MD shows relaxed intermolecular G-C base pairs. The core region of the kissing complex forms a cation-binding pocket with highly negative electrostatic potential. The pocket shows nanosecond-scale breathing motions coupled with oscillations of the whole molecule. Additional simulations were carried out for 6 bp kissing complexes of the DIS HIV-1 subtypes A and B. The simulated structures agree well with the X-ray data. The subtype B forms a novel four-base stack of bulged-out adenines. Both 6 bp kissing complexes have extended cation-binding pockets in their central parts. While the pocket of subtype A interacts with two hexacoordinated Mg²⁺ ions and one sodium ion, pocket of subtype B is filled with a string of three delocalized Na⁺ ions with residency times of individual cations 1–2 ns. The 6 bp complexes show breathing motions of the cation-binding pockets and loop major grooves.     

Exploration of the binding of curcumin analogues to human P450 2C9 based on docking and molecular dynamics simulation

Molecular docking and molecular dynamics (MD) simulations are used to investigate the interactions of curcumin analogues (CAs) with human cytochrome P450 2 C9 (CYP2C9 or 2 C9) and the conformations of their binding sites. In order to examine conformations of CAs/2 C9 and interaction characteristics of their binding sites, RMSDs, RMSFs, and B-factors are computed, and electrostatic and hydrophobic interactions between CAs and 2 C9 are analyzed and discussed. Results demonstrate that the most CAs studied lie 4~15 ? above the heme of CYP2C9. The presence of CAs makes some residues in bound CYP2C9s become more flexible. In the binding sites of A0/2 C9 and C0/2 C9, the formation of H-bond networks (or CA-water-residue bridges) enhances the interactions between CAs and 2 C9. The stronger inhibitory effects of A0, B0, and C0 on 2 C9 can be ascribed to stronger electrostatic and hydrophobic interactions in the binding sites of CAs/2 C9.     

Enzyme kinetic and molecular docking studies on the metabolic interactions of 1-hydroxy-2,3,5-trimethoxy-xanthone, isolated from Halenia elliptica D. Don, with model probe substrates of human cytochrome P450 enzymes

Halenia elliptica D. Don is a Tibetan herb and medicinal preparations containing Halenia elliptica have been commonly used for the treatment of hepatitis B virus infection in China. The metabolism of 1-hydroxy-2,3,5-trimethoxy-xanthone (HM-1) to its metabolites is mediated through cytochrome P450 enzymes. This study aimed to investigate the herb-drug interaction potential of HM-1 by studying its effects on the metabolism of model probe substrates of five major CYP450 isoforms in human liver microsomes. HM-1 showed moderate inhibitory effects on CYP1A2 (IC(50)=1.06μM) and CYP2C9 (IC(50)=3.89μM), minimal inhibition on CYP3A4 (IC(20)=11.94μM), but no inhibition on model CYP2D6 (dextromethorphan) and CYP2E1 (chlorzoxazone) probe substrates. Inhibition kinetic studies showed that the K(i) values of HM-1 on CYP1A2, CYP2C9 and CYP3A4 were 5.12μM, 2.00μM and 95.03μM, respectively. HM-1 competitively inhibited testosterone 6β-hydroxylation (CYP3A4) but displayed mixed type inhibitions for phenacetin O-deethylation (CYP1A2) and tolbutamide 4-hydroxylation (CYP2C9). Molecular docking study confirmed the inhibition modes of HM-1 on these human CYP isoforms.     

A comparison of aroclor 1254-induced and uninduced rat liver microsomes to human liver microsomes in phenytoin O-deethylation, coumarin 7-hydroxylation, tolbutamide 4-hydroxylation, S-mephenytoin 4'-hydroxylation, chloroxazone 6-hydroxylation and testosterone 6beta-hydroxylation

Aroclor 1254-induced rat liver homogenate supernatant (liver S-9) is routinely used as an exogenous metabolic activation system for the evaluation of mutagenicity of xenobiotics. The purpose of this study is to evaluate whether results obtained with Aroclor 1254-induced liver microsomes would be relevant to human. Aroclor 1254-induced and uninduced rat liver microsomes were compared to human liver microsomes in the metabolism of substrates which are known to be selectively metabolized by the major human cytochrome P450 (CYP) isoforms. The activities studied and the major CYP isoforms involved were as follows: phenacetin O-deethylation (CYP1A2); coumarin 7-hydroxylation, (CYP2A6); tolbutamide 4-hydroxylation (CYP2C9), S-mephenytoin 4'-hydroxylation (CYP2C19); dextromethorphan O-demethylation (CYP2D6); chloroxazone 6-hydroxylation (CYP2E1); and testosterone 6beta-hydroxylation (CYP3A4). We found that both induced and uninduced rat liver microsomes were active in all the pathways studied with the exception of coumarin 7-hydroxylation. Coumarin 7-hydroxylation was observed with human liver microsomes but not the rat liver microsomes. Aroclor-1254 was found to induce all activities measured, with the exception of coumarin 7-hydroxylation. Dextromethorphan O-deethylation activity was higher in the rat liver microsomes than the human liver microsomes. Testosterone 6beta-hydroxylation activity was found to be similar between the human liver microsomes and the induced rat liver microsomes. Our results suggest that experimental data obtained with Aroclor 1254-induced rat liver microsomes may not always be relevant to human.     

Characterization of inhibitory effects of perfluorooctane sulfonate on human hepatic cytochrome P450 isoenzymes: focusing on CYP2A6

Perfluorooctane sulfonate (PFOS) is a chemically stable compound extensively used as oil and water repellent, surface active agents in our daily life. Accumulative research evidence gradually appears the toxicity of PFOS against mammals, but the whole figure remains to be elucidated. The present study was conducted to know the effects of PFOS on human hepatic drug metabolizing-type cytochrome P450 (CYP) isoenzymes such as CYP1A2 (7-ethoxyresorufin as a substrate), CYP2A6 (coumarin), CYP2B6 (7-ethoxy-4-trifluoromethylcoumarin), CYP2C8 (paclitaxel), CYP2C9 (diclofenac), CYP2C19 (S-mephenytoin), CYP2D6 (bufuralol), CYP2E1 (chlorzoxazone) and CYP3A4 (testosterone) in human livers employing their typical substrates. Although all of the oxidation reactions tested were more or less inhibited by PFOS, diclofenac 4'-hydroxylation mediated mainly by CYP2C9 was most strongly inhibited (K(i) value of 40 nM), followed by paclitaxel 6α-hydroxylation mediated mainly by CYP2C8 (K(i) value of 4 μM). The substrate oxidation reactions catalyzed by CYP2A6, CYP2B6, CYP2C19 and CYP3A4 were moderately (K(i) values of 35 to 45 μM), and those by CYP1A2, CYP2D6 and CYP2E1 were weakly inhibited by PFOS (K(i) values of 190-300 μM). The inhibition by PFOS for coumarin 7-hydroxylation mainly catalyzed by human liver microsomal CYP2A6 as well as by the recombinant enzyme was found to be enhanced by the preincubation of PFOS with human liver microsomes and NADPH as compared to the case without preincubation. The inhibition of the human liver microsomal cumarin 7-hydroxylation was PFOS concentration-dependent, and exhibited pseudo-first-order kinetics with respect to preincubation time, yielding K(inact) and K(I) values of 0.06 min(-1) and 23 μM, respectively. These results suggest that the metabolism of medicines which are substrates for CYP2C9 may be altered by PFOS in human bodies, and that PFOS is a mechanism-based inhibitor of CYP2A6.     

CYP2A6*6, a novel polymorphism in cytochrome p450 2A6, has a single amino acid substitution (R128Q) that inactivates enzymatic activity

By using the polymerase chain reaction technique combined with restriction enzyme fragment length polymorphism (PCR-RFLP), a novel polymorphism of CYP2A6, CYP2A6*6, was detected in 0.4% of the Japanese population. To study the enzymatic properties of the CYP2A6.6 protein with a single amino acid substitution of arginine 128 to glutamine, both this isozyme and the CYP2A6.1 protein (wild-type) were produced in insect cells using a baculovirus system. Coumarin 7-hydroxylation, which reflects CYP2A6 activity, was significantly reduced (one-eighth of normal) in cell lysate from CYP2A6*6-transfected Sf9 cells compared with that lysate from CYP2A6*1-transfected cells. To clarify the mechanism of inactivation of the CYP2A6.6 enzyme, the heme content and reduced CO difference spectrum were examined. Although CYP2A6.6 retained about one-half the heme content of CYP2A6.1, the reduced CO-bound Soret peak was completely lost. These results suggest that the inactivation of CYP2A6.6 is mainly due to disordering of the holoprotein structure rather than a failure of heme incorporation.     

Analysis of clinically relevant substrates of CYP2B6 enzyme by computational methods

Mounting evidence thus far indicates that human cytochrome P450 2B6 (CYP2B6), an enzyme expressed at a relatively low level functionally, is primarily responsible for the metabolism of several clinically relevant drugs, including propofol, efavirenz, bupropion, mephobarbital, and the propofol analog 2,6-di-sec-butyl phenol. We used molecular dynamics and molecular docking methods to predict such interactions and to compare with experimentally measured metabolisms. Insight II and Discover Studio 2.5 were used to carry out the docking of these substrates into CYP2B6 to explore the critical residues and interaction energies of the complexes. Phe297, Glu301, Thr302 and Val367 were identified as major drug-binding residues, which is consistent with previous data on site-directed mutagenesis, crystallography structure, and from modeling and docking studies. In addition, our docking results suggest that nonpolar amino acid clusters and heme also participate in binding to mediate drug oxidative metabolism. The binding modes of the five clinically relevant substrates mentioned above for metabolism on CYP2B6 are presented.     

Molecular mechanism of ligand recognition by NR3 subtype glutamate receptors

NR3 subtype glutamate receptors have a unique developmental expression profile, but are the least well-characterized members of the NMDA receptor gene family, which have key roles in synaptic plasticity and brain development. Using ligand binding assays, crystallographic analysis, and all atom MD simulations, we investigate mechanisms underlying the binding by NR3A and NR3B of glycine and D-serine, which are candidate neurotransmitters for NMDA receptors containing NR3 subunits. The ligand binding domains of both NR3 subunits adopt a similar extent of domain closure as found in the corresponding NR1 complexes, but have a unique loop 1 structure distinct from that in all other glutamate receptor ion channels. Within their ligand binding pockets, NR3A and NR3B have strikingly different hydrogen bonding networks and solvent structures from those found in NR1, and fail to undergo a conformational rearrangement observed in NR1 upon binding the partial agonist ACPC. MD simulations revealed numerous interdomain contacts, which stabilize the agonist-bound closed-cleft conformation, and a novel twisting motion for the loop 1 helix that is unique in NR3 subunits.     

Effect of Conformational Dynamics on Substrate Recognition and Specificity as Probed by the Introduction of a de Novo Disulfide Bond into Cytochrome P450 2B1

The conformational dynamics of cytochrome P450 2B1 (CYP2B1) were investigated through the introduction of a disulfide bond to link the I- and K-helices by generation of a double Cys variant, Y309C/S360C. The consequences of the disulfide bonding were examined both experimentally and in silico by molecular dynamics simulations. Under high hydrostatic pressures, the partial inactivation volume for the Y309C/S360C variant was determined to be −21 cm³mol⁻¹, which is more than twice as much as those of the wild type (WT) and single Cys variants (Y309C, S360C). This result indicates that the engineered disulfide bond has substantially reduced the protein plasticity of the Y309C/S360C variant. Under steady-state turnover conditions, the S360C variant catalyzed the N-demethylation of benzphetamine and O-deethylation of 7-ethoxy-trifluoromethylcoumarin as the WT did, whereas the Y309C variant retained only 39% of the N-demethylation activity and 66% of the O-deethylation activity compared with the WT. Interestingly, the Y309C/S360C variant restored the N-demethylation activity to the same level as that of the WT but decreased the O-deethylation activity to only 19% of the WT. Furthermore, the Y309C/S360C variant showed increased substrate specificity for testosterone over androstenedione. Molecular dynamics simulations revealed that the engineered disulfide bond altered substrate access channels. Taken together, these results suggest that protein dynamics play an important role in regulating substrate entry and recognition.Liver microsomal cytochromes P450 (CYP or P450)² metabolize a large number of clinically used drugs that have diverse steric and functional moieties. Despite low sequence homology among CYPs from different families, all P450s invariably contain a heme cofactor that is coordinated to a thiolate and catalyze the oxidative metabolism, mostly through hydroxylation, of substrates. However, production of reactive intermediates by P450s is often associated with drug toxicity and carcinogenesis, and inhibition or induction of a specific P450 isoform may lead to adverse drug-drug interactions (1). From a clinical and pharmacological perspective, it is important to understand the structure, function, and dynamics of P450s.Structural studies of P450s by x-ray crystallography in the past decade have provided us with a wealth of information regarding the structural organization, critical active site residues, and proton delivery pathways of P450s (2–4). In particular, these structural analyses have consistently shown that certain regions of the P450 structures such as the F/G and B/B′-C loops are extremely flexible and can undergo large conformational changes to accommodate substrates of various sizes, although the overall folding pattern of all P450s is conserved. For instance, an open conformation was observed in the ligand-free CYP2B4 crystal structure, whereas a closed conformation was reported for the CPI-bound CYP2B4 (3, 5). The open-to-closed conformational change involves large motions of the F- and G-helices and the F/G and B/B′-C loops. Based on comparisons of the crystal structures of CYP2B4 bound with inhibitors of different sizes, Zhao et al. (6) identified five plastic regions in P450s, including the B/B′-C loop (PR2) and F/G loop (PR4). Binding of ketoconazole or erythromycin to CYP3A4 led to a large increase in the active site volume (>80% increase) because of conformational changes primarily in the PR4, but interestingly the F- and G-helices moved in the opposite direction (7). These authors proposed that the extreme flexibility of CYP3A4 accounts for its promiscuity, as CYP3A4 metabolizes nearly ∼50% of all clinically used drugs. The complexity of the conformational flexibility and dynamics are also revealed in an MD simulation study of CYP3A4, 2C9 and 2A6 (8). Importantly, this molecular dynamics (MD) simulation study shows that the three-dimensional structure of P450s is more flexible in solution than was observed in the crystal structure.Despite intensive studies of the crystal structures of microsomal P450s, insights into the conformational dynamics of P450s in solution, particularly in relation to their functional importance, are lacking. A laser flash photolysis study of CO rebinding to CYP2E1 in solution revealed that the binding of substrates such as ethanol, pyrazole, and acetaminophen restricts the conformational flexibility of CYP2E1, as the kinetics for the rebinding of CO to ligand-bound CYP2E1 are significantly slower than those for the ligand-free CYP2E1 (9). A solution thermodynamics study of CYP2B4 supports the notion that CYP2B4 is remarkably flexible, as the entropy substantially decreases upon inhibitor binding resulting from reduction of the hydrophobic surface (10). In this study, a de novo disulfide bond is engineered into CYP2B1 and the consequences resulting from the disulfide bonding are examined both experimentally and in silico using MD simulations. To discern the effect of the de novo disulfide bond apart from the Cys mutagenesis, both the single and double Cys variants were characterized in detail. To our knowledge, this is the first report that investigates the consequences of limiting conformational dynamics in a P450 by incorporating a disulfide bond. Our results demonstrate that protein dynamics play an important role in regulating substrate entry/product egress channels and substrate recognition and provide insights that will be valuable for rational drug design and protein engineering.     

Exploration of the binding of proton pump inhibitors to human P450 2C9 based on docking and molecular dynamics simulation

Human P450 protein CYP2C9 is one of the major drug-metabolizing isomers, contributing to the oxidation of 16% of the drugs currently in clinical use. To examine the interaction mechanisms between CYP2C9 and proton pump inhibitions (PPIs), we used molecular docking and molecular dynamics (MD) simulation methods to investigate the conformations and interactions around the binding sites of PPIs/CYPP2C9. Results from molecular docking and MD simulations demonstrate that nine PPIs adopt two different conformations (extended and U-bend structures) at the binding sites and position themselves far above the heme of 2C9. The presence of PPIs changes the secondary structures and residue flexibilities of 2C9. Interestingly, at the binding sites of all PPI–CYP2C9 complexes except for Lan/CYP2C9, there are hydrogen-bonding networks made of PPIs, water molecules, and some residues of 2C9. Moreover, there are strong hydrophobic interactions at all binding sites for PPIs/2C9, which indicate that electrostatic interactions and hydrophobic interactions appear to be important for stabilizing the binding sites of most PPIs/2C9. However, in the case of Lan/2C9, the hydrophobic interactions are more important than the electrostatic interactions for stabilizing the binding site. In addition, an interesting conformational conversion from extended to U-bend structures was observed for pantoprazole, which is attributed to an H-bond interaction in the binding pocket, an internal π–π stacking interaction, and an internal electrostatic interaction of pantoprazole.     

Prediction of sites of metabolism in a substrate molecule, instanced by carbamazepine oxidation by CYP3A4

In drug discovery process, improvement of ADME/Tox properties of lead compounds including metabolic stability is critically important. Cytochrome P450 (CYP) is one of the major metabolizing enzymes and the prediction of sites of metabolism (SOM) on the given lead compounds is key information to modify the compounds to be more stable against metabolism. There are two factors essentially important in SOM prediction. First is accessibility of each substrate atom to the oxygenated Fe atom of heme in a CYP protein, and the other is the oxidative reactivity of each substrate atom. To predict accessibility of substrate atoms to the heme iron, conventional protein-rigid docking simulations have been applied. However, the docking simulations without consideration of protein flexibility often lead to incorrect answers in the case of very flexible proteins such as CYP3A4. In this study, we demonstrated an approach utilizing molecular dynamics (MD) simulation for SOM prediction in which multiple MD runs were executed using different initial structures. We applied this strategy to CYP3A4 and carbamazepine (CBZ) complex. Through 10 ns MD simulations started from five different CYP3A4-CBZ complex models, our approach correctly predicted SOM observed in experiments. The experimentally known epoxidized sites of CBZ by CYP3A4 were successfully predicted as the most accessible sites to the heme iron that was judged from a numerical analysis of calculated ΔG(binding) and the frequency of appearance. In contrast, the predictions using protein-rigid docking methods hardly provided the correct SOM due to protein flexibility or inaccuracy of the scoring functions. Our strategy using MD simulation with multiple initial structures will be one of the reliable methods for SOM prediction.     

Exploring the structure characteristics and major channels of cytochrome P450 2A6, 2A13, and 2E1 with pilocarpine

The majority of cytochromes P450 play a critical role in metabolism of endogenous and exogenous substrates, some of its products are carcinogens. Therefore, inhibition of P450 enzymes activity can promote the detoxification and elimination of chemical carcinogens. In this study, molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed to explore the structure features and channel dynamics of three P450 isoforms 2A6, 2A13, and 2E1 bound with the common inhibitor pilocarpine. The binding free energy results combined with the PMF calculations give a reasonable ranking of binding affinity, which are consistent with the experimental data. Our results uncover how a sequence divergence of different CYP2 enzymes causes individual variations in major channel selections. On the basis of channel bottleneck and energy decomposition analysis, we propose a gating mechanism of their respective major channels in three enzymes, which may be attributed to a reversal of Phe209 in CYP2A6/2A13, as well as the rotation of Phe116 and Phe298 in CYP2E1. The hydrophobic residues not only make strong hydrophobic interactions with inhibitor, but also act as gatekeeper to regulate the opening of channel. The present study provides important insights into the structure–function relationships of three cytochrome P450s and the molecular basis for development of potent inhibitors.     

Molecular modelling of human CYP1B1 substrate interactions and investigation of allelic variant effects on metabolism

Molecular modelling of human CYP1B1 based on homology with the mammalian P450, CYP2C5, of known three-dimensional structure is reported. The enzyme model has been used to investigate the likely mode of binding for selected CYP1B1 substrates, particularly with regard to the possible effects of allelic variants of CYP1B1 on metabolism. In general, it appears that the CYP1B1 model is consistent with known substrate selectivity for the enzyme, and the sites of metabolism can be rationalized in terms of specific contacts with key amino acid residues within the CYP1B1 heme locus. Furthermore, a mode of binding interaction for the inhibitor, alpha-naphthoflavone, is presented which accords with currently available information. The current paper shows that a combination of molecular modelling and experimental determinations on the substrate metabolism for CYP1B1 allelic variants can aid in the understanding of structure-function relationships within P450 enzymes.     

Multiple molecular dynamics simulations of human p450 monooxygenase CYP2C9: the molecular basis of substrate binding and regioselectivity toward warfarin

To examine the molecular basis of activity and regioselectivity of the clinically important human microsomal cytochrome P450 (CYP) monooxygenase 2C9 toward its substrate warfarin, 22 molecular dynamics simulations (3-5 ns each) were performed in the presence and absence of warfarin. The resulting trajectories revealed a stable protein core and mobile surface elements. This mobility leads to the formation of two surface channels in the region between F-G loop, B' helix/B-B' loop, beta(1)-sheet, and between helices F and I and the turn in the C-terminal antiparallel beta-sheet in the presence of warfarin. Besides the nonproductive state of the CYP2C9 warfarin complex captured in the crystal structure, three additional states were observed. These states differ in the shape of the substrate binding cavity and the position of the warfarin molecule relative to heme. In one of these states, the 7- and 6-positions of warfarin contact the heme with a marked geometrical preference for position 7 over position 6. This modeling result is consistent with experimentally determined regioselectivity (71 and 22% hydroxylation in positions 7 and 6, respectively). Access to the heme group is limited by the core amino acids Ala297, Leu362, Leu366, and Thr301, which therefore are expected to have a major impact on regioselectivity. In addition, modeling predicts that autoactivation of warfarin is sterically hindered. Our study demonstrates how the combination of mobile surface and rigid core leads to interesting properties: a broad substrate profile and simultaneously a high regioselectivity.     

Structural Analysis of the Streptomyces avermitilis CYP107W1-Oligomycin A Complex and Role of the Tryptophan 178 Residue

CYP107W1 from Streptomyces avermitilis is a cytochrome P450 enzyme involved in the biosynthesis of macrolide oligomycin A. A previous study reported that CYP107W1 regioselectively hydroxylated C12 of oligomycin C to produce oligomycin A, and the crystal structure of ligand free CYP107W1 was determined. Here, we analyzed the structural properties of the CYP107W1-oligomycin A complex and characterized the functional role of the Trp178 residue in CYP107W1. The crystal structure of the CYP107W1 complex with oligomycin A was determined at a resolution of 2.6 Å. Oligomycin A is bound in the substrate access channel on the upper side of the prosthetic heme mainly by hydrophobic interactions. In particular, the Trp178 residue in the active site intercalates into the large macrolide ring, thereby guiding the substrate into the correct binding orientation for a productive P450 reaction. A Trp178 to Gly mutation resulted in the distortion of binding titration spectra with oligomycin A, whereas binding spectra with azoles were not affected. The Gly178 mutant’s catalytic turnover number for the 12-hydroxylation reaction of oligomycin C was highly reduced. These results indicate that Trp178, located in the open pocket of the active site, may be a critical residue for the productive binding conformation of large macrolide substrates.     

Three-dimensional model of the human aromatase enzyme and density functional parameterization of the iron-containing protoporphyrin IX for a molecular dynamics study of heme-cysteinato cytochromes

Mammalian cytochromes P450 (CYP) are enzymes of great biological and pharmaco-toxicological relevance. Due to their membrane-bound nature, the structural characterization of these proteins is extremely difficult, and therefore computational techniques, such as comparative modeling, may help obtaining reliable structures of members of this family. An important feature of CYP is the presence of an iron-containing porphyrin group at the enzyme active site. This calls for quantum chemical calculations to derive charges and parameters suitable for classical force field-based investigations of this proteins family. In this report, we first carried out density functional theory (DFT) computations to derive suitable charges for the Fe2+-containing heme group of P450 enzymes. Then, by means of the homology modeling technique, and taking advantage of the recently published crystal structure of the human CYP2C9, we built a new model of the human aromatase (CYP19) enzyme. Furthermore, to study the thermal stability of the new model as well as to test the suitability of the new DFT-based heme parameters, molecular dynamics (MD) simulations were carried out on both CYP2C9 and CYP19. Finally, the last few ns of aromatase MD trajectories were investigated following the essential dynamics protocol that allowed the detection of some correlated motions among some protein domains.