Dataset Information

Ensemble generation and the influence of protein flexibility on geometric tunnel prediction in cytochrome P450 enzymes.

ABSTRACT: Computational prediction of ligand entry and egress paths in proteins has become an emerging topic in computational biology and has proven useful in fields such as protein engineering and drug design. Geometric tunnel prediction programs, such as Caver3.0 and MolAxis, are computationally efficient methods to identify potential ligand entry and egress routes in proteins. Although many geometric tunnel programs are designed to accommodate a single input structure, the increasingly recognized importance of protein flexibility in tunnel formation and behavior has led to the more widespread use of protein ensembles in tunnel prediction. However, there has not yet been an attempt to directly investigate the influence of ensemble size and composition on geometric tunnel prediction. In this study, we compared tunnels found in a single crystal structure to ensembles of various sizes generated using different methods on both the apo and holo forms of cytochrome P450 enzymes CYP119, CYP2C9, and CYP3A4. Several protein structure clustering methods were tested in an attempt to generate smaller ensembles that were capable of reproducing the data from larger ensembles. Ultimately, we found that by including members from both the apo and holo data sets, we could produce ensembles containing less than 15 members that were comparable to apo or holo ensembles containing over 100 members. Furthermore, we found that, in the absence of either apo or holo crystal structure data, pseudo-apo or -holo ensembles (e.g. adding ligand to apo protein throughout MD simulations) could be used to resemble the structural ensembles of the corresponding apo and holo ensembles, respectively. Our findings not only further highlight the importance of including protein flexibility in geometric tunnel prediction, but also suggest that smaller ensembles can be as capable as larger ensembles at capturing many of the protein motions important for tunnel prediction at a lower computational cost.

SUBMITTER: Kingsley LJ

PROVIDER: S-EPMC4067289 | biostudies-literature | 2014

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Ensemble generation and the influence of protein flexibility on geometric tunnel prediction in cytochrome P450 enzymes.

Kingsley Laura J LJ Lill Markus A MA

PloS one 20140623 6

Computational prediction of ligand entry and egress paths in proteins has become an emerging topic in computational biology and has proven useful in fields such as protein engineering and drug design. Geometric tunnel prediction programs, such as Caver3.0 and MolAxis, are computationally efficient methods to identify potential ligand entry and egress routes in proteins. Although many geometric tunnel programs are designed to accommodate a single input structure, the increasingly recognized impor ...[more]

PMID: 24956479

Similar Datasets

Project description:In silico docking studies and quantitative structure-activity relationship analysis of a number of in-house cytochrome P450 inhibitors have revealed important structural characteristics that are required for a molecule to function as a good inhibitor of P450 enzymes 1A1, 1A2, 2B1, and/or 2A6. These insights were incorporated into the design of pharmacophores used for a 2D search of the Chinese medicine database. Emodin, a natural anthraquinone isolated from Rheum emodi and known to be metabolized by cytochrome P450 enzymes, was one of the hits and was used as the lead compound. Emodin was found to inhibit P450s 1A1, 1A2, and 2B1 with IC(50) values of 12.25, 3.73, and 14.89 ?M, respectively. On the basis of the emodin molecular structure, further similarity searches of the PubChem and ZINC chemical databases were conducted resulting in the identification of 12 emodin analogues for testing against P450s 1A1-, 1A2-, 2B1-, and 2A6-dependent activities. 1-Amino-4-chloro-2-methylanthracene-9,10-dione (compound 1) showed the best inhibition potency for P450 1A1 with an IC(50) value of 0.40 ?M. 1-Amino-4-chloro-2-methylanthracene-9,10-dione (compound 1) and 1-amino-4-hydroxyanthracene-9,10-dione (compound 2) both inhibited P450 1A2 with the same IC(50) value of 0.53 ?M. In addition, compound 1 acted as a mechanism-based inhibitor of cytochrome P450s 1A1 and 1A2 with K(I) and K(inactivation) values of 5.38 ?M and 1.57 min(-1) for P450 1A1 and 0.50 ?M and 0.08 min(-1) for P450 1A2. 2,6-Di-tert-butyl-5-hydroxynaphthalene-1,4-dione (compound 8) directly inhibited P450 2B1 with good selectivity and inhibition potency (IC(50) = 5.66 ?M). Docking studies using the 3D structures of the enzymes were carried out on all of the compounds. The binding modes of these compounds revealed the structural characteristics responsible for their potency and selectivity. Compound 1, which is structurally similar to compound 2 with the presence of an amino group at position 1, showed a difference in the mechanism of inhibition toward P450s 1A1 and 1A2. The mechanism-based inhibition seen for compound 1 may be attributed to the presence of the methyl group at the 2-position, in close proximity to the amino group. Compound 2, which is otherwise similar, lacks that methyl moiety and did not show mechanism-based inhibition.

Dataset Information

Ensemble generation and the influence of protein flexibility on geometric tunnel prediction in cytochrome P450 enzymes.

Publications

Ensemble generation and the influence of protein flexibility on geometric tunnel prediction in cytochrome P450 enzymes.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets