Molecular modeling of cytochrome P450 2B1: mode of membrane insertion and substrate specificity

Molecular modeling studies of L-arabinitol 4-dehydrogenase of Hypocrea jecorina: its binding interactions with substrate and cofactor

L-arabinitol 4-dehydrogenase (LAD1; EC 1.1.1.12) is an enzyme in the L-arabinose catabolic pathway of fungi that catalyzes the conversion of L-arabinitol into L-xylulose. The primary objective of this work is to identify the catalytic and coenzyme binding domains of LAD1 from Hypocrea jecorina in order to provide better insight into the possible catalytic events in these domains. The 3D structure of NAD(+)-dependent LAD1 was developed based on the crystal structure of human sorbitol dehydrogenase as a template. A series of molecular mechanics and dynamics operations were performed to find the most stable binding interaction for the enzyme and its ligands. Using the verified model, a docking study was performed with the substrate L-arabinitol, Zn(2+) and NAD(+). This study found a catalytic Zn(2+) binding domain (Cys66, His91, Glu92 and Glu176) and a cofactor NAD(+) binding domain (Gly202, ILeu204, Gly205, Cys273, Arg229 and Val298) with strong hydrogen bonding contacts with the substrate and cofactor. The binding pockets of the enzyme for l-arabinitol, NAD(+), and Zn(2+) have been explicitly defined. The results from this study should guide future mutagenesis studies and provide useful clues for engineering enzymes to improve the utilization of polyols for rare sugar production.

Structural requirements for inhibitors of cytochromes P450 2B: Assessment of the enzyme interaction with diamondoids

The series of diamondoids: adamantane, diamantane, triamantane, 2-isopropenyl-2-methyladamantane and 3-isopropenyl-3-methyldiamantane (3-IPMDIA), were employed to elucidate the molecular basis of their interaction with the active site of cytochromes P450 (CYP) of a 2B subfamily. These potent inhibitors of CYP2B enzymes were docked into the homology model of CYP2B4. Apparent dissociation constants calculated for the complexes of CYP2B4 with docked diamandoids agreed closely with the experimental data showing inhibition potency of the compounds and their binding affinity to CYP2B4. Superimposed structures of docked diamondoids mapped binding site residues. As they are mainly non-polar residues, the hydrophobicity plays the major role in the binding of diamondoids. Overlapping structure of diamondoids defined an elliptical binding cavity (5.9A inner diameter, 7.9A length) forming an angle of approximately 43 degrees with the heme plane. CYP2B specific diamondoids, namely 3-IPMDIA, showing the highest binding affinity, should be considered for a potential clinical use.

Homology modeling of cytochrome P450scc and the mutations for optimal amperometric sensor

A new homology model of bovine cytochrome P450scc is obtained starting from the recently determined crystal structure of mammalian cytochrome P450 2B4. The new emerging structure appears compatible with recent diffraction patterns of bovine P450scc microcrystals as obtained at the Microfocus Beamline of the European Synchrotron Radiation in Grenoble and here reported for the first time. The same atomic structure is utilized thereby to predict the mutations needed for modifying redox potential. A comprehensive comparison is finally carried out with the previous model present in the RCSB Protein DataBank also in terms of the alternative mutations being predicted for the same functional modification. The implication of these studies for optimal sensor construction is discussed.

Progress in cytochrome P450 active site modeling

Models capable of predicting the possible involvement of cytochromes P450 in the metabolism of drugs or drug candidates are important tools in drug discovery and development. Ideally, functional information would be obtained from crystal structures of all the cytochromes P450 of interest. Initially, only crystal structures of distantly related bacterial cytochromes P450 were available-comparative modeling techniques were used to bridge the gap and produce structural models of human cytochromes P450, and thereby obtain some useful functional information. A significant step forward in the reliability of these models came four years ago with the first crystal structure of a mammalian cytochrome P450, rabbit CYP2C5, followed by the structures of two human enzymes, CYP2C8 and CYP2C9, and a second rabbit enzyme, CYP2B4. The evolution of a CYP2D6 model, leading to the validation of the model as an in silico tool for predicting binding and metabolism, is presented as a case study.

The functional role of threonine-205 in the mechanism-based inactivation of P450 2B1 by two ethynyl substrates: the importance of the F helix in catalysis

We have previously demonstrated that substituting Val for Thr-205 in P450 2B1 abolishes the 16beta-hydroxylation of testosterone and markedly decreases the ability of 2-ethnylnaphthalene (2EN) and 17alpha-ethynylestradiol (17EE) to inactivate P450 2B1. The role of Thr-205 has been further investigated by measuring the kinetics of the mechanism-based inactivation of the 7-ethoxy-(trifluoromethyl)coumarin deethylation activity of 2B1 by 2EN and 17EE in wild-type (WT) and mutant P450s. In general, the kinetics of the inactivation of the Ser and Ala mutants was not significantly altered compared with WT. In contrast, the efficiency of the inactivation of the Val mutant decreased by approximately 6- and approximately 30-fold for 2EN and 17EE, respectively. High-pressure liquid chromatography (HPLC) analysis and SDS gel electrophoresis demonstrated the covalent binding of radiolabeled 2EN- and 17EE-reactive intermediates to the WT apoprotein, but not the Val mutant. The Val mutant was able to metabolize 2EN to 2-naphthylacetic acid, except the initial rate was slower than the WT. HPLC analysis of the 17EE incubation mixtures revealed three major metabolites and showed a correlation between the efficiency of inactivation and the generation of one of the major metabolites (C). Metabolite C was generated by the WT, Ser mutant, and Ala mutant. Metabolite C may be formed by the oxidation of the ethynyl group, and this reactive intermediate contributes to the inactivation of P450 2B1 by 17EE. The site-specific mutation of one residue, Thr-205 to Val, is sufficient to alter the profile of products formed during 17EE metabolism, such that very low levels of metabolite C are formed and inactivation is essentially abolished.

Threonine-205 in the F helix of p450 2B1 contributes to androgen 16 beta-hydroxylation activity and mechanism-based inactivation

Four mutants of Thr-205 in cytochrome p450 2B1 were constructed and expressed in Escherichia coli. The Ser-, Ala-, and Val-mutants displayed stable reduced CO difference spectra and were able to metabolize 7-ethoxy-4-(trifluoromethyl)coumarin, testosterone, androstenedione, and benzphetamine. The Arg-mutant displayed an unstable reduced CO difference spectrum at 450 nm, was concomitantly converted to a denatured form with a peak at 422 nm, and showed no catalytic activity with any of the four substrates tested. The Ser-mutant displayed activity and metabolite profiles for testosterone and androstenedione similar to those of the wild-type p450 2B1 (WT). Substitution of Thr-205 with Ala or Val markedly suppressed the 16 beta-hydroxylation activity but exhibited little effect on the 16 alpha-hydroxylation activity for testosterone and androstenedione. Because 16 beta-hydroxylation activity of androgens is a specific p450 2B subfamily marker and residue 205 is located in the F helix, which forms the ceiling of the active site, we postulate that the gamma-hydroxyl side chain of Thr may play an important role in directing the 16 beta-face of testosterone and androstenedione toward the active site. Surprisingly, the Val-mutant retained full activity for benzphetamine demethylation. When mechanism-based inactivators for p450 2B1 were used to evaluate the susceptibility to inactivation, the Val-mutant was resistant to inactivation by 17 alpha-ethynylestradiol and less sensitive to inactivation by 2-ethynylnaphthalene compared with the WT enzyme. Our results demonstrate the importance of Thr-205 in determining substrate specificity and product formation as well as in influencing the susceptibility of p450 2B1 to mechanism-based inactivators.

Analysis of differential substrate selectivities of CYP2B6 and CYP2E1 by site-directed mutagenesis and molecular modeling

