Quantitative structure-activity relationships (QSARs) within cytochromes P450 2B (CYP2B) subfamily enzymes: the importance of lipophilicity for binding and metabolism

A general linear free energy relationship for predicting partition coefficients of neutral organic compounds

Predicting the effects of organic compounds on environments and biological systems is an important objective for environmental chemistry and human health. The logarithm (to base 10) of the n-octanoll-water partition coefficient has been widely used to predict the mentioned properties. However, the suitability of this parameter for the purpose has been questioned, since the environments relating to the properties may be quite different from that of bulk n-octanol. In this study, we used a theoretical derivation approach to develop a model for predicting the partition coefficients of solutes between water and an organic solvent that may be similar to n-octanol or quite different from it. Our model relies on solute descriptors that can be calculated based on solute structures. It was used to predict the n-octanoll-water, hexadecanel-water and chloroforml-water partition coefficients of solutes. The calculated values of the examined parameters agreed well with their experimental counterparts. The model can find application in the accurate prediction of the effects of organic compounds on environments and the physicochemical properties of organic compounds by a full in-silico approach and can provide useful guidance for improving some properties of organic compounds.

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.

Molecular characterization of CYP2B6 substrates

CYP2B6 has not been as fully characterized at the molecular level as other members of the human cytochrome P450 family. As more widely used in vitro probes for characterizing the involvement of this enzyme in the metabolism of xenobiotics have become available, the number of molecules identified as CYP2B6 substrates has increased. In this study we have analyzed the available kinetic data generated by multiple laboratories with human recombinant expressed CYP2B6 and along with calculated molecular properties derived from the ChemSpider database, we have determined the molecular features that appear to be important for CYP2B6 substrates. In addition we have applied 2D and 3D QSAR methods to generate predictive pharmacophore and 2D models. For 28 molecules with K(m) data, the molecular weight (mean +/- SD) is 253.78+/-74.03, ACD/logP is 2.68+/-1.51, LogD(pH 5.5) is 1.51+/-1.43, LogD(pH 7.4) is 2.02+/-1.25, hydrogen bond donor (HBD) count is 0.57 +/-0.57, hydrogen bond acceptor (HBA) count is 2.57+/-1.37, rotatable bonds is 3.50+/-2.71 and total polar surface area (TPSA) is 27.63+/-19.42. A second set of 15 molecules without K(m) data possessed similar mean molecular property values. These properties are comparable to those of a set of 21 molecules used in a previous pharmacophore modeling study (Ekins et al., J Pharmacol Exp Ther 288 (1), 21-29, 1999). Only the LogD and HBD values were statistically significantly different between these different datasets. We have shown that CYP2B6 substrates are generally small hydrophobic molecules that are frequently central nervous system active, which may be important for drug discovery research.

Structure-activity relationships for the inhibition of recombinant human cytochromes P450 by curcumin analogues

Inhibition of cytochrome P450 (CYP) is a major cause of drug-drug interactions. In this work, inhibitory potentials of 33 curcumin analogues, i.e. 2,6-dibenzylidenecyclohexanone (A series), 2,5-dibenzylidenecyclopentanone (B series) and 1,4-pentadiene-3-one (C series) substituted analogues of curcumin towards recombinant human CYP1A2, CYP3A4, CYP2B6, CYP2C9 and CYP2D6, all important for drug metabolism, were studied in vitro. Fluorescence plate reader and high performance liquid chromatography (HPLC) assays were used to evaluate CYP-inhibitory activities. MOE-based Quantitative structure-activity relationship (QSAR) analysis suggested that electrostatic and hydrophobic interactions and lipophilicity are important factors for CYP inhibition. Apart from insights in important molecular properties for CYP inhibition, the present results may also guide further design of curcumin analogues with less susceptibility to drug-drug interactions.