Rational design, discovery, and synthesis of a novel series of potent growth hormone secretagogues.

In the joint experimental and computational efforts reported here to obtain novel chemical entities as growth hormone secretagogues (GHSs), a small database of peptides and non-peptides known to have GHS activity was used to generate and assess a 3D pharmacophore for this activity. This pharmacophore was obtained using a systematic and efficient procedure, "DistComp", developed in our laboratory. The 3D pharmacophore identified was then used to search 3D databases to explore chemical structures that could be novel GHSs. A number of these were chosen for synthesis and assessment of their ability to release growth hormone (GH) from rat pituitary cells. Among the compounds tested, those with a benzothiazepin scaffold were discovered with micromolar activity. To facilitate lead optimization, a second program, a site-dependent fragment QSAR procedure was developed. This program calculates a library of chemical and physical properties of "fragments" or chemical components in a known pharmacophore and determines which, if any, of these properties are important for the observed activity. The combined use of the 3D pharmacophore and the results of the site-dependent fragment QSAR analysis led to the discovery and synthesis of a novel series of potent GHSs, a number of which had nanomolar in vitro activity.

Influence of benzodiazepine binding site ligands on fear-conditioned contextual memory.

Eight compounds that bind to the benzodiazepine binding site on the gamma-amino butyric acid(A) (GABA(A)) receptor were assessed for their influence on contextual memory, an aspect of memory affected in various cognitive disorders including Alzheimer's disease. Using a Pavlovian fear-conditioning paradigm, each ligand was evaluated in C57Bl/6 mice in regards to its direct affect on contextual memory and whether the ligand could attenuate scopolamine-induced contextual memory impairment. Of the eight ligands tested, one impaired contextual memory (agonist), six attenuated scopolamine-induced contextual memory impairment (inverse agonists), and one antagonized the ability of an inverse agonist to attenuate scopolamine-induced contextual memory impairment. Hence, further demonstrating the bi-directional influence benzodiazepine binding site ligands are able to exert on memory modulation. This study serves as an initial starting point in the development of pharmacological tools to be used in deciphering how GABA(A) receptors influence contextual memory.

Characterization of novel benzodiazepine ligands in Spodotera frugiperda (Sf-9) insect cells.

The goals of the present work were to characterize the binding profile of nine benzodiazepine ligands in Spodotera frugiperda (Sf-9) insect cells expressing specific gamma aminobutyric acid (A) (GABA(A)) receptor subunit combinations and compare the affinities to those for the receptors in the rat cerebellum. Three recombinant baculovirus constructs, each harboring a different GABA(A) receptor subunit, were introduced into insect cells by simultaneous infection. Saturation and competition binding assays were carried out in membranes from Sf-9 cells infected with either alpha1beta2gamma2 or alpha6beta2gamma2 subunit combinations. The affinities of the ligands to the alpha1beta2gamma2 or alpha6beta2gamma2 receptors expressed in Sf-9 cells were similar to the affinities previously determined for the alpha1 or alpha6 subunit-containing GABA(A) receptors in the rat cerebellum, respectively, thus confirming the previously assigned receptor types in the cerebellum.

Differentiation of delta, mu, and kappa opioid receptor agonists based on pharmacophore development and computed physicochemical properties.

Compounds that bind with significant affinity to the opioid receptor types, delta, mu, and kappa, with different combinations of activation and inhibition at these three receptors could be promising behaviorally selective agents. Working on this hypothesis, the chemical moieties common to three different sets of opioid receptor agonists with significant affinity for each of the three receptor types delta, mu, or kappa were identified. Using a distance analysis approach, common geometric arrangements of these chemical moieties were found for selected delta, mu, or kappa opioid agonists. The chemical and geometric commonalities among agonists at each opioid receptor type were then compared with a non-specific opioid recognition pharmacophore recently developed. The comparison provided identification of the additional requirements for activation of delta, mu, and kappa opioid receptors. The distance analysis approach was able to clearly discriminate kappa-agonists, while global molecular properties for all compounds were calculated to identify additional requirements for activation of delta and mu receptors. Comparisons of the combined geometric and physicochemical properties calculated for each of the three sets of agonists allowed the determination of unique requirements for activation of each of the three opioid receptors. These results can be used to improve the activation selectivity of known opioid agonists and as a guide for the identification of novel selective opioid ligands with potential therapeutic usefulness.

Molecular determinants of non-specific recognition of delta, mu, and kappa opioid receptors.

Identification of the molecular determinants of recognition common to all three opioid receptors embedded in a single three-dimensional (3D) non-specific recognition pharmacophore has been carried out. The working hypothesis that underlies the computational study reported here is that ligands that bind with significant affinity to all three cloned opioid receptors, delta, mu, and kappa, but with different combinations of activation and inhibition properties at these receptors, could be promising behaviorally selective analgesics with diminished side effects. The study presented here represents the first step towards the rational design of such therapeutic agents. The common 3D pharmacophore developed for recognition of delta, mu, and kappa opioid receptors was based on the receptor affinities determined for 23 different opioid ligands that display no specificity for any of the receptor subtypes. The pharmacophore centers identified are a protonated amine, two hydrophobic groups, and the centroid of an aromatic group in a geometric arrangement common to all 23, non-specific, opioid ligands studied. Using this three-dimensional pharmacophore as a query for searching 3D structural databases, novel compounds potentially involved in non-specific recognition of delta, mu, and kappa opioid receptors were retrieved. These compounds can be valuable candidates for novel behaviorally selective analgesics with diminished or no side effects, and thus with potential therapeutic usefulness.

Characterization of benzodiazepine receptors in the cerebellum.

1. The goals of the work reported here were to further characterize benzodiazepine/GABA(A) (BDZR) receptor heterogeneity in the cerebellum and to measure the affinities and selectivities of structurally diverse benzodiazepines at each site identified. 2. Five chemical families were included in these studies. These were 1,4-benzodiazepines (flunitrazepam), imidazobenzodiazepines (RO15-1788 and RO15-4513 and RO16-6028), beta-carbolines (Abecarnil) and pyrazoloquinolines (CGS 8216, CGS 9895 and CGS 9896). 3. Saturation and competition binding assays were combined with powerful data analysis software developed in our laboratory. Among the capabilities of this software is the identification of multiple binding sites for a cold ligand using a non-selective labeled ligand that binds with equal, but high, affinity to all the binding sites 4. Saturation binding assays using either [3H]-RO15-1788 or [3H]-RO15-4513 revealed only one apparent binding site, with a higher affinity for RO15-4513 than for RO15-1788. However, using [3H]-RO15-4513 for the competition binding studies in the cerebellum, together with our data analysis software, led to the identification of two distinct binding sites with equal densities for the diverse benzodiazepines studied. 5. In rat cerebellum one of the sites identified corresponds to GABA(A) receptors exhibiting alpha1 subunit pharmacology and the other to GABA(A) receptors exhibiting alpha6 subunit pharmacology. In general, the diverse families of BDZR ligands studied had much lower affinities for the alpha6 containing receptors.

Unraveling the identity of benzodiazepine binding sites in rat hipppocampus and olfactory bulb.

The goals of the work reported here were (i) to identify distinct GABA(A)/benzodiazepine receptors in the rat hippocampus and olfactory bulb using receptor binding assays, and (ii) to determine the affinities and selectivities of benzodiazepine receptor ligands from structurally diverse chemical families at each site identified. These studies were aided by the use of software AFFINITY ANALYSIS SYSTEM, developed in our laboratory for analysis of receptor binding data that allows the determination of receptor heterogeneity using non-selective radioligands. Saturation binding assays using [3H]RO15-4513 (ethyl 8-azido-6-dihydro-5-methyl-6-oxo-4H-imidazo[1, 5-a]-[1,4]benzodiazepine-3-carboxylate) revealed two binding sites in each of these two tissues. The higher affinity site corresponds to alpha(5) subunit-containing GABA(A) receptor and the lower affinity site to a combination of alpha(1), alpha(2), and alpha(3) subunit-containing receptors. These results should be useful in the challenging task of identifying the various functional GABA(A) receptors in the central nervous system, and in providing a link between receptor affinities and in vivo activities of the GABA(A)/benzodiazepine receptor ligands studied.

Development of a 3D pharmacophore for nonspecific ligand recognition of alpha1, alpha2, alpha3, alpha5, and alpha6 containing GABA(A)/benzodiazepine receptors.

Transfected cells containing GABA(A) benzodiazepine receptors (BDZRs) have been utilized to systematically determine the affinity of ligands at alpha1, alpha2, alpha3, alpha5 and alpha6 subtypes in combination with beta2 and gamma2. All but a few of the ligands thus far studied have relatively high affinities for each of these alpha subtype receptors. Thus, these ligands must contain common stereochemical properties favorable for recognition by each of the subtype combinations. In the present work, such a common three-dimensional (3D) pharmacophore for recognition of alpha1, alpha2, alpha3, alpha5 and alpha6 containing GABA(A)/BDZRs types of receptors has been developed and assessed, using as a database receptor affinities measured in transfected cells for 27 diverse compounds. The 3D-recognition pharmacophore developed consists of three proton accepting groups, a hydrophobic group, and the centroid of an aromatic ring found in a common geometric arrangement in the 19 nonselective ligands used. Three tests were made to assess this pharmacophore: (i) Four low affinity compounds were used as negative controls, (ii) Four high affinity compounds, excluded from the pharmacophore development, were used as compounds for pharmacophore validation, (iii) The 3D pharmacophore was used to search 3D databases. The results of each of these types of assessments provided robust validation of the 3D pharmacophore. This 3D pharmacophore can now be used to discover novel nonselective ligands that could be activation selective at different behavioral end points. Additionally, it may serve as a guide in the design of more selective ligands, by determining if candidate ligands proposed for synthesis conform to this pharmacophore and selecting those that do not for further experimental assessment.

Determinants of recognition of ligands binding to benzodiazepine receptor/GABA(A) receptors initiating sedation.

Complementary behavioral and computational studies of 21 structurally diverse, gamma-amino butyric acid (GABA)(A) benzodiazepine receptor ligands that influence spontaneous locomotor activity have been performed in this work. This behavioral endpoint is a well-accepted indicator of sedation particularly for GABA(A)/benzodiazepine receptor ligands. The goal of the work presented here is the identification and assessment of the minimum requirements for ligand recognition of GABA(A)/benzodiazepine receptors leading to activity at the sedation endpoint embedded in a common 3D pharmacophore for recognition. Using the experimental results, together with a systematic computational procedure developed in our laboratory, a five-component 3D pharmacophore for recognition of the GABA(A) receptor subtypes associated with the sedative behavioral response has been developed consisting of: two proton-accepting moieties, a hydrophobic region, a ring with polar moieties and an aromatic ring in a common geometric arrangement in all ligands having an effect at the sedation endpoint. To provide further evidence that the 3D pharmacophore developed embodied common requirements for receptor recognition, a pharmacophore analysis was performed for agonists, inverse agonists and antagonists separately. Each of the resulting pharmacophores contained the same five moieties at comparable distances to those found for the pharmacophore generated using all of them together. This result confirms that this pharmacophore constitutes a recognition pharmacophore representing required features in the overlapping portion of their binding sites. The reliability of this 3D pharmacophore was then assessed in several ways. First, it was determined that ligands that had no effect at the sedation endpoint did not comply with the pharmacophore requirements. Second, four benzodiazepine receptor ligands known to have an effect at the sedation endpoint, but not used in the pharmacophore development were found to satisfy the requirements of this pharmacophore. Third, the geometric and chemical requirements embedded in this pharmacophore were used to search 3D databases resulting in the identification of benzodiazepine receptor ligands known to affect sedation, but not included in the pharmacophore development. Finally, a 3D-quantitative structure analysis procedure (QSAR) model was developed based upon the ligands in the training set superimposed at their sedation pharmacophore points. The 3D-QSAR model shows good predictivity for binding of these ligands to receptor subtypes containing alpha1 but not alpha5 subunits. The pharmacophore developed for the sedation endpoint thus provides a predictive binding model for diverse ligand binding to alpha1 containing receptor subtypes.

Benzodiazepine-induced hyperphagia: development and assessment of a 3D pharmacophore by computational methods.

Benzodiazepine receptor (BDZR) ligands are structurally diverse compounds that bind to specific binding sites on GABA(A) receptors and allosterically modulate the effect of GABA on chloride ion flux. The binding of BDZR ligands to this receptor system results in activity at multiple behavioral endpoints, including anxiolytic, sedative, anticonvulsant, and hyperphagic effects. In the work presented here, a computational procedure developed in our laboratory has been used to obtain a 3D pharmacophore for ligand recognition of the GABA(A)/BDZRs initiating the hyperphagic response. To accomplish this goal, 17 structurally diverse compounds, previously assessed in our laboratory for activity at the hyperphagic endpoint, were used. The result is a four-component 3D pharmacophore. It consists of two proton acceptor atoms, the centroid of an aromatic ring and the centroid of a hydrophobic moiety in a common geometric arrangement in all compounds with activity at this endpoint. This 3D pharmacophore was then assessed and successfully validated using three different tests. First, two BDZR ligands, which were included as negative controls in the set of seventeen compounds used for the pharmacophore development, did not fit the pharmacophore. Second, some benzodiazepine ligands known to have activity at the hyperphagia endpoint, but not included in the pharmacophore development, were used as positive controls and were found to fit the pharmacophore. Finally, using the 3D pharmacophore developed in the present work to search 3D databases, over 50 classical benzodiazepines were found. Among them, were benzodiazepine ligands known to have an effect at the hyperphagic endpoint. In addition, the novel compounds also found in this search are promising therapeutic agents that could beneficially affect feeding behavior.