Human CYP2B6 and CYP2E1 were used to investigate the extent to which differential substrate selectivities between cytochrome P450 subfamilies reflect differences in active-site residues as opposed to distinct arrangement of the backbone of the enzymes. Reciprocal CYP2B6 and CYP2E1 mutants at active-site positions 103, 209, 294, 363, 367, and 477 (numbering according to CYP2B6) were characterized using the CYP2B6-selective substrate 7-ethoxy-4-trifluoromethylcoumarin, the CYP2E1-selective substrate p-nitrophenol, and the common substrates 7-ethoxycoumarin, 7-butoxycoumarin, and arachidonic acid. This report is the first to study the active site of CYP2E1 by systematic site-directed mutagenesis. One of the most intriguing findings was that substitution of CYP2E1 Phe-477 with valine from CYP2B6 resulted in significant 7-ethoxy-4-trifluoromethylcoumarin deethylation. Use of three-dimensional models of CYP2B6 and CYP2E1 based on the crystal structure of CYP2C5 suggested that deethylation of 7-ethoxy-4-trifluoromethylcoumarin by CYP2E1 is impeded by van der Waals overlaps with the side chain of Phe-477. Interestingly, none of the CYP2B6 mutants acquired enhanced ability to hydroxylate p-nitrophenol. Substitution of residue 363 in CYP2E1 and CYP2B6 resulted in significant alterations of the metabolite profile for the side chain hydroxylation of 7-butoxycoumarin. Probing of CYP2E1 mutants with arachidonic acid indicated that residues Leu-209 and Phe-477 are critical for substrate orientation in the active site. Overall, the study revealed that differences in the side chains of active-site residues are partially responsible for differential substrate selectivities across cytochrome P450 subfamilies. However, the relative importance of active-site residues appears to be dependent on the structural similarity of the compound to other substrates of the enzyme.

Quantitative structure-activity relationship (QSAR) and three-dimensional QSAR analysis of a series of xanthates as inhibitors and inactivators of cytochrome P450 2B1

1. Various xanthates (R-OCS2) were found to be mechanism-based inactivators of cytochrome P450 2B1 (CYP2B1) and CYP2B6 via formation of reactive metabolites. 2. In the present study, quantitative structure-activity relationships (QSARs) were derived with inhibitory and inactivation potencies of 15 xanthates (R = two to 20 methylene groups, allyl, cyclohexyl or O-tricyclo[5.2.1.0(2,6)]dec-9-yl (D609)) against purified, reconstituted rat liver CYP2B1. Factor, regression and comparative molecular field analyses (CoMFA) were used. 3. The compounds formed two groups whose activities depended on different structural features: the first group consisted of compounds with ethyl, propyl, allyl, cyclohexyl and D609 substituents; the second involved compounds with eight to 20 methylene groups. 4. High correlation between the molecular volume and inhibitory potency of the xanthates of the second group was found. The inactivation potency in the first group correlated with the charge of the first carbon atom of R, identifying this atom as a potential target for metabolic attack. A decrease in the inactivation potency with an increase in the size of R was observed in the second group. This finding could be explained by a decreased rate of metabolism of the long alkyl chain compounds and/or by difficulty in binding of the resulting metabolite(s) to the enzyme molecule.

Comparative modelling of cytochromes P450

The superfamily of enzymes known as the cytochromes P450 (P450s) comprises a wide-ranging class of proteins with diverse functions. They are known, amongst other things, to be involved in the hormonal regulation of metabolism and reproduction, as well as having a major clinical significance through their association with diseases such as cancer, diabetes and hepatitis. Knowledge of the three-dimensional (3D) structure of a protein gives insight into its function. The 3D structures of P450s are therefore of considerable scientific interest. A number of high-resolution structures of P450s have been determined by X-ray crystallography and studies of these structures have provided valuable insights into the mechanism of these enzymes. Only one of these structures is mammalian and as yet there is no structural information on human P450s in the public domain. Until such a structure is solved it is necessary to employ alternative methods to gain structural insight into how human P450s perform their biological function. Here we report on the use of comparative modelling to predict the structure of human P450s based on knowledge of their amino acid sequences plus the 3D structures of other (not human) P450s. As an illustrative example of these techniques we have modelled the structure of P450 2C5 using five bacterial P450 structures as templates. We examine the importance of selecting suitable templates, obtaining a good amino acid sequence alignment, and evaluating the models generated. To improve the quality of the models an iterative cycle of sequence alignment, model building, and model evaluation is employed. The result is a model with excellent stereochemistry, good amino acid side chain environment properties, and a Calpha trace similar to the crystal structure.

A truncation of 2B subfamily cytochromes P450 yields increased expression levels, increased solubility, and decreased aggregation while retaining function

The hydrophobic membrane-spanning domain in four cytochromes P450 2B was removed (Delta3-21) and several positive charges were substituted at the N-terminus to increase expression and solubility. Histidine residues were appended to the C-terminus to simplify purification. The truncated proteins were highly expressed in Escherichia coli, could be released from the membrane using high salt conditions, and were purified from this fraction to specific contents up to 19 nmol P450/mg protein using a single Ni(2+)-agarose column. Gel filtration revealed that truncated P450 2B1 forms a monodisperse solution of hexamers in the absence of detergent and >95% monomers in 0.25% sodium cholate. All truncated proteins, including human 2B6, were active with 7-ethoxy-4-trifluoromethylcoumarin, and truncated 2B1 was shown to retain the native regio- and stereospecificity of testosterone hydroxylation. These data demonstrate that modification of the N-terminus yields high levels of properly folded P450s 2B with increased solubility, which are suitable for functional and structural analysis.

