ABSTRACT
The human cytochrome P450 enzymes (CYPs) are heme-protein monooxygenases, which catalyze oxidative reactions of a broad spectrum of substrates. Consequently, they play a critical role in the metabolism of xenobiotics, such as drugs and carcinogens, and the catabolism of endogenous lipophilic factors. Bioavailability and toxicity, both of which can be related to CYPs, continue to pose problems in the development of new drugs. The isoform which metabolizes over one-third of drugs, CYP 3A4, was investigated employing ensemble-docking experiments of a 195-substrate library with induced fit and GOLD docking algorithms and a number of scoring functions. Enzyme conformations included three currently available CYP 3A4 crystal structures. All docking experiments were performed in duplicates with and without inclusion of crystallographic waters. Resultant poses were assessed based on accuracy of site of metabolism prediction. Analyses of the docked solutions pertaining to ranking efficacy, ligand molecular properties, stabilizing residues in the ligand-enzyme complexes, and metabolic reactions are discussed. Our analyses suggest that certain residues make favorable interactions with the bound substrates. Employing multiple receptor conformations enhances the accuracy of catalytic prediction, while ligand size and flexibility impact docking performance. The presence of waters observed in crystal complexes does not necessarily lead to improved performance.
Subject(s)
Cytochrome P-450 CYP3A/metabolism , Pharmaceutical Preparations/metabolism , Binding Sites , Cytochrome P-450 CYP3A/chemistry , Drug Discovery , Humans , Ligands , Molecular Docking Simulation , Pharmaceutical Preparations/chemistry , Substrate SpecificityABSTRACT
It has been reported that somatostatin receptor subtypes 4 and 5 would be high-impact templates for homology modeling if their 3D structures became available. We have generated a homology model of the somatostatin receptor subtype 4 (sst4), using the newest active state ß(2) adrenoreceptor crystal structure, and subsequently docked a variety of agonists into the model-built receptor to elucidate the binding modes of reported agonists. Using experimental restraints, we were able to explain observed activity profiles. We propose two binding modes that can consistently explain findings for high-affinity agonists and reason why certain structures display low affinities for the receptor.