3D modeling, ligand binding and activation studies of the cloned mouse delta, mu; and kappa opioid receptors.

Refined 3D models of the transmembrane domains of the cloned delta, mu and kappa opioid receptors belonging to the superfamily of G-protein coupled receptors (GPCRs) were constructed from a multiple sequence alignment using the alpha carbon template of rhodopsin recently reported. Other key steps in the procedure were relaxation of the 3D helix bundle by unconstrained energy optimization and assessment of the stability of the structure by performing unconstrained molecular dynamics simulations of the energy optimized structure. The results were stable ligand-free models of the TM domains of the three opioid receptors. The ligand-free delta receptor was then used to develop a systematic and reliable procedure to identify and assess putative binding sites that would be suitable for similar investigation of the other two receptors and GPCRs in general. To this end, a non-selective, 'universal' antagonist, naltrexone, and agonist, etorphine, were used as probes. These ligands were first docked in all sites of the model delta opioid receptor which were sterically accessible and to which the protonated amine of the ligands could be anchored to a complementary proton-accepting residue. Using these criteria, nine ligand-receptor complexes with different binding pockets were identified and refined by energy minimization. The properties of all these possible ligand-substrate complexes were then examined for consistency with known experimental results of mutations in both opioid and other GPCRs. Using this procedure, the lowest energy agonist-receptor and antagonist-receptor complexes consistent with these experimental results were identified. These complexes were then used to probe the mechanism of receptor activation by identifying differences in receptor conformation between the agonist and the antagonist complex during unconstrained dynamics simulation. The results lent support to a possible activation mechanism of the mouse delta opioid receptor similar to that recently proposed for several other GPCRs. They also allowed the selection of candidate sites for future mutagenesis experiments.

Homology modeling and substrate binding study of human CYP2C9 enzyme.

The CYP2C subfamily of human liver P450 isozymes is of major importance in drug metabolism. The most abundant 2C isozyme, CYP2C9, regioselectively hydroxylates a wide variety of substrates. A major obstacle to understanding this specificity in human CYP2C9 is the absence of a 3D structure. A 3D model of CYP2C9 was built, assessed, and used to characterize explicit enzyme-substrate complexes using methods previously developed in our laboratory. The 3D model was assessed by determining its stability to unconstrained molecular dynamics and by comparison of specific properties with those of known protein structures. The CYP2C9 model was then used to characterize explicit enzyme complexes with three structurally and chemically diverse substrates: (S)-naproxen, phenytoin, and progesterone. Each substrate was found to bind to the enzyme with a favorable interaction energy and to remain in the binding site during unconstrained molecular dynamics. Moreover, the mode of binding of each substrate led to calculated preferred hydroxylation sites consistent with experiment. Binding-site residues identified for the models included Arg 105 and Arg97 as key cationic residues, as well as Asn 202, Asp 293, Pro 101, Leu 102, Gly 296, and Phe 476. Site-specific mutations are proposed for further integrated computational and experimental study.

Homology modeling and substrate binding study of human CYP2C18 and CYP2C19 enzymes.

It is well established that the variable binding-site architecture and composition of the P450 metabolizing heme proteins are major modulators of substrate and product specificity. Even the three closely related human liver isozymes, CYP2C9, CYP2C18, and CYP2C19, do not share all substrates and do not always lead to the same preferred hydroxylation products. The lack of knowledge of their three-dimensional (3D) structures has hindered efforts to understand the differences in their specificities. Building on previous work for the CYP2C9 enzyme, 3D models of CYP2C18 and 2C19 have been constructed and validated by computational methods developed and tested in our laboratory. They were used to characterize explicit enzyme-substrate complexes using the isoform-specific substrates progesterone and (S)-mephenytoin for 2C19 and 2-[2,3-dichloro-4-(3-hydroxypropyloxy)benzoyl]thiophene for 2C18. The results allowed both common and unique binding-site residues to be identified in each model. The calculated preferred hydroxylation site was obtained for each substrate and was found to be consistent with experimental observation. Comparisons were made among the 2C9, 2C18, and 2C19 model binding sites to investigate the subtle differences among them. These models can be used as structure-based guides for mutagenesis studies and screening of potential pharmaceuticals or toxins.