Protein engineering of cytochromes P-450

The cytochromes P-450 are an immensely important superfamily of heme-containing enzymes. They catalyze the monooxygenation of an enormous range of substrates. In bacteria, cytochromes P-450 are known to catalyze the hydroxylation of environmentally significant substrates such as camphor, phenolic compounds and many herbicides. In eukaryotes, these enzymes perform key roles in the synthesis and interconversion of steroids, while in mammals hepatic cytochromes P-450 are vital for the detoxification of many drugs. As such, the cytochromes P-450 are of considerable interest in medicine and biotechnology and are obvious targets for protein engineering. The purpose of this article is to illustrate the ways in which protein engineering has been used to investigate and modify the properties of cytochromes P-450. Illustrative examples include: the manipulation of substrate selectivity and regiospecificity, the alteration of membrane binding properties, and probing the route of electron transfer.

Theoretical investigation of substrate specificity for cytochromes P450 IA2, P450 IID6 and P450 IIIA4

Three-dimensional models of the cytochromes P450 IA2, P450 IID6 and P450 IIIA4 were built by means of comparative modeling using the X-ray crystallographic structures of P450 CAM, P450 BM-3, P450 TERP and P450 ERYF as templates. The three cytochromes were analyzed both in their intrinsic structural features and in their interaction properties with fifty specific and non-specific substrates. Substrate/enzyme complexes were obtained by means of both automated rigid and flexible body docking. The comparative analysis of the three cytochromes and the selected substrates, in their free and bound forms, allowed for the building of semi-quantitative models of substrate specificity based on both molecular and intermolecular interaction descriptors. The results of this study provide new insights into the molecular determinants of substrate specificity for the three different eukaryotic P450 isozymes and constitute a useful tool for predicting the specificity of new compounds.

Synthetic peptide mimics of a predicted topographical interaction surface: the cytochrome P450 2B1 recognition domain for NADPH-cytochrome P450 reductase

In order to identify the cytochrome P450-binding domain for NADPH-cytochrome P450 reductase, synthetic peptide mimics of predicted surface regions of rat cytochrome P450 2B 1 were constructed and evaluated for inhibition of the P450-reductase interaction. A peptide corresponding to residues 116-134, which includes the C helix, completely inhibited reductase-mediated benzphetamine demethylation by purified P450 2B1. Replacement of Arg-125 by Glu yielded a noninhibitory peptide, suggesting that this residue significantly contributes to the reductase-P450 interaction. Additional P450 peptides were prepared which correspond to combinations of regions distant in primary sequence, but predicted to be spatially proximate. A peptide derived from segments of the C and L helices was a more potent inhibitor than peptides derived from either segment alone. This topographically designed peptide not only inhibited P450 2B1 in its purified form, but also when membrane-bound in rat liver microsomes. The peptide also inhibited microsomal aryl hydrocarbon hydroxylase, aniline hydroxylase, and erythromycin demethylase activities derived from other P450s. These results indicate that the C and L helices contribute to a reductase-binding site common to multiple P450s, and present a peptide mimic for this region that is useful for inhibition of P450-mediated microsomal activities.

Inhibition of human cytochrome P450 1A2 by flavones: a molecular modeling study

Cytochrome P450 1A2 metabolizes a number of important drugs, procarcinogens, and endogenous compounds. Several flavones, a class of phytochemicals consumed in the human diet, have been shown to differentially inhibit human P450 1A2-mediated methoxyresorufin demethylase. A molecular model of this P450 was constructed in order to elucidate the molecular basis of the P450-flavone interaction. Flavone and its 3,5,7-trihydroxy and 3,5,7-trimethoxy derivatives were docked into the active site to assess their mode of binding. The site is hydrophobic and includes several residues that hydrogen bond with substituents on the flavone nucleus. The binding interactions of these flavones in the modeled active side are consistent with their relative inhibitory potentials, namely 3,5,7-trihydroxylflavone > flavone > 3,5,7-trimethoxylflavone, toward P450 1A2-mediated methoxyresorufin demethylation.