Molecular dynamics simulations of P450 BM3--examination of substrate-induced conformational change.

Cytochrome P450 BM3, of bacterial origin, is one of only five isozymes of the ubiquitous family of over 400 metabolizing heme proteins with a known crystal structure and only one of two with both substrate-free and substrate-bound forms determined. P450 BM3 is of particular interest since it has a similar function and similar substrates as mammalian P450s particularly of the 4A subfamily. Thus, the extent to which the substrate-free form of P450 BM3 undergoes a conformational change upon binding of a typical fatty acid substrate, palmitoleic acid, has been the subject of recent active experimental effort. Surprisingly, direct examination of the substrate-free (pdb2hpd.ent and pdb2bmh.ent) and substrate-bound (pdb1fag.ent) forms do not provide a clear answer to this question. The main reason for this ambiguity is that the two substrate-free monomers reported in the crystal structures themselves have significantly different conformations from each other, one with a more open substrate-access channel than the other. Since there is no way to tell to which substrate-free form the substrate binds, the effect of substrate binding cannot be deduced directly from comparisons of the experimental substrate-bound and substrate-free forms. The computational studies reported here have been designed to more robustly establish the effect of substrate binding on this isozyme. Specifically, molecular dynamics simulations were performed for each of the two substrate-free forms found in the asymmetric unit of the X-ray structure and for the two corresponding substrate-bound forms, constructed by docking palmitloeic acid into each of them. Comparisons of the results showed that palmitoleic acid binding had little effect on the conformation of the more closed substrate-free form of P450 BM3. By contrast, in the more open substrate-free form, this same substrate induced a closing of the entrance to the substrate-binding channel. The MD averaged structure of these two complexes obtained from docking of pamitoleic acid into the two asymmetric units of the substrate-free form were also compared to that obtained starting with the X-ray structure of the substrate-bound form. These results taken together led to the conclusion that, if indeed the substrate induces conformational changes in P450 BM3, the mouth of the substrate-access channel first closes down in response to the presence of the substrate, followed by rotation of the F-G domain to further optimize the P450 BM3-substrate interaction that would occur at a later stage.

Homology modeling and substrate binding study of human CYP4A11 enzyme.

Although both bacterial CYP102 (P450BM3) and mammalian CYP4A isozymes share a common function as fatty acid hydroxylases, distinctly different preferred sites of oxidation are observed with the CYP102 performing the usual non-terminal hydroxylation or epoxidation and the CYP4A enzymes performing the unusual and enigmatic terminal hydroxylation. The origin of this unique product specificity in human CYP4A11 has been explored in this work, focusing on possible differences in the binding site architecture of the two isozymes as the cause. To this end, 3D model structures of the human CYP4A11 enzyme were built and compared to the X-ray structure of CYP102. The substrate-binding channel identified in CYP4A11 was found to have a much more sterically restricted active site than that in CYP102 that could cause limited access of long-chain fatty acid to the ferryl oxygen leading to the preferred omega-hydroxylation. Results of docking of a common substrate, lauric acid, into the binding site of both CYP4A11 and CYP102 and molecular dynamics simulations provided additional support for this hypothesis. Specifically, in the CYP4A11-lauric acid simulations, the omega hydrogens were closest to the ferryl oxygen most of the time. By contrast, in the CYP102-lauric acid complex, the substrate could penetrate further into the active site providing access of the non-terminal (omega-1, omega-2) positions to the ferryl oxygen. These results, taken together, have elucidated the origin of the unusual product specificity of CYP4A11 and illustrated the central role of binding site architecture in subtle modulation of function.

Construction of a 3D model of cytochrome P450 2B4.

A three-dimensional structural model of rabbit phenobarbital-inducible cytochrome P450 2B4 (LM2) was constructed by homology modeling techniques previously developed for building and evaluating a 3D model of the cytochrome P450choP isozyme. Four templates with known crystal structures including cytochrome P450cam, terp, BM-3 and eryF were used in multiple sequence alignments and construction of the cytochrome P450 2B4 coordinates. The model was evaluated for its overall quality using available protein analysis programs and found to be satisfactory. The model structure was stable at room temperature during a 140 ps unconstrained full protein molecular dynamics simulation. A putative substrate access channel and binding site were identified. Two different substrates, benzphetamine and androstenedione, that are metabolized by cytochrome P450 2B4 with pronounced product specificity were docked into the putative binding site. Two orientations were found for each substrate that could lead to the observed preferred products. Using a geometric fit method three regions on the surface of the model cytochrome P450 structure were identified as possible sites for interaction with cytochrome b5, a redox partner of P450 2B4. Residues that may interact with the substrates and with cytochrome b5 have been identified and mutagenesis studies are currently in progress.

Thermodynamics of ligand binding to the cloned delta-opioid receptor.

The goal of this study was to determine the relative contribution of entropy and enthalpy to the free energies of binding to recombinant mouse delta-opioid receptors for the peptide agonist, DPDPE ([D-Pen2,D-Pen5]enkephalin), the peptide antagonist, TIPP(psi) (Tyr-Tic(psi)[CH2NH]Phe-Phe-OH), the nonpeptide agonist, SNC80 ((+)-4-[(alphaR)-alpha-((2S,5R)-4-allyl-2,5-dimethyl-1-piperazinyl )-3-methoxybenzyl]-N,N-diethylbenzamide), and the nonpeptide antagonist, naltrindole. Competitive binding studies were carried out using [3H]naltrindole at 0 degrees C, 12 degrees C, 25 degrees C and 37 degrees C, the affinities calculated and van't Hoff plots constructed for each ligand. The temperature dependence of binding and van't Hoff plots indicated that the entropy contribution is the major component of the free energy, for all four ligands, independent of its activity or chemical nature.

Refinement of 3D models of horseradish peroxidase isoenzyme C: predictions of 2D NMR assignments and substrate binding sites.

In this study, two alternative three-dimensional (3D) models of horseradish peroxidase (HRP-C)-differing mainly in the structure of a long untemplated insertion-were refined, systematically assessed, and used to make predictions that can both guide and be tested by future experimental studies. A key first step in the model-building process was a procedure for multiple sequence alignment based on structurally conserved regions and key conserved residues, including those side chains providing ligands to the two Ca2+ binding sites. The model refinements reported here include (1) optimization of side-chain conformations; (3) addition of structural waters using a template-independent procedure; (2) structural refinement of the untemplated 34 amino acid insertion located between the F and G helices, using both energy criteria and NMR data; (4) unconstrained energy optimization of the refined models. Using these procedures, two refined structures of HRP-C were obtained, differing mainly in the conformation of this long insertion. The presence of residues in this insertion that could potentially interact with bound substrates suggests a functional role that may be related to the general ability of class III peroxidases to form stable 1:1 complexes with a variety of substrates. The structural validity of the models was systematically assessed by a variety of criteria. Most notably, the ProsaII z scores and Profiles 3D scores of the two HRP-C models indicated that they are significantly better than would be obtained by simple amino acid replacement, using any of the known structures as a template. These two 3D HRP-C models, were then used to predict candidate residues for the assignment of NOESY cross-peaks previously noted in 2D-NMR studies. Specifically, the residues known as Ile X, Phe A, Phe B, aliphatic residue Q, and Ile T. Candidate substrate binding sites were also identified and compared with experimentally based predictions. This work is timely because new X-ray structures are anticipated that will facilitate the validation of these procedures.

Construction and evaluation of a three-dimensional structure of cytochrome P450choP enzyme (CYP105C1).

The purpose of this work was to develop and carefully evaluate improved strategies for constructing reliable 3-D models of P450 isozymes. To this end, a unique combination of steps for building and evaluating a model structure was used to build a homology model of the P450choP isozyme, based on knowledge of the X-ray structures of P450cam, P450terp, P450BM-3 and P450eryF. Specifically, the reliability of this model was examined by systematic comparisons of its conformational, energetic, environmental and packing properties and those of the four reference proteins with corresponding properties from the database of proteins with known structures. The results showed that the examined properties of this model structure are well within the criteria established for reliable structures and are of nearly as good quality as those of the reference proteins. In addition, the result from a 120 ps unconstrained MD simulation of the model with structural waters provided evidence that the model is stable at room temperature. This 3-D model can now be reliably used for explicit characterization of substrate and inhibitor complexes. Most importantly, although it is envisioned that building models for mammalian P450s will be even more challenging, the steps described here should be very useful in future construction of 3-D models of mammalian P450 isozymes.